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EXECUTIVE  SUMMARY 


The  ability  to  estimate  the  gas  capacity  of  carbon-based  filters  in  industrial  and 
military  applications  is  important  to  the  engineering  and  design  communities.  These  design 
aspects  factor  into  the  configuration,  performance  characteristics,  and  service  aspects  of  the 
filter.  A  common  goal  of  these  efforts  is  to  be  able  to  quantify  and  track  the  transmission  of  a 
compound  (typically  a  contaminant)  as  it  traverses  the  adsorbent  bed  from  the  inlet  to  the 
effluent. 


Within  the  field  of  adsorption  modeling  and  science,  extensive  efforts  have  been 
undertaken  and  continue  to  be  pursued  in  developing  analytical  tools  to  enable  the  engineering 
community  a  means  to  gauge  filter  performance.  The  analytical  models  are  founded  on  the  basic 
engineering  principles  of  mass,  energy,  and  momentum  conservation.  The  essential  element  of 
these  approaches  is  to  establish  a  fundamental  conservation  relationship  over  a  differential 
volume  element  and  then  integrate  the  resultant  differential  expressions  over  the  entire  volume  of 
the  reactor  vessel  to  characterize  the  performance  of  the  adsorber.  Although  the  development  of 
the  basic  differential  conservation  relationships  is  relatively  straightforward,  the  development  of 
closed- form  solutions  to  these  resultant  conservation  expressions  can  be  quite  formidable. 
Advances  in  modem  analytical  capability  enable  numeric  solutions  to  be  obtained;  however, 
closed-form  analytical  solutions  for  even  relatively  simplified  adsorption  processes  can  be  quite 
complex.  Because  of  these  substantial  analytical  obstacles,  engineering  and  design  engineers 
must  often  rely  on  correlative  approaches  to  characterize  filter  performance. 

The  following  analytical  model  has  been  used  (and  continues  to  be  used)  to 
estimate  the  breakthrough  capacity  of  carbon-based  filters  against  low-volatility  compounds  of 
military  interest: 


tb  =  [(We  A  p b)/(Q  Co)]  (L  -  AT)  (ES-1) 

where 

AT  =  { (cp,  dp)/[  6  ( 1  -  8)]  }  [(cp,  dp  u  p)/p]a41  [p/(p  B)f67  In  (CJCb)  (ES-2) 


The  development  of  these  relationships  seems  to  trace  back  to  efforts  conducted 
by  G.  S.  Bohart  and  E.  Q.  Adams  during  the  1920s  through  their  investigations  with  the 
transmission  of  chlorine  through  carbon  beds.  These  relationships  are  based  on  the  assumption 
that  bulk  diffusion  from  the  vapor  phase  to  the  adsorbed  phase  represents  the  rate -controlling 
removal  process. 

The  purpose  of  the  current  effort  is  to  provide  a  developmental  basis  and  rationale 
for  the  relationships  given  by  eqs  ES-1  and  ES-2.  The  approach  taken  in  this  review  effort  was  to 
demonstrate  the  complexity  of  the  basic  mass  balance  relationships,  the  need  to  employ 
simplifying  assumptions,  and  the  use  of  accepted  correlative  relationships  to  develop  the 
relationships  given  by  eqs  ES-1  and  ES-2.  This  review  is  intentionally  very  detailed  in  the 
mathematical  developments.  A  further  intent  of  this  effort  was  to  provide  a  historical  lineage  of 
the  currently  used  breakthrough  relationships  to  the  analytical  developments  of  the  past. 


iii 


During  the  conduct  of  this  review  effort,  there  were  analytical  developments  that 
were  not  thoroughly  understood  by  this  author.  These  knowledge  gaps  are  attributable  in  part  to 
mathematical  as  well  as  experience  limitations.  The  areas  that  were  not  well  understood  included 
the  detailed  form  of  eq  ES-3  (eq  195  within  the  report)  and  the  analytical  development  and 
progression  that  followed  eq  ES-4  (eq  232  within  the  report).  It  was  surmised  that  a  pressure 
factor  was  inadvertently  omitted  from  the  Dole  and  Klotz  relationship  given  by  ES-3. 
Dimensional  consistency  arguments  and  an  analytical  development  were  conducted  to  support 
this  contention. 


l(Po  Pim.ii)/ApJ  =  In (C0ICb) 


ES-3 


tb  =  [(pi,  we)/(u  C0  (MW),)]  {l 


Cb  L 

/[I  -  (Wi/We)]  (dCi/dzy'  d Ci  -  / (Wo/We)  dz. 
Co  0 
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+  [( u  Co  (MW)i)/(p,,  We)]  f  [Ci(L)/Co\  ( dCi/dty 1  dC,} 
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The  energetic  aspects  and  prior  loading  history  (or  service  environment)  of  the 
adsorbent  are  important  considerations  with  adsorption  processes.  Each  of  these  aspects  would 
be  expected  to  contribute  to  the  characterization  of  the  transmission  of  a  contaminant  through  the 
carbon  bed.  Because  of  the  added  analytical  complexity  associated  with  the  energy  balances  and 
prior  loading  considerations,  this  review  only  focused  on  fundamental  mass  balance 
requirements.  Energy  considerations  would  likely  have  to  be  considered  with  design  efforts 
requiring  a  high  degree  of  refinement.  These  issues  were  beyond  the  scope  of  this  review. 
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The  work  described  in  this  report  was  started  in  May  1995  and  completed 
in  September  2011. 

The  use  of  either  trade  or  manufacturers'  names  in  this  report  does  not 
constitute  an  official  endorsement  of  any  commercial  products.  This  report  may  not  be 
cited  for  purposes  of  advertisement. 

This  report  has  been  approved  for  public  release. 
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ANALYTICAL  ASPECTS  RELATING  TO  THE  ESTIMATION 
OF  CARBON  FILTER  PERFORMANCE  FOR  MILITARY  APPLICATIONS 


1.  BACKGROUND 

The  estimation  of  filtration  capabilities  for  carbon-based  military  filters  is 
typically  based  on  the  characteristics  of  the  compound  that  is  targeted  for  removal  from  a 
contaminated  airstream.  An  important  metric  relating  to  the  removal  of  a  compound  from 
a  contaminated  airstream  is  the  room-temperature  vapor  pressure  of  the  material.  In  many 
instances,  compounds  with  room-temperature  vapor  pressures  below  10  torr  are 
efficiently  removed  on  activated  carbonaceous  adsorbents  through  a  process  known  as 
physical  adsorption.  This  grouping  of  materials  covers  the  known  H,  G,  and  V  classes  of 
chemical  warfare  compounds.  With  vapor  pressures  greater  than  approximately  100  torr, 
impregnants  are  incorporated  into  the  pore  structure  of  the  activated  carbon,  and  the 
removal  mechanisms  are  generically  classified  as  chemisorption. 

This  report  will  review  the  analytical  foundations  upon  which  the  carbon 
filter  performance  for  military  applications  is  based.  As  indicated  earlier,  considerable 
effort  has  already  been  put  forth  in  this  area,  and  to  that  end,  the  content  presented  here 
represents  more  of  a  compilation  of  past  efforts  rather  than  anything  new  or  novel.  The 
review  does  expand  upon  mathematical  developments,  to  provide  sufficient  background 
so  as  to  be  understood  by  the  design  community. 

The  foundation  upon  which  the  performance  estimations  are  based  rests 
with  fundamental  mass  balance  relationships  that  are  typically  developed  in  differential 
form.  Solution  to  the  resultant  differential  mass  balances  would  yield  vapor-phase 
concentrations  as  a  function  of  both  time  and  position  within  an  adsorbent  bed.  As  is 
demonstrated,  closed-form  analytical  solutions  to  the  resultant  differential  mass  balance 
relationships  are  complex;  therefore,  simplifying  assumptions  are  necessary.  Even  when 
the  simplifying  assumptions  are  invoked,  analytical  solutions  to  the  differential  mass 
balance  expressions  can  be  formidable. 


2.  ANALYTICAL  CONSIDERATIONS  OF  DIFFERENTIAL  VOLUME 

ELEMENTS  AND  DIRECTIONAL  GRADIENTS 

The  intent  of  this  section  is  to  consider  differential  volume  elements  and 
the  flow  fields  associated  with  such  elements  for  varying  geometries.  For  rectilinear 
geometries,  an  orthogonal  volume  element  can  often  be  used,  whereas  for  radial  flow 
geometries,  a  volume  element  having  cylindrical  characteristics  would  be  appropriate. 
The  M48A1  and  M98  gas  filters  used  with  collective  protection  applications  are 
examples  of  radial  flow  filters. 
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A  note  relating  to  the  nomenclature  convention  used  throughout  this 
review  is  needed  at  this  stage.  The  following  nomenclature  terminology  and 
interpretation  were  used  throughout  this  report: 

•  Nix(x)  represents  the  molar  flux  of  component  i  in  the  v-direction 
evaluated  at  position  x. 

•  Ntx(x  +  Ax)  represents  the  molar  flux  of  component  i  in  the  v-direction 
evaluated  at  position  (x  +  Ax). 

•  Nj(t,  r,  0,  z)  represents  the  molar  flux  of  component  i  evaluated  at 
time  t  and  coordinate  position  (r,  0,  z). 

•  Nj(t  +  At,  r  +  A r,  0  +  A0,  z  +  A z)  represents  the  molar  flux  of 
component  i  evaluated  at  time  (t  +  At)  and  coordinate  position 
(r  +  A r,  0  +  A0,  z  +  Az). 

This  nomenclature  convention  was  chosen  in  an  effort  to  explicitly  define 
the  variable  dependency  as  well  as  the  evaluation  conditions  of  the  affected  parameters. 
In  cases  where  the  functional  dependency  and/or  evaluation  conditions  are  unambiguous, 
this  explicit  descriptor  notation  is  omitted. 

The  development  of  the  mass  balance  relationships  in  this  review  effort  is 
based  on  the  characterization  of  differential  volume  elements  with  respect  to  a  fixed 
reference  frame.  This  is  in  contrast  with  the  characterization  of  the  volume  elements 
relative  to  a  reference  frame  that  undergoes  positional  displacements. 


2.1  Differential  Volume  Element  Considerations  with  Three  Mutually 

Perpendicular  Axes 

Figure  1  depicts  an  orthogonal  differential  volume  element  in  free  space 
having  three  mutually  perpendicular  axes. 
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Niy(y  +  Ay) 


Figure  1 .  Three  mutually  perpendicular  molar  fluxes 
in  a  Cartesian  coordinate  system. 


The  volume  of  the  differential  element  depicted  in  Figure  1  is  Ax  Ay  A z. 

A  generic  mass  balance  for  component  i  based  on  a  reacting  system  can  be 

written  as 

(dMi/dt)  =  2(777,) in  -  2(771,) out  +  (MW),  R,  Ax  Ay  A z  (1) 

where  R,  represents  the  molar  rate  of  generation  (or  consumption)  of  component  i  within 
the  volume  element  expressed  on  a  per-volume  basis.  The  sum  of  the  mass  flows  of 
component  i  into  the  differential  volume  element  can  be  expressed  by 

2(777, )in  =  (MW);  Nix(x)  Ay  Az  +  (MW),-  Niy(y)  Ax  A z  +  (MW),-  Niz(z)  Ax  Ay 

(2) 

=  (MW),-  [Ay  Az  Nix(x)  +  Ax  A z  Niy(y )  +  Ax  Ay  Niz(z )] 

The  sum  of  the  mass  flows  of  component  i  exiting  the  differential  volume 
element  can  be  expressed  by 

2(777, )out  =  (MW),-  Nix(x  +  Ax)  Ay  A z  +  (MW),-  Niy(y  +  Ay)  Ax  A z 

+  (MW),  Niz(z  +  Az)  Ax  Ay 

(^) 

=  (MW),-  [Ay  Az  Nix(x  +  Ax)  +  Ax  A z  Niy(y  +  Ay) 

+  Ax  Ay  Niz(z  +  Az)] 
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Substituting  the  relationships  given  by  eqs  2  and  3  into  eq  1  gives  the 
following  expression  after  some  simplification: 


(dMi/dt)  =  -(MW);  {  Ay  Az  \Nix(x  +  Ax)  -  /V,x(*)]  +  Ax  Az  [Niy(y  +  Ay)  -  Niy(y)] 

+  Ax  Ay  [Niz(z  +  Az)  -  Niz(z )]  }  +  (MW),  R,  Ax  Ay  A z 


(4) 


The  vapor  phase  concentration  of  component  i  within  the  differential 
volume  element  is  taken  to  be  homogeneous  and  characterized  by  a  uniform  molar 
concentration  of  Q.  The  mass  of  component  i  within  the  differential  volume  element  is 
then  given  by 


Mi  =  (MW),-  Q  Ax  Ay  A z 


(5) 


Equation  5  can  be  differentiated  with  respect  to  time,  which  yields 

(dMj/dt)  =  (MW),-  C;  [o(Ax  Ay  A z)/dt]  +  Ax  Ay  Az  {  c[(MW),  C,]/S?} 

=  (MW);  C;  [d(Ax  Ay  A z)/dt]  +  Av  Ay  Az  (MW),  (dQ/dt) 


The  full  expansion  of  eq  6  will  yield  terms  of  the  form 

Ay  Az  [o(Av')/o/J,  At  Az  [5(A y)/dt\,  and  Ar  Ay  [c(Az)/o/] 

In  the  following  development,  it  is  assumed  that  the  derivatives  of  the 
differentials  are  negligible  relative  to  the  first-order  derivatives.  Mathematically,  this 
requires  that  the  following  holds  true: 

[3(A x)ldt\,  [S(Av)/g7],  and  [o(Az)/o/J  «  ( dCJdt ) 

With  this  assumption,  eq  6  is  approximated  by 

(dMi/dt)  ~  (MW);  Ar  Ay  Az  (dQ/dt)  (7) 

Substituting  the  relationship  given  by  eq  7  into  eq  4  and  dividing  through 
by  the  quantity  [Ax  Ay  Az  (MW),]  gives 

(dQ/dt)  =  -  {  [Nix (x  +  Ax)  -  Nix(x)]/Ax  }  -  {  [Niy(y  +  Ay)  -  Niy(y)]/Ay  } 

(8) 

-  {  [Niz(z  +  Az)  -  Niz(z)\/Az  }  +Ri 

In  the  limit  as  Ax,  Ay,  and  Az  collapse  to  zero,  eq  8  takes  on  the  form 

(dCi/dt)  =  ~  [( dNJdx )  +  (dNiy/dy)  +  ( dNjJdz )]  +  Ri  (9) 


4 


Each  of  the  three  orthogonally  directed  flux  terms  given  in  eq  9  is 
comprised  of  a  bulk  flow  component  as  well  as  a  dispersive  (or  diffusive)-based 
contribution.  The  dispersive-based  contribution  is  typically  based  on  a  Fickian  diffusion- 
type  relationship  (i.e.,  Fick’s  first  law  of  diffusion).  The  dispersive  contribution  attempts 
to  capture  deviations  from  ideal  plug-based  flow.  The  following  generic  relationship 
characterizes  the  molar  flux  through  each  of  the  orthogonally  directed  axes: 


(Nj y)total 


=m  convective  bulk  flow  “1“  (-^ (^diffusion  contribution 

=  uj  Q  -  D ij  (dCj/dj) 


(10) 


Differentiation  of  eq  10  with  respect  to  a  generic  directional  axis  gives  the 
following  series  of  relationships: 


mm = [8(uj  c,)/dj\  -  {s[By  (dQ/em/dj)  ai) 

If  it  can  be  assumed  that  the  directionally  oriented  velocities  are  constant, 
then  eq  1 1  can  be  expressed  as 

(dNij/dj)  =  Uj  [ dCi/dj ]  -  {  8[Bij  (dQ/8j)]/8j  }  (12) 

Equation  12  can  be  further  simplified  if  it  can  be  assumed  that  a  constant- 
dispersion  coefficient  can  be  used  to  characterize  the  molar  flux  of  i  in  each  of  the  three 
orthogonal  directions.  The  directional  flux  gradients  can  then  be  expressed  as  follows  for 
each  of  the  three  coordinate  axes: 


(8NjX/8x)  =  ux  (8Cj/8x)  -  B/v  ( d2Ci/dx 2)  (13) 

(8Niy/8y)  =  uy  (dQ/dy)  -  B(V  (S2C,%2)  (14) 

(8Nlz/8z)  =  uz  ( dQ/dz )  -  B/z  ( 82Q/8z 2)  (15) 


Substituting  the  relationships  given  by  eqs  13-15  into  eq  9  gives 

(i 8Ci/8t )  +  ux  (8Ci/8x)  +  uy  ( 8Ci/8y )  +  uz  ( 8Ci/8z ) 

=  B ix  ( 82Ct/8x 2)  +  B iy  ( 82Ci/8y 2)  +  Biz  (, d2Q/8z 2)  +  R, 


(16) 


Equation  16  can  be  further  simplified  by  assuming  that  the  directional 
dispersion  coefficients  are  equal  to  one  another.  If  a  uniform,  nondirectional  dispersion 
coefficient  depicted  by  Bl  can  be  used  to  characterize  the  transport  of  i  within  the 
system,  eq  16  can  be  expressed  as 
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(17) 


(i dCi/dt )  +  ux  ( dCi/dx )  +  uy  ( dCi/dy )  +  uz  ( dCi/dz ) 

=  [(d^/ar)  +  (S2C,/dy2)  +  (< 82Q/dz 2)]  +  /?,- 

Equation  17  represents  the  differential  mass  balance  for  component  i  for  a 
system  that  is  devoid  of  an  adsorbent  packing.  The  expression  was  developed  from  a 
simple  geometry  with  inherent  assumptions.  One  can  observe  that  the  relationship  is 
rather  complex  if  the  gradients  are  retained  for  each  of  the  three  orthogonal  directions.  It 
is  often  assumed  that  the  superficial  flow  velocities  in  the  non-principal  flow  directions 
can  be  considered  as  negligible  relative  to  the  primary  flow  direction.  Taking  the 
principle  flow  direction  to  be  along  the  "-axis  and  assuming  that  the  concentration 
gradients  in  the  x-  and  y-directions  are  negligible,  eq  17  simplifies  to 

(dCi/dt)  +  uz  ( dCi/dz )  =  Bl  ( d2Ci/dz 2)  +  Ri  (18) 

Although  eq  18  represents  a  simplified  form  of  the  differential  mass 
balance,  it  is  still  quite  formidable  for  many  practical  applications.  The  following  section 
addresses  the  analytical  development  for  a  differential  volume  element  having  cylindrical 
characteristics.  Many  of  the  same  assumptions  that  were  used  for  the  orthogonally  based 
coordinate  system  are  also  used  for  the  cylindrical  coordinate  system. 


2.2  Differential  Volume  Element  Considerations  for  a  Cylindrical 

Coordinate  System 

Figure  2  depicts  a  differential  volume  element  in  a  cylindrical  coordinate 

system. 


Figure  2.  The  coordinates  associated  with  a  cylindrical  system. 
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The  three  coordinate  directions  of  concern  in  a  cylindrical  coordinate 
system  are  the  radial  (r),  angular  (0),  and  axial  (z)  directions.  The  volume  of  the 
differential  element  that  is  considered  with  the  cylindrical  coordinate  system  is  given  by 


AVcyl  =  A3  Az  Ar 

(19) 

AVcyi  =  r  A0  Az  Ar 

(20) 

The  mass  of  component  i  within  the  differential  volume  element  is  given  by 


Mi  =  (MW),  Q  A  ycyi 


(21) 


Differentiating  eq  21  with  respect  to  time  gives 

(dMJdt)  =  (MW),  {  C,  [S(AVcyi )/dt\  +  AVcyl  (, dQ/dt ) }  (22) 

As  was  the  case  for  Cartesian  coordinate  system  development,  the 
derivative  of  the  differential  is  assumed  to  be  negligible  relative  to  the  first-order 
derivative,  such  that 


[S(AVcyi )/dt]  «  ( dCj/dt ) 

With  this  assumption,  eq  22  is  approximated  as 

(dMi/dt)  ~  (MW),-  AVcyi  (dQ/dt)  =  (MW),  r  A0  Az  Ar  (dQ/dt)  (23) 

where  the  relationship  given  by  eq  20  has  been  substituted  for  (AVcyi)  in  eq  23.  The  net 
mass  flows  associated  with  each  of  the  three  directions  are  then  considered  as  follows. 


2.2. 1  Net  Rate  of  Mass  Flow  of  Component  i  Associated  with  the 

Radial  Direction 

The  mass  flow  associated  with  the  radial  direction  is  taken  to  be  in  the 
outward  direction.  The  area  associated  with  this  flow  is  (r  A0)(Az).  The  cross-sectional 
area  associated  with  this  directional  flux  varies  with  radial  position.  The  following 
expression  describes  the  rate  of  mass  flow  of  component  i  entering  the  differential 
volume  element  at  position  r: 

mdr  =  (MW);  Nir(r)  [r(r)  A0  A z]  (24) 

An  expression  for  the  rate  of  mass  flow  of  component  i  exiting  the 
differential  volume  element  at  position  (r  +  Ar)  is  then  given  by 

(771;  W  =  (MW);  Nir(r  +  Ar)  [r(r  +  Ar)  A0  A z]  (25) 
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A  note  of  clarification  with  respect  to  the  nomenclature  convention  used  in 
eqs  24  and  25  and  throughout  this  review  effort  is  emphasized  here.  Through  the 
examples  that  follow,  this  nomenclature  approach  is  demonstrated.  The  following 
nomenclature  and  interpretation  was  used  throughout: 

•  Njr(r ):  represents  the  molar  flux  of  component  i  in  the  radial  direction 
evaluated  at  radial  position  r. 

•  Nir(j  +  A r):  represents  the  molar  flux  of  component  i  in  the  radial 
direction  evaluated  at  radial  position  (r  +  Ar). 

•  r(r):  represents  the  radial  position  from  a  central  reference  evaluated  at 
radial  position  r. 

•  r(r  +  Ar):  represents  the  radial  position  from  a  central  reference 
evaluated  at  radial  position  (r  +  Ar). 

•  Nu-(r)  r(r):  represents  the  product  of  the  molar  flux  of  component  i  in 
the  radial  direction  and  the  radius,  both  evaluated  at  radial  position  r. 

•  /V,>(r  +  Ar)  r(r  +  Ar):  represents  the  product  of  the  molar  flux  of 
component  i  in  the  radial  direction  and  the  radius,  both  evaluated  at 
radial  position  (r  +  Ar). 

The  net  rate  of  mass  flow  of  component  i  associated  with  the  r-direction  is 
the  difference  between  the  influent  and  effluent  rates  or  the  difference  between  eqs  24 
and  25.  The  result  is  given  by 


[(m,)dnet  =  (fUdr  ~  (HlfW 

=  -  (MW);  A9  Az  [r(r  +  Ar)  Nu{  r  +  Ar)  -  r(r)  Nir( r)  \ 


(26) 


2.2.2  Net  Rate  of  Mass  Flow  of  Component  i  Associated  with  the  0-Direction 

The  area  associated  with  this  angular  flow  is  (Ar  A z).  The  cross-sectional 
area  associated  with  this  directional  flux  is  invariant  with  respect  to  the  angular  position. 
The  following  expression  describes  the  rate  of  mass  flow  of  component  i  entering  the 
differential  volume  element  at  position  0: 


(Tflj){)  =  (MW);  Ar  Az  /V;o(0) 


(27) 


An  expression  for  the  rate  of  mass  flow  of  component  i  exiting  the 
differential  volume  element  at  position  (0  +  A9)  is  then  given  by 


(m,)0+Ae  =  (MW);  Ar  Az  A;0(0  +  A0)  (28) 
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The  net  rate  of  mass  flow  of  component  i  associated  with  the  0-direction  is 
the  difference  between  the  influent  and  effluent  rates  or  the  difference  between  eqs  27 
and  28.  The  result  is  given  by 


[(m,)eJnet  =  (m,)o  -  (m,)e+A0 


=  -(MW);  Ar  Az  [NdQ  +  A0)  -  /V,o(0)] 


(29) 


2.2.3  Net  Rate  of  Mass  Flow  of  Component  i  Associated  with  the  z-Direction 

The  area  associated  with  the  axially  directed  flow  is  (r  A0  Ar).  The  cross- 
sectional  flow  area  associated  with  this  directional  flux  does  not  vary  in  the  axial 
direction.  The  following  expression  describes  the  rate  of  mass  flow  of  component  i 
entering  the  differential  volume  element  at  position  z: 


(jn,i)z  -  (MW);  r  A0  A rNiz(z) 


(30) 


An  expression  for  the  rate  of  mass  flow  of  component  i  exiting  the 
differential  volume  element  at  position  (z  +  Az)  is  then  given  by 


(Vldz+Az  =  (MW),-  r  A0  Ar  Niz(z  +A z)  (3 1 ) 


The  net  rate  of  mass  flow  of  component  i  associated  with  the  z-direction  is  the  difference 
between  the  influent  and  effluent  rates  or  the  difference  between  eqs  30  and  31.  The 
result  is  given  by 


L(m().dnet  =  (m,)z  -  an,)z-Az 

=  -(MW),-  r  A0  Ar  [Niz(z  +Az)  -  Niz(z)] 


(32) 


2.2.4  Composite  Mass  Balance  for  Component  i:  Cylindrical  Volume  Element 

A  generic  mass  balance  for  component  i  based  on  the  cylindrical 
coordinate  system  depicted  in  Figure  2  for  a  reacting  system  can  be  expressed  in  an 
analogous  manner  to  that  for  the  orthogonal  coordinate  system.  The  generic  mass  balance 
is  given  by 

(oMJdt)  =  E(m()in  -  E(/tt()out  +  (MW),  R,  r  AO  Az  Ar  (33) 

As  expressions  for  the  directional  mass  flows  of  component  i  have  been 
developed,  eq  33  can  be  equivalently  expressed  through  the  following  relationship: 
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(dMi/dt)  =  2(771/) net  +  (MW);  R,  r  AO  Az  A r 

=  mdrU  +  [(m,)0]net  +  [(  771;)  Jnet  ( 

+  (MW),-  /?,  r  A9  Az  Ar 

Substituting  expressions  given  by  eqs  23,  26,  29,  and  32  into  the  mass 
balance  relationship  given  by  eq  34  yields 


(< dCi/dt )  =  -(l/r)  {  [r(r  +  Ar)  Njr(r  +  A r)  -  r(r)  Af,y(r)]/Ar  } 

-(1/r)  {  [iV,e(0  +  A0)  -  /V,o(0)]/A0  }  (35) 

-{[W(z  +Az)-W4z)]/Az]  +Ri 

In  the  limit  as  Ar,  A9,  and  A z  collapse  to  zero,  eq  35  takes  on  the  following  form: 

(dCi/dt)  +  (1/r)  [8(r  Nir)/8r]  +  (1/r)  (dNi()/dQ)  +  ( dNiz/8z )  =  Ri  (36) 

A  useful  comparison  of  eq  36  can  be  made  with  respect  to  a  similarly 
developed  relationship  appearing  in  a  standard  reference  text.  Equation  36  is  consistent 
with  the  relationship  given  in  Table  18.2-1  (eq  B)  of  the  textbook  Transport  Phenomena 
(ref  1). 


2.2.5  Relating  Molar  Flux  and  Directional  Flux  Gradients  to  Concentration  and 

Diffusivity  Parameters 

The  molar  flux  and  directional  flux  gradients  have  forms  consistent  with 
the  expressions  given  by  eq  11.  These  relationships  expressed  in  terms  of  a  cylindrical 
coordinate  system  are  as  follows: 


[d(r  Nir)/8r ]  =  [8{r  ur  Ct)/8r ]  -  D„.  {  8[r  ( 8Ci/dr)]/dr  }  (37) 

(1/r)  (8N ,-()/ 50)  =  (1/r)  [8(u()  Q)/8Q \  -  D/e  [(1/r2)  (52C/502)|  (38) 

and 

(8Niz/8z)  =  [8(uz  Ci)/8z]  -  D„  (52C,/oz2)  (39) 

The  relationship  for  the  angular  flux  gradient  (eq  38)  is  addressed  more 
fully  in  Appendix  A.  With  the  flux  gradient  relationships  given  by  eqs  37,  38,  and  39,  it 
has  been  assumed  that  the  directional  dispersion  coefficients  are  constant-valued. 
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Substitution  of  the  relationships  given  by  eqs  37,  38,  and  39  into  eq  36 
yields  the  following  expression: 

(dCj/dt)  +  [G(r  ur  C,)/5r]  +  (1/r)  [G(uq  C,)/50]  +  [G(uz  C,)/oz] 


(40) 

=  Ri  +  D„  { d[r  (dC,/dr)]/dr}  +  B(0  [(1/r2)  (52CyS02)]  +  Biz  (G2Ci/Gz2) 

If  the  directional  dispersion  coefficients  are  assumed  to  be  equal  to  one  another,  such  that 

B;>  =  B,-e  =  Biz  =  B/_ 


then  eq  40  simplifies  to 

(GQ/Gt)  +  [G(r  ur  C,)/Gr]  +  (1/r)  [G(uq  C,)/50]  +  [G(uz  C,)/ozJ 

=  Ri  +  Bl{  {  G[r  (GCj/dr)\/Gr  j  +  [(1/r2)  (o2C,/502)]  +  (< 82Q/dz 2)  } 


The  full  expansion  of  eq  41  gives  the  following  relationship: 

(< dCi/dt )  +  ur  ( GCi/dr )  +  (we/r)  (<9C(7o0)  +  uz  ( dCi/dz ) 

+  Ci  [( ur/r )  +  ( dujdr )  +  (1/r)  (due/60)  +  (duz/8z)\  (42) 

=  Ri  +  (B dr)  (GCGGr)  +  B/.  [(G2Ci/Gr)  +  (1/r2)  (52C,/502)  +  (o2Cl/oz2)\ 

With  the  assumption  of  a  constant  mass  density,  continuity  requires  that 
the  following  relationship  holds  (see  Appendix  B  for  development)  relative  to  a  fixed 
spatial  reference  frame: 

(ujr)  +  ( dujdr )  +  (1/r)  (d«e/50)  +  ( dujdz )  =  0  (43) 


With  the  relationship  given  by  eq  43,  eq  42  simplifies  to 


(dCi/dt)  +  ur  ( dCJdr )  +  (we/r)  (5C//50)  +  uz  ( dCJdz ) 

=  Ri  +  (B//r)  ( GQ/dr )  +  BL  [( d2Ci/dr 2)  +  (1/r2)  (S2C//d02)  +  (d2Cl/dz2)\ 


(44) 


Recognizing  that  the  following  relationship  holds: 


(B//r)  {  8[r  (dQ/dr)]/dr}  =  (B Jr)  [r  (d2Ci/dr2)  +  (GQ/dr)] 

=  Bl  (GrCJGr)  +  (B//r)  (GQ/dr) 


(45) 


Equation  44  can  be  rewritten  as 


11 


(dCi/dt)  +  ur  ( dCj/dr )  +  (we/r)  (oC,/o0)  +  w,  (dC/dz) 


(46) 


=  /?,-  +  DL  {  { (1/r)  {  8[r  {dCtldr)\ldr }  }  +  (1/r2)  (S2C,/c02)  +  (< 32C;/& 2)  } 

On  a  term-by-term  basis,  eq  46  is  consistent  with  eq  (B)  in  Table  18.2-2 
(page  559)  of  Transport  Phenomena  (ref  1).  The  analytical  complexity  of  eq  46  is  readily 
apparent.  The  relationship  was  developed  with  the  following  assumptions: 

•  the  directional  dispersion  coefficients  are  constant-valued  and  equal  to 
one  another, 

•  the  overall  mass  density  of  the  system  is  constant,  and 

•  the  reference  frame  of  the  differential  volume  element  is  spatially 
fixed  with  respect  to  its  coordinate  axes. 


2.3  Closing  Commentary 

The  analysis  conducted  on  the  differential  volume  elements  in  the 
Cartesian  and  cylindrical  coordinate  systems  was  performed  to  demonstrate  the  analytical 
complexity  that  arises  when  one  considers  the  directional  fluxes  associated  with  each 
coordinate  system.  The  preceding  analyses  were  also  simplified  in  that  the  volume 
elements  were  devoid  of  adsorbent  packing.  With  an  adsorbent  packing,  considerations 
relating  to  the  void  fraction  and  bulk  density  would  have  to  be  incorporated  into  the 
analytical  expressions.  The  reaction  rate  terms  appearing  in  the  relationships  would  also 
have  to  be  considered  in  a  rigorous  analytical  analysis  of  the  mass  balance  expressions. 
Depending  on  the  mechanistic  details  associated  with  the  chemical  reactionary  process, 
these  rate  expressions  would  add  analytical  complexity  to  the  relationships. 

A  goal  of  the  preceding  discussion  was  to  provide  justification  and 
rationalization  for  using  analytical  expressions  based  on  gradients  in  a  singular  direction. 
For  fixed  cross-sectional  area  geometries  in  which  an  orthogonal  coordinate  system  may 
be  appropriate,  the  principle  directional  gradient  would  be  in  the  direction  of  primary 
flow  consistent  with  an  orthogonally  oriented  axis.  The  gradients  and  velocities 
associated  with  the  non-principle  flow  axes  are  often  assumed  to  be  negligible  relative  to 
the  primary  flow  axis.  Taking  z  to  be  the  primary  flow  axis  with  a  Cartesian  coordinate 
system,  eq  18  would  typically  be  used: 

(dCj/dt)  +  uz  (dQ/dz)  =  BL  ( d2CJdz 2)  +  R,  (18) 

The  relationship  given  by  eq  18  would  also  apply  to  a  cylindrical 
coordinate  system  in  which  the  cross-sectional  flow  area  normal  to  the  principle  flow 
direction  was  constant  (i.e.,  axial  flow). 

Through  similar  reasoning  with  a  radial  flow  system  in  a  cylindrical 
coordinate  system,  gradients  and  velocities  in  the  angular  and  axial  directions  are  often 
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considered  to  be  negligible  relative  to  the  gradients  and  velocity  associated  with  the 
primary  or  radially  directed  flow.  Assuming  that  the  directional  gradients  associated  with 
the  angular  and  axial  directions  are  negligible  relative  to  the  gradients  in  the  radial 
direction,  eq  46  simplifies  to 

(dCi/dt)  +  ur  ( dCt/dr )  =  R,  +  (D Jr)  {  8[r  (dCl/Gr)]/8r  }  (47) 

Despite  these  simplifying  assumptions,  eq  47  retains  significant  analytical 
complexity.  A  further  simplification  to  eq  47  can  be  made  if  the  dispersive  effects 
associated  with  the  flow  of  i  can  be  neglected.  In  that  case,  BL  =  0,  and  eq  47  becomes 

(dCi/dt)  +  ur  (dCj/dr)  =  Rj  (48) 


3.  CONSIDERATIONS  ON  THE  PARAMETERS  USED  IN 

ADSORBENT-BASED  FILTRATION  SYSTEMS 

3.1  Void  Fraction 

The  void  fraction,  depicted  by  8,  is  simply  the  ratio  of  the  void  volume  to 
the  total  volume  of  an  adsorbent  packing.  The  total  volume  is  comprised  of  the  bulk 
material  volume  of  the  adsorbent  as  well  as  the  interstitial  volume  about  the  individual 
bulk  particles.  The  void  fraction  does  not  include  the  intraparticle  volume  associated  with 
the  pore  structure  within  the  individual  granular  particles.  The  intraparticle  volumes 
associated  with  the  porous  networked  structure  contribute  significantly  to  the 
characteristics  of  carbon-based  adsorption;  however,  the  void  fraction  is  more  of  a  bulk¬ 
packing  characteristic.  Some  useful  relationships  associated  with  the  void  fraction  are: 

8  =  void  volume/total  volume 


Vvoid/Vtotal 


(49) 


—  Vvoid/(Vads  +  Vvoid) 


Other  groupings  involving  the  void  fraction  that  factor  into  subsequent 
analytical  developments  are 


(Vads/Vtot)  =1-8  (50) 

and 

(VadS/Vvoid)  =  (l-£)/8  (51) 
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3.2 


Dense  Packing  of  Ideal  Spheres 


A  void  fraction  of  0.39  is  typically  assumed  for  the  adsorbent  loadings  of 
12  x  30  granular  activated  carbon.  For  very  high-efficiency,  vapor-phase  contaminant 
reductions,  dense  packing  of  the  adsorbent  is  required. 

Three  typical  lattice  arrangements  are  often  used  to  model  atomic  and 
molecular  structures.  These  include  the  body-centered-cubic  (BCC),  the  hexagonal-close- 
packed  (HCP),  and  the  face-centered-cubic  (FCC)  arrangements.  Of  the  three,  the  HCP 
and  FCC  arrangements  provide  a  more  dense  packing  than  the  BCC  arrangement.  The 
HCP  arrangement  can  be  likened  to  an  ABAB. . .  sequential  layering  of  ideal  spheres, 
whereas  the  FCC  arrangement  can  be  likened  to  an  ABCABC. . .  sequential  layering  of 
ideal  spheres,  in  which  additional  layerings  are  placed  in  the  interstitial  depressions 
created  by  the  preceding  coplanar  layering  of  spheres. 

In  the  following  analysis,  it  was  assumed  that  the  carbon  particles  were 
comprised  of  ideal,  uniform  spheres,  packed  in  an  FCC  arrangement. 


3.2. 1  Unit  Cell  Considerations 

Figure  3  provides  a  description  of  the  lattice  sites  that  were  considered  in 

this  analysis. 


Comer  Site  Location  (typical) 


Facial  Site  Location  (typical) 


Figure  3.  Lattice  site  depictions  of  an  FCC  unit  cell. 


The  solids  volume  associated  with  a  singular  unit  cell  of  an  ideal  FCC 
lattice  arrangement  (depicted  in  Figure  3)  is  determined  as  follows: 

Corner  spheres:  (8  comers/unit  cell)  x  [(1/8  sphere)/corner)] 

=  1  sphere/unit  cell 

Facial  spheres:  (6  faces/unit  cell)  x  [(1/2  sphere)/face] 

=  3  spheres/unit  cell 

The  net  number  of  spheres  associated  with  the  corner  sites  of  an  FCC  unit 
cell  is  1,  and  the  net  number  of  spheres  associated  with  the  facial  planes  of  an  FCC  unit 
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cell  is  3.  Therefore,  a  total  of  four  spheres  are  associated  with  the  unit  cell  of  an  ideal 
FCC  lattice  arrangement: 


Number  of  spheres  associated  with  FCC  unit  cell  =  4 

The  solids  volume  associated  with  a  unit  cell  is  therefore  4x  the  volume  of 
a  singular  sphere,  or 

Vsoiids  =  4  x  (4/3)  7i  r3  =  ( 1 6/3)  n  (d/2)3  =  (2/3)  n  d3  (52) 

The  volume  of  an  FCC  unit  cell  is  determined  through  consideration  of  a 
face  plane.  A  representative  face  plane  is  depicted  in  Figure  4. 


Figure  4.  Representative  face  plane  of  an  FCC  unit  cell  arrangement. 


The  spheres  located  along  the  diagonals  of  Figure  4  are  in  direct  contact  in 
an  FCC  unit  cell  arrangement.  These  diagonal  spheres  form  the  hypotenuse  of  a  right 
triangle  formed  by  the  sides  of  unit  cell.  With  this  direct  contact,  the  hypotenuse  of  the 
right  triangle  formed  has  a  length  of  2d.  The  sides  of  the  FCC  unit  cell  have  the 
dimension  given  by  x,  which  can  be  solved  in  terms  of  the  sphere  diameter  to  yield 

1/2  i 

x  =  2  a 

The  total  volume  of  the  FCC  unit  cell  can  be  expressed  in  terms  of  the 
diameter  of  the  sphere,  yielding 


Vfcc  unit  cell  =^3  =  23/2J3 


(53) 


3.2.2  Void  Fraction  Considerations 

In  this  section,  the  definition  of  a  void  fraction  (8)  is  used  in  conjunction 
with  the  relationships  developed  for  the  FCC  unit  cell  arrangement.  The  void  fraction,  as 
used  in  adsorbent  packings,  was  provided  by  eq  49,  and  can  be  expressed  as 

8  =  void  volume/total  volume 
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Because  the  summation  of  the  void  and  solids  volume  constitutes  the  total 
volume  of  the  unit  cell,  the  void  volume  can  be  expressed  as 


kvoid  volume —  l^total  volume  F solids  volume 


(54) 


Substituting  the  relationship  given  by  eq  54  into  eq  49  gives  the  following  expressions: 


8  —  1  F total  volume  ^solids  volume)/F total  volume 


1  [(^solids  volume)/!^ total  volume)] 


(55) 


With  the  volume  of  the  solids  given  by  eq  52  and  the  total  volume  of  the 
FCC  unit  cell  given  by  eq  53,  these  relationships  can  be  substituted  into  eq  55,  which 
yields  the  following: 


8  =  1  -  { [(2/3)  7i  d3]  /  (23/2  d3) } 

=  1  -  (2‘1/2  3'1 7t)  =  1  -  [tt  /  (3  21/2)] 
=  0.2595 


3.2.3  Ideal  FCC  Void  Fraction  Considerations 

The  following  observations  are  made  with  respect  to  this  modeling 

consideration: 

a.  Assuming  an  ideal  spherical  geometry  for  the  lattice  locations  of  the 
FCC  unit  cell,  the  void  fraction  is  independent  of  the  diameter  of  the 
sphere. 

b.  The  FCC  packing  is  assumed  to  be  ideal  and  represents  a  most 
efficient  packing  arrangement  of  ideal  homogeneous  spheres. 

c.  The  0.2595  void  fraction  for  this  idealized  case  compares  favorably  to 
the  0.39  void  fraction,  which  is  used  in  the  analytical  predictions  for 
12  x  30  mesh  carbon  granules. 


3.2.4  Geometric  Parameter  (a) 

The  estimation  of  a  mass  transfer  zone  (MTZ)  thickness  per  the  Gamson, 
Hougen,  and  Thodos  correlation  (ref  2)  contains  the  parameter  a.  The  Dole  and  Klotz 
article  (ref  3)  indicates  that  a  is  the  “superficial  surface  of  the  granule  (ignoring  pore 
structure)  per  unit  volume  of  bed”.  Similarly,  in  Perry’s  Chemical  Engineering  Handbook 
(ref  4;  pages  16-19),  a  similar  parameter,  ap,  is  given,  where  ap  is  the  “outer  surface  area 
of  particles  per  unit  volume  of  contacting  system  (solid  plus  fluid),  equal  to  6  (1  -  £)/dp.” 
The  analysis  that  follows  continues  with  the  ideal  packing  of  homogeneous  spheres  in  a 
close-packed  FCC  lattice  arrangement. 
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The  following  relationship  can  be  obtained  from  eq  55: 


ktotal 


volume 


=  [1  /(I  -  8)]  Vsoiids 


volume 


(56) 


The  surface  area  of  the  ideal,  uniform  spheres  within  an  FCC  unit  cell 
arrangement  is  determined  as 

(4  spheres/unit  cell)  x  (4  n  r2/sphere)  =  16  n  r2  =  16  n  (cl/2)2 


or  (57) 

A  Soiids  =  4  71  <r/2/unit  cell 

The  volume  of  the  solids  in  an  FCC  lattice  arrangement  was  previously  determined  and  is 
given  by  eq  52  as 


V* solids  =  (2/3)  7i  d3/unit  cell 

Substituting  eq  52  into  eq  56  yields  the  following  relationship: 

V total  volume  —  [1  /(I  -  8)]  { (2/3)  7i  d3/unit  cell }  (58) 

The  expressions  given  by  eqs  57  and  58  are  then  substituted  into  the 
definition  of  parameter  a  (or  ap)  as  given  in  the  Klotz  reference  and  the  Perry  and 
Chilton’s  Chemical  Engineers  ’  Handbook ,  to  yield  the  following  series  of  relationships: 

3  =  3p  =  ^fsolids/Ftotal  volume 

=  {4  7i  r/2/unit  cell }/{  [1  /(I  -  8)]  {(2/3)  71  d3/unit  cell}  }  (59) 

=  [4x3(1-  £)]/(2  d)  =  6  (1  -  £)/d 
Equation  59  is  the  targeted  relationship. 


3.3  Bulk  Density  Considerations 

There  is  a  need  to  differentiate  between  the  density  of  the  adsorbent  and 
the  bulk  density  of  the  adsorbent  packing.  The  two  are  related  through  the  void  fraction. 
The  intent  of  this  section  is  to  detail  the  specific  analytical  relationship. 

The  density  of  the  adsorbent  is  simply  the  mass  of  the  adsorbent  granules 
divided  by  the  volume  that  the  particles  would  displace  exclusive  of  the  void  volume. 
Because  the  adsorbent-loading  parameters  (which  are  described  in  the  following  sections) 
are  based  on  the  amount  of  the  adsorbent,  packing  considerations  must  be  considered. 
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Activated  carbon  adsorbents  are  comprised  of  a  heterogeneous  mixture  of 
irregularly  shaped  particles.  Each  of  these  particles  in  turn  has  an  intraparticle  network  of 
pores.  As  used  in  the  following  development,  the  adsorbent  density  is  taken  to  be  the 
mass  of  the  individual  particle  divided  by  the  apparent  volume  of  the  individual  particle. 
This  “apparent”  density,  as  used  here,  does  not  take  into  account  the  intraparticle  volume 
associated  with  the  porous  framework  that  exists  within  the  individual  coarse  particles. 
With  this  understanding,  an  apparent  adsorbent  density  can  be  expressed  through  the 
relationship 


Pads  =  mass  of  adsorbent/volume  of  adsorbent 

=  MadJVads 


(60) 


The  bulk  density  of  the  adsorbent  represents  the  volume  occupied  by  the 
adsorbent  as  well  as  the  void  volume  associated  with  the  packing  arrangement.  It  can  be 
described  through  the  relationship 


p b  =  mass  of  adsorbent/total  volume 


—  Tfads/Vtot 


(61) 


Equating  the  total  mass  of  adsorbent  in  eqs  60  and  61  yields 


(Pfe/ Pads)  =  Vads/Vto,  (62) 

From  eq  50,  the  ratio  of  the  adsorbent  volume  to  the  total  volume  goes  as 
(1  -  8).  Substituting  the  relationship  given  by  eq  50  into  eq  62  then  yields 


Pb  =  (1  -  S)  Pads 


(63) 


3.4  Equilibrium  Adsorbate  Loading  Parameters 

Adsorbate  loading  parameters  are  often  specified  on  a  per- volume  or  per- 
mass  basis.  The  two  are  related  through  bulk  property  characteristics  through  the 
following  development. 

The  adsorbate  loading  on  an  adsorbent  volume  basis  can  be  expressed  by 


cji  =  mass  of  adsorbate/volume  of  adsorbent 

=  Mi/V ads 


(64) 


The  adsorbate  loading  on  an  adsorbent  mass  basis  can  be  expressed  as 


Wi  =  mass  of  adsorbate/mass  of  adsorbent 
=  Mj/Mads 


(65) 
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Dividing  eq  64  by  eq  65  gives  the  following  relationship: 

C/i/Wi  =  MadJVads  =  Pads  (66) 

Substituting  the  relationship  given  by  eq  63  into  eq  66  yields 

qt/Wi  =  p/,/0  -  8) 

or  (67) 

Wi  =  qt  [(1  -  8)/p b] 

Differentiating  eq  67  with  respect  to  time  gives 

(dWj/dt)  =  o{qt  [(1  -  8)/p/d  }/dt  (68) 

If  the  void  fraction  and  the  bulk  density  are  taken  to  be  invariant,  then  eq  68  simplifies  to 

(dW/dt)  =  [(1  -  S)/p,]  {dq-Jdt)  (69) 

3.5  Uniformity  of  Adsorption  Wave  Contour  for  Strongly  Adsorbed  Materials 

Figure  5  depicts  the  contours  of  an  advancing  adsorption  wavefront  with 
the  progression  of  time  along  the  length  of  an  adsorbent  bed  for  a  strongly  adsorbed 
material.  The  contours  represent  snapshots  in  time  of  the  advancing  wavefront  where  the 
elapsed  time  t\  <  t2  <  h  <  f4.  The  position  Lr  represents  the  length  (or  depth)  of  the 
adsorbent  bed. 


Figure  5.  Depiction  of  advancing  adsorption  wavefront  with  time  for  a  strongly 

adsorbed  compound  (from  ref  5). 
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The  adsorption  wave  contour  depicted  by  t\  is  observed  to  differ  from 
those  depicted  by  to,  h,  and  /4.  This  is  because  the  adsorption  wave  has  not  had  sufficient 
time  and  duration  within  the  column  to  fully  develop.  The  contours  depicted  by  to,  h,  and 
f4  are  more  uniform  and  consistent  with  respect  to  one  another.  These  contours  reflect  a 
more  fully  developed  and  propagating  wavefront. 

A  hypothetical  breakthrough  concentration  is  depicted  as  the  dashed  line 
just  above  the  concentration  ratio  C/CQ  =  0.  The  thickness  of  the  MTZ  encompasses  that 
region  of  the  bed  length  where  the  concentration  ratio  C/C0  departs  from  a  value  of  1.0 
and  extends  to  the  dashed  line  of  the  breakthrough  concentration  ratio  along  a  fixed  time 
contour.  As  depicted  in  the  Figure  5  plot,  a  breakthrough  concentration  would  have  been 
reached  between  the  elapsed  time  periods  of  and  t4.  As  time  extends  beyond  t4,  the 
adsorbent  bed  will  eventually  reach  a  state  of  full  saturation  as  its  residual  adsorptive 
capacity  diminishes  with  time. 

There  are  three  regions  of  note  in  the  Figure  5  plot.  These  include  the 
region  in  which  a  state  of  saturation  has  been  achieved  on  the  adsorbent,  the  region  in 
which  active  mass  transfer  is  occurring,  and  the  region  beyond  the  MTZ  that  retains  full 
residual  adsorptive  capacity.  The  saturated  region  of  the  advancing  adsorption  wavefront 
is  depicted  by  the  relatively  flat  region  of  the  contour  extending  along  the  length  of  the 
column  where  the  concentration  ratio  C/C0  =  1.  This  region  of  the  adsorbent  has  reached 
its  maximum  adsorption  loading  capacity  (equilibrium  limit)  and  can  no  longer  effect  any 
reductions  in  the  vapor  phase  concentration.  The  second  or  active  MTZ  region  is 
characterized  by  the  downward-sloping  portions  of  the  time  contours.  The  width  of  the 
time  contours  along  the  length  axis  correlates  to  the  thickness  of  the  MTZ  that  is  bounded 
by  the  relative  concentration  ratios  established  by  C/C0  <  1  and  C/,/C(l  =  1 .  This  is  the 
region  in  which  there  is  an  active  transfer  of  adsorbate  from  the  vapor  phase  to  a 
condensed  or  adsorbed  phase.  The  region  to  the  right  of  the  MTZ  represents  the  segment 
of  the  adsorbent  bed  thickness  that  retains  residual  adsorption  capacity  provided  that  a 
concentration  breakthrough  condition  has  not  been  reached. 

In  the  following  analytical  development,  it  is  assumed  that  there  is  a 
singular  MTZ  thickness  that  fully  characterizes  the  propagation  of  the  adsorption 
wavefront  as  it  traverses  the  thickness  of  the  bed.  This  assumption  is  taken  to  apply  to 
adsorbent  beds  of  a  constant  cross-sectional  flow  area  as  well  as  for  cylindrical  systems 
with  a  positionally  dependent  cross-sectional  flow  area.  This  notion  of  a  constant  MTZ 
thickness  as  the  adsorbent  bed  is  traversed  represents  an  idealized  conceptualization.  It 
further  assumes  that  sufficient  bed  depth  is  available  in  which  a  uniform  adsorption 
wavefront  can  fully  develop. 
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4. 


DEVELOPMENT  OF  A  COMPONENT  MASS  BALANCE  FOR 
A  REACTING  SYSTEM  IN  A  RADIAL  FLOW  SYSTEM 


In  Section  2,  it  was  demonstrated  that  analytical  complexity  arises  upon 
retention  of  the  flow  characterizations  associated  with  each  directional  flow  with 
commonly  encountered  geometries.  In  Section  3,  a  foundation  was  provided  for  the 
parameters  often  encountered  in  an  adsorbent -based  filtration  system.  In  Section  4,  a 
differential  mass  balance  is  developed  for  a  reacting  system  with  radially  directed  flow  in 
a  cylindrical  coordinate  system.  This  development  builds  upon  the  approaches 
established  earlier  and  incorporates  the  parameters  associated  with  an  adsorbent  fill. 
Figure  6  depicts  the  system  that  was  considered  in  this  development. 


Co 


A  r 


n  r  r  +  Ar  ri 


Figure  6.  Radial  flow  adsorbent  system. 


Referring  to  Figure  6,  the  bed  depth  of  the  radial  flow  filter  is  (j2  -  ri),  and 
a  challenge  concentration  of  species  i  (C0)  is  introduced  into  the  filter  at  radial  position  r\. 
An  effluent  concentration  of  species  i  (Ce)  exits  the  filter  at  radial  position  r2.  It  is 
assumed  that  species  i  is  consumed  within  the  filter  at  a  molar  rate  depicted  by  where 
Ri  is  defined  on  a  per- volume  basis. 

The  height  of  the  cylinder  taken  along  the  axial  coordinate  is  h.  The  cross- 
sectional  area  associated  with  the  principle  flux  direction  is  therefore  A  =  2n  r  h. 


The  volume  of  the  differential  element  depicted  in  Figure  6  is  given  as 


AVcyl  =  (AS)  (Az)  (Ar) 

=  (r  A9)  (Az)  (Ar) 


(70) 


In  accordance  with  the  arguments  presented  in  Section  2.3,  it  is  assumed 
that  the  flows  associated  with  the  angular  (9)  and  axial  (z)  directions  are  negligible 
relative  to  the  radially  directed  flow,  which  is  taken  to  be  the  principle  flow  direction. 
Given  this  premise,  the  volume  of  the  differential  element  depicted  in  Figure  6  and 
expressed  by  eq  70  can  be  expressed  by 
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AVcyi =  (271  r  h)  A r 


(71) 


or 


AVcyi  =  A  A r 


4.1  Rate  of  Mass  Flow  of  Species  i  Entering  the  Differential  Volume 

Element  at  Radial  Position  r 

An  analytical  relationship  for  the  rate  of  mass  flow  of  species  i  entering 
the  differential  volume  element  will  be  developed  in  terms  of  a  mass  flux.  The  analytical 
expression  is 


md in  =  (MW);  INi(r)  A(r)]adsorbent  (72) 

The  mass  flux  given  by  eq  72  represents  a  total  or  effective  flux  that  exists 
within  the  adsorbent  packing.  The  area  that  appears  in  eq  72  is  more  appropriately 
associated  with  the  area  through  which  species  i  is  transferred  within  the  vapor  phase 
through  the  packing  in  the  radial  direction.  As  such,  this  area  represents  a  restricted  area 
due  to  the  presence  of  the  solid  adsorbent  packing.  The  appropriate  area  expression  in 
eq  72  should  be  the  restricted  area.  With  this  notion,  eq  72  is  more  appropriately 
expressed  as 


(ITldin  =  (MW);  (V;  A)adsorbent  =  (MW);  N,{r)  8  A(r) 


(73) 


=  (MW);  8  Ni(r)  A(r) 


4.2  Rate  of  Mass  Flow  of  Species  i  Exiting  the  Differential  Volume  Element 

at  Radial  Position  r  +  A r 

By  analogy,  with  the  development  of  eq  73,  the  rate  of  mass  flow  of 
species  i  which  exits  the  differential  volume  element  at  radial  position  (r  +  A r)  can  be 
expressed  as 


(m,)out  =  (MW);  \Nj(r  +  A r)  A(r  +  Ar)]  adsorbent  (74) 

Adjusting  the  area  to  account  for  the  packing  of  adsorbent,  eq  74  takes  on  the  following 
forms: 


(llld out  =  (MW);  [N,{r  +  Ar)  8  A(r  +  Ar)] 
=  (MW);  8  [V;(r  +  Ar)  A(r  +  A/-)] 


(75) 
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4.3 


Mass  of  Species  i  within  the  Differential  Volume  Element 


The  mass  of  species  i  within  the  differential  volume  element  is  comprised 
of  the  mass  in  the  adsorbed  state  and  the  vapor  phase.  An  analytical  expression  that 
characterizes  the  mass  of  species  i  within  the  differential  volume  element  is  developed 
through  the  following  relationships: 

Mf  —  (M i) adsorbed  phase  (Af,')v apor  phase 

=  Vtot  (^ads/l^tot)  [(4/;)vap/ EvoidJ  Cvoid 

(76) 

=  [A(r)  Ar]  (1  -  8)  q,  +  (MW),  C,  (8  Vtot) 

=  [. A(r )  Ar]  (1  -  8)  qt  +  (MW),  C,  8  [A(r)  Ar] 

The  molar  concentration  of  species  i  in  the  vapor  phase  is  based  on  the 
void  or  interstitial  volume  within  the  adsorbent  packing. 


4.4  Rate  of  Change  of  Mass  of  Species  i  within  the  Differential  Volume 

Element 

An  expression  relating  the  rate  of  change  in  mass  of  species  i  within  the 
differential  volume  element  can  be  obtained  by  differentiating  eq  76  with  respect  to  time. 
This  differentiation  is  accomplished  as  follows: 


(dMi/dt)  =  d[A  Ar  (1  -  8)  q<]ldt  +  d[(MW),  C,  8  A  Ar]/8t 

=  (1  -  8)  A r  [6(A  qi)/dt]  +  (MW),  8  Ar  [d(A  Q)/dt] 


(77) 


The  relationship  given  by  eq  77  was  obtained  by  assuming  that  the  void 
fraction  is  constant  and  the  derivatives  of  the  differentials  are  negligibly  small  relative  to 
the  first-order  derivatives.  Expanding  the  derivatives  in  eq  77  yields  the  following 
equivalent  form  of  the  expression: 


( dMi/dt )  =  (1  -  8)  Ar  [A  ( dqjdt )  +  q,  (dAldt)\ 

+  (MW),  8  Ar  [A  (6C,/ot)  +  C,  (QA/dt)] 


(78) 


4.5  Rate  of  Consumption  of  Species  i  Due  to  a  Chemical  Reaction  on  the 

Adsorbent  Surface 

A  generic  reaction  rate  term  is  used  to  describe  the  rate  of  disappearance 
of  species  i  on  the  adsorbent  surface.  This  reaction  rate  term  is  based  on  a  per- volume 
basis  and  is  a  function  of  several  variables.  A  suggested  form  of  the  rate  of  reaction 
relationship  is  given  by  eq  79  as 
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(m,)rxn  =  (MW),  R!  (A  Ar) 


(79) 


4.6  Analytical  Mass  Balance  Development  for  Species  i 

The  component  mass  balance  expression  can  be  expressed  as 


(dM-Jct)  =  -  Wi) out  -  (m,)rxn  (8°) 

Equation  80  is  a  relationship  relating  the  accumulation  of  mass  of  species  i 
within  the  differential  volume  element  to  the  mass  input,  the  mass  output,  and  a  mass 
consumption  term.  The  dynamic  loading  and  desorption  of  species  i  to  and  from  the 
surface  of  the  adsorbent,  respectively,  are  addressed  through  the  surface-loading  term 
appearing  in  eq  78. 

Substitution  of  eqs  73,  75,  78,  and  79  into  eq  80  yields  the  following 

relationship: 

(1-8)  A r  [A  ( dqjdt )  +  c/,  (i dA/dt )]  +  (MW);  8  A r  [A  (dCJdt)  +  C,  (dA/dt)] 

(81) 

=  (MW);  8  Nj(r)  A(r)  -  (MW);  8  [A,(r  +  Ar)  A(r  +  Ar)]  -  (MW),  R,  (A  Ar) 


Division  of  eq  81  by  [(MW);  8  Ar)]  yields 

[1/(MW),]  [(1  -  S)/8]  [A  (oqjdl)  +  qt  (dA/dt)]  +  [A  (dQ/dt)  +  C,  (dA/dt)] 

=  ~  {  { [ Nj( r  +  Ar)  A(r  +  Ar)]  -  N) r)  A(r)  }/A r  }  -  (R,/E)  A(r) 
In  the  limit  as  Ar  collapses  to  zero,  eq  82  takes  on  the  form 
[1/(MW ),]  [(1  -  8)/S]  [A  (dqi/dt)  +  qt  (dA/dt)]  +  [A  (dQ/dt)  +  C,  (dA/dt)] 


(82) 


(83) 


=  -[o(N,  A)/ dr ]  ~  (Ri/ 8)  A(r) 

Noting  the  following  relationships  for  a  cylinder  with  a  constant  height  h , 

A  =  2n  h  r 

(dA/dt)  =  271  h  ( dr/dt ) 

[d(Nj  A)/ dr]  =  271  h  [A;  +  r  (dN)dr)] 

Equation  83  can  therefore  be  expressed  in  the  following  equivalent  form: 


and 
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[1/(MW),J  [(1  -  8)/8]  [r  (dqjdt)  +  q,  (Gr/ot)\  +  [r  (, dCJdt )  +  C,  (dr/dt)\ 


(84) 


=  -Ni  -  r  ( dNi/dr )  -  r  (/?,/£) 
or 

[l/(MW)z]  [(1  -  8)/8]  [d(r  q)/dt]  +  [d(r  C,)/dt\  +  [d(r  N,)/dr\  +  r  (RJE)  =  0  (85) 


4.7  Comparison  of  Equation  85  to  Reference  Relationship 

It  is  instructive  at  this  point  to  compare  eq  85  to  a  reference  relationship 
that  is  provided  in  a  modern  text  on  the  analytical  aspects  associated  with  adsorption 
processes.  Such  a  reference  expression  is  provided  by  eq  7.5  in  the  text  Principles  of 
Adsorption  and  Adsorption  Processes  by  D.  M.  Ruthven  (ref  6).  This  equation  served  as 
a  comparative  reference  for  many  of  the  analytical  developments  contained  in  this  review 
effort.  The  Ruthven  relationship  (ref  6)  is  consistent  with  the  following  relationship: 

-D L  ( &CJdz 2)  +  [d(v  Ci)/8z\  +  ( dQ/dt )  +  [1/(MW),J  [(1  -  s)/s]  (dqjdt)  =  0  (86) 

The  Ruthven  expression  is  based  on  plug  flow  in  a  singular  direction, 
which  is  characterized  by  a  uniform  dispersion  coefficient  D/.  The  relationship  is  further 
based  on  a  nonreacting  system  having  a  constant  cross-sectional  flow  area. 

For  a  constant  cross-sectional  area,  the  r  terms  within  eq  85  can  be 
considered  as  constants  and  can  therefore  be  removed  from  the  derivative  terms.  The  r 
variables  in  eq  85  would  essentially  divide  out,  leaving  the  following  relationship: 

[1/(MW),J  [(1  -  eye]  (dqjdt)  +  ( dc/dt )  +  (dNi/dr)  +  (RJ 8)  =  0  (87) 

In  a  nonreacting  system,  Rj  =  0,  and  eq  87  becomes 

[1/(MW),J  [(1  -  S)/S]  (dqjdt)  +  ( dCJdt )  +  ( GNJGr )  =  0  (88) 

The  molar  flux  gradient  (dNJdr)  is  replaced  with  (dNJdz)  because  an  axial 
system  with  a  constant  cross-sectional  area  is  being  considered.  Equation  88  therefore 
becomes 


[1/(MW)/]  [(1  -  8)/8]  (dqjdt)  +  (dCJdt)  +  (dNJdz)  =  0  (89) 

The  relationship  represented  by  eq  89  corresponds  to  an  axial  flow  equivalent  to  eq  85  for 
a  nonreactive  system  having  a  constant  cross-sectional  area. 

The  total  flux  with  respect  to  the  axial  direction  is  based  on  eq  10.  For  the 
particular  system  under  consideration,  the  following  relationship  applies: 
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Niz  =  vz  Ci  ~  B(-  (dC-Jdz) 


(90) 


Assuming  that  a  homogeneous  coefficient  (Bz,)  can  be  used,  and  removing 
the  directional  subscript  from  the  molar  flux,  eq  90  can  be  re-expressed  as 

Ni  =  vzCi-BL(dCi/dz)  (91) 

Bifferentiating  eq  91  with  respect  to  z  and  assuming  that  the  dispersion 
coefficient  is  constant-valued  yields 

(dNi/dz)  =  [d(vz  Ci)/dz ]  -  Bl  (d2Ci/dz2)  (92) 

Substituting  the  expression  given  by  eq  92  into  eq  89  yields 

-Bl  (d2Cj/dz2)  +  [d(v-  Q)/dz \  +  ( dCJdt )  +  [1/(MW),J  [(1  -  8)/8]  (dqjdt)  =  0  (93) 

Equation  93  is  the  Ruthven  relationship  as  given  by  eq  86. 


4.8  Molar  Flux  Term  of  Equation  85 

This  term  can  be  expanded  through  the  generic  flux  relationship  given  by 
eq  10.  From  eq  10,  the  following  expression  applies  for  the  radially  directed  molar  flux: 

Ni  =  vr  Ci  -  B L  ( dCi/dr )  (94) 

Using  the  molar  flux  relationship  given  by  eq  94,  the  following  partial 
derivative  is  obtained: 


[G(r  Nj)!or\  -  vr  Ci  +  r  vr  ( dCJdr )  +  r  Ci  ( dvjdr ) 

-  r  Bl  ( d2Ci/dr2 )  ~  Bl  ( dCi/dr ) 


(95) 


The  relationship  given  by  eq  95  can  be  substituted  into  eq  85,  which 
would  yield  a  rather  complex  expression. 


4.9  Basis  of  the  Adsorbent  Loading 

It  is  often  useful  to  express  the  adsorbent  loading  term  on  a  mass  basis 
such  as  Wj.  This  is  done  because  gravimetric  methods  are  often  more  precise  than 
volumetric-based  determinations.  The  relationship  between  the  loading  terms  qt  and  Wi 
was  developed  in  eq  67.  Using  eq  67  as  a  basis,  the  following  relationship  can  be  written: 

r  q,  =  lPb  /  (1  ^S)J  (rWi)  (96) 
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Assuming  that  the  bulk  density  and  the  void  fraction  are  constant,  eq  96  is 
differentiated  with  respect  to  time  to 

[d{r  qi)/dt ]  =  [p*  /  ( 1  -  £)]  [d{r  W)/dt]  (97) 

Substituting  the  relationship  given  by  eq  97  into  eq  85  yields 

[1/(MW),J  (p/,/8)  [c(r  Wi)/dt\  +  [5(r  C,)/8t ]  +  [d(r  N)/Gr\  +  r  (7 ?,-/£)  =  0 


4.10  Notes  on  idA/dt),  idr/dt ),  and  Positional  References 

The  rate  of  change  of  mass  within  the  differential  volume  element  can  be 
considered  from  two  reference  perspectives.  In  one  perspective,  the  volume  element  is 
spatially  displaced  from  position  r  to  (r  +  Ar).  The  analytical  development  with  this 
frame  of  reference  gives  rise  to  the  idA/dt )  terms  of  eq  78.  The  rate  of  change  of  mass 
within  the  volume  element  with  respect  to  time  and  displacement  would  be  given  by 
eq  78. 


With  the  second  reference  frame,  one  is  interested  in  the  rate  of  change  of 
mass  relative  to  a  fixed  spatial  reference.  This  situation  would  arise  if  one  were 
monitoring  the  effluent  properties  relative  to  a  fixed  reference  location.  With  respect  to  a 
fixed  spatial  reference,  the  area  would  be  constant-valued  and  the  derivative  idA/dt )  =  0. 
With  [dA/dt )  =  0,  it  follows  that  idr/dt)  =  0.  With  idr/dt)  =  0,  eq  84  can  be  expressed  as 

[1/(MW)(]  [(1  -  £)/£]  [r  idqjdt )]  +  [r  (dQ/dt)]  +  [d(r  N)/dr ]  +  r  (R,/£)  =  0  (98) 


or 


[1/(MW)/]  (p/,/8)  [r  idWJdt)\  +  [r  (dC,/dt)\  +  [d[r  N)/dr\  +  r  (/?,/£)  =  0  (99) 

Information  relating  to  partial  and  total  derivatives  as  well  as  reference 
frames  is  expanded  upon  in  Appendix  C. 


4. 1 1  Closing  Commentary 

The  utility  of  the  various  forms  of  eqs  85  and  99  for  practical  engineering 
design  purposes  is  formidable.  The  causes  of  this  analytical  complexity  stem  from  the 
need  to  characterize  the  molar  flux  with  position  and  loading  parameters  with  respect  to 
time,  the  need  to  characterize  the  reaction  rate  in  terms  of  its  independent  variables,  and 
the  relative  motion  of  the  volume  element  with  respect  to  a  fixed  spatial  reference  (if 
applicable).  The  dynamics  associated  with  the  changing  cross-sectional  area  can  be 
simplified  somewhat  by  consideration  of  the  reference  frame;  however,  the  resultant 
expressions  still  retain  a  degree  of  analytical  complexity. 
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It  must  be  recalled  that  the  analytical  development  contained  in  this 
section  addressed  the  radial  components  of  the  flux  and  assumed  that  the  molar  fluxes 
and  gradients  associated  with  the  nonprinciple  flow  directions  were  negligible  relative  to 
their  counterparts  in  the  radial  direction. 

For  very  precise  efforts,  computational  algorithms  can  be  developed  to 
arrive  at  numerical  solutions  to  equations  of  the  form  of  eqs  85  and  99.  For  practical 
engineering  purposes,  much  simpler  analytical  approaches  are  often  used;  however,  these 
simplified  approaches  are  necessarily  limited  in  their  ability  to  truly  model  the 
performance  characteristics  of  actual  filters.  It  must  also  be  remembered  that  the 
energetic  aspects  associated  with  physical  adsorption  and  chemisorption  were  not 
considered  in  this  review  effort. 


5 .  EFFORT  BY  BOHART  AND  ADAMS 

An  early  effort  to  characterize  the  performance  of  carbon  against  an 
airstream  contaminated  by  chlorine  was  undertaken  by  Bohart  and  Adams  (ref  7).  This 
early  work  was  extensive,  as  Bohart  and  Adams  conducted  numerous  experimental  trials 
and  developed  analytical  expressions  in  an  effort  to  characterize  the  breakthrough 
characteristics  of  chlorine  through  carbon  beds.  Their  analytical  expressions  provided  a 
means  to  predict  the  concentration  of  chlorine  with  both  position  and  time.  With  a  room 
temperature  vapor  pressure  considerably  higher  than  10  torr,  chlorine  is  removed  through 
a  chemisorption  process  within  the  adsorbent  system. 

Bohart  and  Adams  (ref  7)  assumed  a  simplistic  rate  expression  in  which 
the  rate  of  loss  of  residual  capacity  of  the  carbon  for  the  chlorine  challenge  vapor  ( da/dt ) 
was  proportional  to  the  product  of  a  first-order  rate  constant  (k),  the  localized  residual 
capacity  of  the  carbon  for  the  chlorine  (a),  and  the  localized  vapor  phase  concentration  of 
the  chlorine  with  axial  position  (C,).  Their  fundamental  rate  expression  took  on  the  form 

(da/dt)  =  —k  a  Ci  (100) 

As  the  ability  of  the  carbon  to  remove  chlorine  diminished  with  time,  the 
residual  loading  capacity  of  the  adsorbent  would  decrease  in  accordance  with  eq  100.  At 
a  condition  of  full  consumption,  or  saturation  of  the  adsorbent,  the  residual  loading 
capacity  would  have  a  value  of  zero. 

To  describe  the  change  in  chlorine  concentration  within  the  vapor  phase  as 
a  given  fluid  element  traversed  the  bed  in  an  axial  direction,  Bohart  and  Adams  proposed 
the  following  rate  relationship: 


(dCi/dz)  =  -(k/v)aCi  (101) 

where  v  represents  the  rate  of  supply  of  chlorine  that  saturates  the  carbon  bed.  As  used  in 
this  definition,  v  can  be  likened  to  the  propagation  of  the  adsorption  wave  as  it  traverses 
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the  adsorbent  column.  It  is  not  directly  related  to  the  superficial  velocity  of  the  airflow 
that  enters  the  column. 

Some  distinction  on  the  use  of  adsorbate  loading  parameters  should  be 
emphasized  here.  The  loading  parameters  <y,  and  Wu  as  used  in  previous  sections,  related 
to  the  loading  of  the  adsorbent  with  time.  Their  values  would  generally  increase  with 
time.  Bohart  and  Adams  based  their  analytical  development  on  the  notion  of  a  residual 
loading  capacity  that  would  tend  to  decrease  with  time.  As  such,  the  parameter  a  is  to  be 
differentiated  from  the  parameters  and  Wu 

The  relationship  given  by  eq  101  relates  the  change  in  concentration  of  the 
chlorine  concentration  in  the  vapor  phase  to  an  axial  position  (constant  cross-sectional 
area  is  assumed).  Inspection  of  eq  101  indicates  a  similarity  to  the  steady-state  design 
equation  used  in  simplified  expressions  for  an  ideal  plug  flow  reactor  (PFR)  system  with 
first-order  reaction  kinetics.  With  ideal  PFR  systems,  the  flow  effects  associated  with 
dispersion  and  deviations  from  ideal  flow  behavior  (i.e.,  axial  mixing  of  adjacent  volume 
elements)  are  not  considered.  Further  assumptions  of  an  ideal  PFR  system  are  that  each 
volume  element  takes  the  same  amount  of  time  to  traverse  the  reactor  from  the  inlet  to  the 
outlet  and  that  there  is  uniformity  throughout  the  element  in  the  radial  direction.  With  the 
ideal  steady-state  design  relationship  for  such  systems,  composition  changes  in  the  vapor 
phase  concentration  of  i  are  spatially  dependent  and  time  invariant.  A  representative 
steady-state  PFR  design  relationship  can  be  expressed  by 

(dC/dz)  =  ~(A/Q)  Ri  (102) 

Assuming  that  the  reaction  rate  term  is  first-order  in  concentration,  then  /?,  =  k  C„  and 
eq  102  takes  on  the  form 


(dCi/dz)  =  ~(A/Q)  k  Ct  (103) 

The  generic  similarities  between  the  ideal  PFR  design  expression  and  eq  101  are  noted. 

Through  a  series  of  variable  transformations,  magnitude  estimations,  and 
rationalizing  arguments,  Bohart  and  Adams  give  the  following  relationships  to 
characterize  the  transmission  of  chlorine  through  carbon  beds  of  constant  cross-sectional 


area: 

(C,IC0)=  1 0'  /(  1 0Z  -  1  +  10'") 

(104) 

and 

(< a/a0 )  =  10z"/(10z"  -  1  +  10'") 

(105) 

where 

z"  =  0.4343  k  a0  (z/v) 

(106) 

and 

t"  =  0.4343  k  C0  t 

(107) 
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Equations  104  and  105  provide  a  means  to  calculate  the  concentration  of 
chlorine  vapor  in  the  airstream  as  a  function  of  both  axial  position  and  time,  provided  that 
the  rate  parameters  can  be  determined.  The  rate  constant  would  be  determined  through 
experimentation. 


5.1  Analytical  Development  of  the  Bohart  and  Adams  Breakthrough 

Relationships 

The  following  development  details  the  analytical  methodology  taken  by 
Bohart  and  Adams  (ref  7)  in  arriving  at  the  relationships  given  by  eqs  104  and  105.  The 
development  is  based  on  content  contained  within  the  appendix  of  the  cited  reference. 

The  following  variable  transformations  are  defined: 


a'  =  a/aQ 

(108) 

Ci'  =  Ci/C0 

(109) 

'  =  kaQ  (z/v) 

(HO) 

II 

(HI) 

With  these  variable  transformations,  the  following  relationships  are 
developed  for  {dCj'/dz1)  and  (da'/dt'): 

0 dCi'/dz ')  =  { dCi'/dCi )  {dQ/dz)  (dz/dz')  =  (1  /C0)  { dQ/dz )  [v/(k  a0)\ 

=  (1  /C0)  [(-k  a  Q)/v]  [v/(k  a0)\  =  ~{a/a0)  { Q/Q )  (112) 

=  -a'c; 

{da'/dt')  =  {da' /da)  {da/dt)  {dt/dt’)  =  (1  /a0)  {8a/ dt)  [1  /{k  C„) \ 

=  { Had)  {~k  a  Cd  [  \/(k  C0)]  =  ~{ata0)  {CJC0)  (113) 

=  -a'Q ' 


The  following  expressions  result  directly  from  eqs  112  and  113: 

[{dCi'/Qydz']  =  [d{ In  Ci')/dz']  =  -a'  (114) 

and 

\{da7a')/df]  =  [5(ln  a’)/ of]  =  -Cf  (115) 
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The  further  development  of  eqs  114  and  1 15  is  based  on  their  integrated 
forms.  Bohart  and  Adams  rationalize  the  integration  of  these  relationships  as  follows:  At 
time  t'  =  0  (corresponding  to  l  =  0  for  nonzero  values  of  k  and  C0),  the  residual  loading 
capacity  of  the  carbon  for  chlorine  is  taken  to  be  uniform  and  constant  valued  at  a0 
throughout  the  entire  length  of  the  carbon  bed.  This  uniformity  of  the  residual  loading 
capacity  would  apply  for  an  adsorbent  system  that  does  not  have  any  prior  contaminant 
loading  history  or  operational  service.  In  other  words,  the  adsorbent  is  fresh  and  assumed 
to  have  uniform  loading  capabilities  throughout  the  column.  The  adsorbent  possesses 
maximum  loading  capability  in  this  condition.  Therefore,  a'  =  1  over  the  entire  length  of 
the  adsorbent  column.  The  transformation  length,  z',  as  defined  by  the  relationship  given 
by  eq  1 10,  is  observed  to  be  a  function  of  axial  position  only  (the  other  parameter  values 
being  multiplicative  constants).  Because  a' is  assumed  to  be  constant  valued,  eq  1 14  can 
be  integrated  over  the  axial  length  of  the  carbon  bed  in  accordance  with  the  following 
development: 


Ci@  z'  z' 

fd(lnC’)  =  -fa'dz'  (116) 

Cj'@z'= 0  z'=  0 

The  integral  on  the  left  side  of  the  equality  in  eq  1 16  is  integrated  as 

In  (C’@  z')  :  An  (C/@  z'  =  0)  =  In  {CJC0)  I  -  In  (C,IC0)  I 

■z'  V  =  0 

At  position  z  =  0  (corresponding  to  z'  =  0),  the  concentration  of  chlorine 
(or  component  i)  in  the  vapor  phase  is  CQ.  Therefore, 

In  (Ci/C0)  I  =  In  1  =  0 

h'  =  0 

The  integral  on  the  right  side  of  the  equality  of  eq  1 16  is  simply  (-a1  z')  for 
a  constant  value  of  a'.  At  t  =  0  (corresponding  to  t'  =  0  for  nonzero-valued  parameters  k 
and  Co),  then  a'  =  1.  With  a'  =  1,  the  integrated  form  of  eq  1 16  is  given  as 

Ci  =  e~z'  (117) 

At  position  z'  =  0  (corresponding  to  z  =  0  for  nonzero  values  of  k,  aQ,  and 
v),  the  inlet  concentration  of  the  chlorine  is  uniform  and  constant  valued  at  C0  with  time. 
Therefore,  C/  =  1  over  any  time  duration.  The  transformation  time  t\  as  defined  by  the 
relationship  given  by  eq  1 1 1,  is  observed  to  be  a  function  of  time  only  (the  other 
parameter  values  being  multiplicative  constants).  Since  C/  is  a  constant,  eq  1 15  can  be 
integrated  over  time  in  accord  with  the  following  development; 

a'@f  t' 

fdilna')  =  -fCi'dt'  (118) 

a'@t'  =  0  t'  =  0 
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The  integral  on  the  left  side  of  the  equality  in  eq  1 18  is  integrated  as 


In  (a1  @  t’)  -  In  (a'  @  t'  =  0)  =  In  (a/a0)  -  In  ( a/aQ ) 

V  Y  =  0 

At  time  t  =  0  (corresponding  to  f  =  0),  the  residual  capacity  of  the 
adsorbent  for  chlorine  is  aQ .  Therefore 

In  ( a/aD )  I  =  In  ( aja0 )  =  In  (1)  =  0 

V  =  0 

The  integral  on  the  right  side  of  eq  1 18  is  simply  (C/  t')  for  a  constant 
value  of  C{.  At  z  =  0  (corresponding  to  z!  =  0  for  nonzero- valued  parameters  k,  a,,,  and  v), 
then  C{  =  1.  With  C  '  =  1,  the  integrated  form  of  eq  1 18  is  given  as 

a'  =  e~f  (119) 


5.2  Analytical  Manipulations  with  Bohart  and  Adams  Relationships 

Through  the  imposition  of  boundary  conditions,  independent  relationships 
for  the  concentration  of  chlorine  as  functions  of  axial  position  and  time  were  developed 
in  Section  5.1.  These  are  given  by  eqs  117  and  119.  The  aim  of  this  section  is  to  further 
detail  the  approach  pursued  by  Bohart  and  Adams  (ref  7)  in  the  development  of  an 
analytical  expression  that  couples  time  and  positional  dependency  to  a  common 
expression. 


Provided  that  C{  does  not  equal  0,  the  following  expression  results 
directly  from  eq  112: 


{  [dC’KC’fydz’ }  =  -ia'/CO  (120) 

Substitution  of  the  relationships  developed  for  C/  (eq  1 17)  and  a'  (eq  1 19) 
into  eq  120  gives  the  expression 

{ \dC’l(C  ’i]ldz! }  =  -  (l/er'-z')  (121) 

Separation  of  the  variables  in  eq  121  sets  up  the  following  integral 

expression: 

(C/)2  z'i 

f  (Ci)~2  dCi  =  /-(l /e'-z)dz'  (122) 

(C/)i  z'  i 
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The  1  and  2  subscript  designations  appearing  in  the  integration  limits  of  eq  122  represent 
the  dimensionless  concentration  and  axial  position  at  specific  locations.  These 
designations  are  defined  as  follows: 

Subscript  1  =  parameter  value  at  entrance  location  or  z!  =  0. 

Subscript  2  =  parameter  value  at  axial  position  z'  within  the  carbon  bed. 

The  integral  on  the  left  side  of  the  equality  given  by  eq  122  is  evaluated  as 

(CO  2  (C’h 

/ (C/)'2  dC{  =  -  (i ICO  =  [i /(COi]  -  [i/(c/)2] 

(CO  i  (CO  i 

The  integral  on  the  right  side  of  the  equality  given  by  eq  122  is  evaluated 
as 

z'l  z!l  z'  2 

/ -  (l/ef'“z)  dz'  =  ~  f  e_(('_z)  dz!  =  ~  f  e(z"°  dz' 

rj'  rjt 

Z  l  Z  l  Z  l 

-t  -f 

Z  2  Z  1 

=  -  f  ez  e~‘  dz'  =  e_/  /  ez  dz' 

Z'l  z’ 2 

=  e~f  (e^1  -  e;’2)  =  efe'1_0  -  e(z'2'° 

Therefore,  eq  122  can  be  expressed  by  the  integrated  form 

[1  l(C0i\  ~  [1  KC0i\  =  e(z  l_f)  -  e(z'2_/,)  (123) 

The  parameter  values,  which  are  positionally  dependent,  are  determined  as 

follows: 

(COi  =  (CJCoh  =  ( Co/Co )  =  1  at  z'i  =  0 

and 

(C’)2  =  (CJC0)  at  z’i  =  z! 

Substitution  of  these  parameter  values  into  eq  123  gives 

1  -  [1/(C/)J  =  e(0_f,)  -  e(z'“°  =  e~f  -  hz'~n 

or 

[l/(C/)]  =  l-e“f'  +  e(zW,)  (124) 

Equation  124  can  be  solved  explicitly  for  C/,  which  yields 
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Ci'=  ef  /  (ef  -  1  +  ez  ) 


(125) 


Recognizing  that  ee  =  lO(0/ln  10)  on  a  term-by-term  basis,  eq  125  can  be 
expressed  in  terms  of  common  logarithms  as 

d  =  10('71n  10)  /  [10('71n  10)  -  1  +  10(z71n  1())]  (126) 

Letting  t"  =  fl In  10  =  0.4343  f  and  z"  =  z'/ln  10  =  0.4343  z',  eq  126  can  be 

rewritten  as 

c;  =  (Q/C0)  =  10r  /  (10f"  -  1  +  10z")  (127) 

Therefore,  eq  127  is  the  same  as  eq  104. 

The  development  of  an  analytical  relationship  to  relate  the  residual 
adsorptive  capacity  of  the  carbon  (a)  as  a  function  of  time  and  position  is  similarly 
generated.  For  completeness,  the  development  is  as  follows. 

Provided  that  the  carbon  retains  residual  adsorptive  capacity  (i.e.,  a  ±  0), 
the  following  relationship  results  directly  from  eq  113: 

(da'/df)  =  -a'  C{ 


or 


(1  la'2)  (da'/df)  =  -Q’/a’ 


(128) 


From  the  relationships  developed  for  C{  (eq  1 17)  and  a'  (eq  1 19),  eq  128  can  be 
expressed  as 

(l/a'2)  (da'/df)  = -e(t'~z'}  (129) 

Separation  of  variables,  and  establishing  integration  limits,  eq  129  takes 
on  the  integral  form 


(a'h  fi 

f(a')~2  da'  =  f-  e(t'~z,)  df  (130) 

(a')  i  fi 

The  1  and  2  subscript  designations  in  the  integration  limits  of  eq  130  represent  the 
dimensionless  residual  chlorine  capacity  of  the  carbon  at  specific  times.  These 
designations  are  defined  as  follows: 

Subscript  1  =  parameter  value  at  time  f  =  0. 

Subscript  2  =  parameter  value  at  time  f  within  the  carbon  bed. 

The  integral  on  the  left  side  of  the  equality  of  eq  130  is  evaluated  as 
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/ 

(fl')l 


(a’h 


(a')  2  da’  =  -(1/a') 


(«')2 


(a')i 


[1/(«,)1]  -  [l/(a')2] 


The  integral  on  the  right  side  of  the  equality  of  eq  130  is  evaluated  as 
t'i  t'2 

f-e(t'~z'}  8t'  =  ~e~z'  / /  dt'  =  -oz'  (er'2  -  erl) 


t\ 


t'  i 


gCrt-zO  _  e(f'2"z') 


Therefore,  eq  130  can  be  expressed  by  the  integrated  form 


[l/(a')1]-[l/(a')2]  =  el 


1  _  Jt' l-z')  Jt'2-z') 


-  e 


(131) 


The  following  parameter  values,  which  are  assumed  to  be  spatially 
independent,  are  calculated  based  on  the  specific  times  established  previously: 


and 


(a')i  =  1  at  t'\  =  0 
(a')2  =  a'  at  t'2  =  t' 

Substitution  of  these  parameter  values  into  eq  131  gives  the  relationship 

1  -  (1/a')  =  e(0_z  )  -  e(t'~z) 

=  Q~Z'  -  e(t'~  Z"> 


or 


(1/a')  =  e(t'~z  )  +  1  -  q~z' 

Equation  132  can  be  solved  explicitly  for  a'  to  yield 


a'  =  e"  /  (ef  -  1  +  e; ) 


(132) 


(133) 


With  e9  =  10(9/ln  10)  on  a  term-by-term  basis,  eq  133  can  be  expressed  in 
terms  of  common  logarithms  as 

a’  =  10fz7ln  10)  /  [10a71n  10)  -  1  +  10(z71n  10)]  (134) 

Letting  t"  =  t'l In  10  =  0.4343  t',  and  z"  =  z'/ln  10  =  0.4343  z\  eq  134  can  be 


rewritten  as 


a'  =  (a/a0)  =  10z"  /  (K)'"  -  1  +  10z") 
Equation  135  is  the  same  as  eq  105. 


(135) 
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5.3 


Simplifications  to  the  Bohart  and  Adams  Breakthrough  Relationships 


The  breakthrough  relationships  developed  by  Bohart  and  Adams,  given  by 
eqs  104  and  105,  can  be  approximated  by  the  following  relationships  for  situations  in 
which  z"  and  t"  are  appreciably  large: 


and 


c;  =  acjCo)  =  i  /  [io(z"_r)  + 1] 


(136) 


a'  =  (a/a0)  =  1  —  C{ 


(137) 


The  mathematical  rationale  behind  these  simplifications  is  that  10z  and 
10f  would  have  values  appreciably  larger  than  1.  As  such,  the  1  appearing  in  the 
denominators  of  eqs  104  and  105  could  be  neglected. 


Very  small  values  of  both  z"  and  t"  would  infer  a  location  within  the 
carbon  adsorbent  that  is  very  close  to  the  inlet  and  a  very  brief  time  (based  on  the 
definitions  given  for  these  parameters).  These  conditions  would  generally  not  allow  for 
the  formation  of  a  consistent  adsorption  wavefront. 


5.4  Closing  Commentary  on  the  Bohart  and  Adams  Analytical  Development 


The  Bohart  and  Adams  analytical  development  is  based  on  a  first-order 
adsorption  rate  expression  in  which  ideal  plug  flow  is  assumed.  With  ideal  plug  flow, 
discrete  fluid  elements  are  assumed  to  traverse  the  reactor  length  from  inlet  to  outlet,  and 
there  is  no  interaction  among  the  adjacent  fluid  elements.  The  effects  of  axial  dispersion 
are  neglected  in  this  idealization.  At  steady  state  (steady  in  that  the  vapor  phase 
concentration  of  i  remains  constant  relative  to  a  fixed  spatial  reference),  the  concentration 
of  species  i  is  positionally  dependent  and  time  invariant.  This  aspect  is  analytically 
reflected  in  eq  101.  Bohart  and  Adams  attempt  to  couple  this  ideal  plug  flow  behavior 
with  the  idea  that  the  residual  capacity  of  the  adsorbent  decreases  with  time,  which  is 
reflected  by  the  relationship  given  by  eq  100.  A  major  assumption  with  their  development 
is  that  the  residual  loading  capacity  of  the  carbon  for  the  chlorine  is  perfectly  uniform 
throughout  the  adsorbent  bed.  The  authors  identify  this  aspect  as  a  limitation  with  their 
modeling  expressions. 


As  indicated,  the  fundamental  Bohart- Adams  rate  relationships  are  of  a 
form  consistent  with  an  ideal  design  relationship  for  a  PFR.  An  inherent  assumption  with 
the  ideal  PFR  relationship  is  that  a  pseudo-steady-state  condition  exists  in  that  the 
conversion  is  spatially  dependent  (i.e.,  the  degree  of  conversion  is  a  function  of  axial 
location  only)  and  time  invariant.  An  ideal  PFR  can  be  likened  to  a  closed  batch  reactor 
(closed  in  the  sense  that  there  are  no  flow  streams  associated  with  the  reactor  vessel)  that 
is  conveyed  along  a  tube.  The  degree  of  conversion  in  such  a  system  is  dependent  upon 
the  space  (or  duration)  time  within  the  tube. 
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In  later  efforts  at  modeling  the  breakthrough  characteristics  of  carbon 
adsorbent  systems,  more  complex  rate  expressions  have  been  assumed.  Many  of  these 
rate  relationships  incorporate  parameters  to  differentiate  between  the  various  mechanisms 
that  contribute  to  the  removal  of  components  from  a  contaminated  airstream.  These 
removal  mechanisms  encompass  diffusive  considerations  as  well  as  reactionary  rate 
expressions.  The  analytics  associated  with  the  characterization  of  breakthrough  can 
become  formidable  even  with  simplistic  rate  mechanisms. 


6.  THE  MECKLENBURG  RELATIONSHIP 

A  conceptual  model,  attributed  to  Mecklenburg  (circa  1925),  is  based  on 
the  subtraction  of  a  “dead  zone  thickness”  from  the  actual  physical  length  of  an  adsorbent 
column  to  account  for  the  penetration  of  contaminants  far  below  the  challenge 
concentration  levels.  If  an  adsorbent  system  is  exposed  indefinitely  to  a  steady  vapor- 
phase  challenge  concentration,  the  effluent  concentration  will  eventually  rise  to  that  of 
the  inlet  challenge  concentration  as  the  adsorbent  capacity  is  diminished  with  time.  The 
assumptions  here  are  that  the  contaminant  is  irreversibly  adsorbed  onto  the  carbon 
substrate,  and  there  is  no  subsequent  off-gassing  of  the  contaminant  from  the  adsorbed 
state.  The  intent  of  this  section  is  to  review  the  basis  of  the  Mecklenburg  mass  balance. 


6.1  Graphical  Interpretations  for  a  Strongly  Adsorbed  Compound 

Figure  7  (ref  8)  depicts  an  adsorption  wave  profile  that  is  characteristic  of 
compounds  having  a  strong  affinity  for  the  carbon  adsorbent.  The  plot  represents  a 
snapshot  of  the  concentration  profile  throughout  the  bed  at  that  moment  in  time  in  which 
a  trace  or  threshold  breakthrough  concentration  of  contaminant  has  reached  the  effluent 
end  of  the  adsorbent  system.  In  most  cases  of  interest,  this  breakthrough  concentration  is 
several  orders  of  magnitude  less  than  the  influent  challenge  concentration. 
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Figure  7.  Depiction  of  adsorption  wave  profile  at  breakthrough  (from  ref  8). 


37 


A  review  of  Figure  7  indicates  two  distinct  regions  of  interest.  The  first 
region  of  interest  covers  the  adsorbent  depth  encompassed  within  the  region  bounded  by 
A  and  B.  In  this  region,  the  localized  vapor  phase  concentration  of  the  adsorbate  is  equal 
to  C0  (or  the  challenge  concentration).  This  region  represents  a  saturated  loading  region 
in  which  the  adsorbent  has  reached  its  maximum  loading  capacity.  Such  a  maximum  (for 
a  physically  adsorbed  contaminant)  would  be  bounded  by  equilibrium  constraints. 

The  second  region  of  interest  concerns  that  section  of  the  adsorbent  bed 
encompassed  within  the  region  bounded  by  B  and  C.  In  this  region,  the  adsorbent  retains 
adsorptive  capacity  because  it  has  not  reached  a  state  of  saturation,  and  the  concentration 
of  contaminant  in  the  vapor  phase  decreases  with  axial  position.  The  width  of  this  region 
is  due  to  the  finite  mass-transfer  resistances  that  exist  between  the  bulk  vapor  phase  and 
the  adsorbed  state.  These  resistances  affect  the  rate  of  adsorption  within  the  system. 

With  infinitely  rapid  mass  transfer  from  the  vapor  phase  to  an  adsorbed 
state,  the  downward  sloping  region  of  the  adsorption  wave  depicted  in  Figure  7  would  be 
a  vertical  line,  and  the  MTZ  thickness  would  be  zero.  The  characteristics  and  breadth  of 
the  decreasing  slope  that  constitutes  the  MTZ  are  associated  with  the  finite  rate  effects 
associated  with  the  adsorptive  process. 

The  time  required  for  a  breakthrough  concentration  to  be  observed  at  the 
effluent  end  (i.e.,  at  position  z  =  L)  of  the  adsorbent  bed  is  termed  as  the  breakthrough 
time  and  is  depicted  by  the  symbol  tb .  The  depth  of  the  adsorbent  corresponding  to  the 
region  of  Figure  7  bounded  by  B  and  C  is  termed  the  MTZ  thickness  and  is  depicted  by 
the  symbol  Mz.  The  breadth  and  characteristics  of  the  MTZ  region  of  the  adsorption  wave 
contour  are  rate  dependent. 


6.2  Development  of  the  Mecklenburg  Breakthrough  Relationship 

In  the  development  that  follows,  it  is  assumed  that  the  removal  mechanism 
is  due  solely  to  physical  adsorption,  and  that  the  adsorptive  process  is  irreversible. 

The  mass  of  contaminant  that  the  adsorbent  can  retain  at  full  saturation  is 
given  by  the  expression 


(A7/)adsorbed  @  saturation  —  Vtot  Pb  We 


(138) 


The  expression  given  by  eq  138  represents  an  equilibrium  bounding 
maximum  that  is  not  rate  dependent.  If  the  entire  volume  of  adsorbent  could  be  loaded  to 
its  saturation  capacity,  then  the  mass  of  adsorbate  contained  within  the  system  would  be 
described  by  the  relationship  given  by  eq  138  and  graphically  depicted  by  the  rectangular 
region  bounded  by  A-C-Ce-C0  with  reference  to  Figure  7.  The  total  volume  of  the 
adsorbent  system  includes  the  volume  of  the  adsorbent  as  well  as  the  associated  void 
volume.  The  volume  is  directly  correlated  to  the  depth  of  the  adsorbent  in  a  system 
having  a  constant  cross-sectional  area.  With  an  imposed  low-level  breakthrough 
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constraint,  the  entire  volume  of  the  adsorbent  cannot  be  used  to  a  condition  of  full 
saturation. 


The  mass  of  contaminant  that  is  introduced  into  the  adsorbent  system  as  a 
vapor  can  be  expressed  by 


(A/fjchallenge  —  Q  C0  t 


(139) 


If  it  is  then  assumed  that  there  is  a  zone  or  region  within  the  carbon 
adsorbent  bed  that  cannot  be  fully  used  for  capacity  due  to  a  breakthrough  constraint, 
then  the  total  volume  of  the  adsorbent  bed  must  be  adjusted  (reduced)  to  account  for  this. 
This  breakthrough  constraint  (concentration  based)  gives  rise  to  an  adjustment  factor  to 
the  total  volume  appearing  in  eq  138.  Using  eq  138  as  a  basis,  a  suggested  relationship  to 
describe  the  mass  of  contaminant  retained  by  the  adsorbent  at  a  breakthrough  condition  is 
given  as 


(iUrjadsorbed  @  breakthrough  —  8  Vtot  P b  We  (140) 

where  5  represents  that  fraction  of  the  total  adsorbent  system  volume  that  can  be  used  for 
loading  until  a  breakthrough  criteria  is  reached. 

The  mass  of  contaminant  that  enters  the  system  at  that  specific  time  at 
which  a  breakthrough  criteria  is  reached  at  the  effluent  can  be  simply  expressed  through 
the  following  modification  to  eq  139: 


(A//)challenge  @  breakthrough  —  QCotb  (141) 

Equating  the  mass  of  contaminant  introduced  into  the  system  eq  141  to  the 
mass  of  contaminant  retained  by  the  adsorbent  at  the  breakthrough  criteria  (eq  140)  yields 
the  following  relationship: 


Q  C0  lh  =  8  Vtot  p  b  We 
or 

tb  =  Kpb  We)/(Q  C0)\  8  Vlol  (142) 

Equation  142  represents  a  very  simplified  mass  balance  expression 
developed  through  the  imposition  of  a  breakthrough  criteria.  An  inherent  assumption  to 
this  development  is  that  the  adsorbent  has  a  characteristic  and  uniform  equilibrium 
loading  capacity  throughout  the  bed. 

If  the  cross-sectional  area  normal  to  the  principle  flow  axis  is  a  constant, 
then  Vtot  =  A  L,  and  eq  142  can  be  expressed  as 

tb  =  Kpb  we)/(Q  C0)\  Sal  (143) 
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A  further  modification  to  eq  143  can  be  made  by  noting  that  the 
volumetric  throughput  can  be  related  to  the  cross-sectional  flow  area  and  the  superficial 
flow  velocity  through  Q  =A  u.  Substitution  of  this  relationship  into  eq  143  yields  the 
expression 


tb  =  Kp/,  We)/(u  C(J)\  8  L  (144) 

The  various  forms  of  eq  142  represent  variations  on  a  simplified  mass 
balance  concept. 

It  is  hypothesized  that  Mecklenberg  related  the  fraction  of  the  usable  bed 
depth  in  a  manner  consistent  with  the  following  expression: 

5  L  =  L-MZ  (145) 

where  Mz  represents  that  thickness  (or  segment)  of  the  total  bed  depth  that  cannot  be  fully 
used  for  adsorption  capacity.  Substitution  of  eq  145  into  eq  144  yields 

tb  =  [(p*  We)/(u  Co)]  (L  -  M-)  (146) 

Equation  146  is  consistent  with  the  analytical  approaches  that  are  currently 
used  to  estimate  the  performance  characteristics  of  carbon-based  filter  systems  in  military 
applications  against  contaminants  that  are  strongly  adsorbed  onto  the  adsorbent  through  a 
physical  adsorption  mechanism.  Several  parameters  within  eq  146  can  be  readily 
evaluated.  The  difficulty  in  the  application  of  eq  146  rests  with  the  evaluation  of  M-  and, 
to  a  lesser  extent,  We. 

It  is  surmised  that  Mecklenburg  pursued  an  analogy  similar  to  this  in  the 
mass  balance  relationships  that  are  attributed  to  him.  These  simplified  relationships  are 
thought  to  be  consistent  with  the  various  forms  of  eq  142. 


7.  EFFORTS  BY  I.  M.  KLOTZ  AND  M.  DOLE 

7.1  Background 

Within  the  open  technical  literature,  numerous  articles  relate  to  the 
adsorption  of  contaminant  compounds  on  porous  solids.  This  effort  will  focus  on  two 
such  efforts  because  of  their  relevance  to  the  current  review. 

Irving  M.  Klotz  (and  other  investigators)  authored  several  technical 
articles  relating  to  the  analytical  aspects  of  estimating  the  performance  characteristics  of 
activated  carbon  against  vapor-phase  contaminants.  Over  the  years,  others  have  also 
contributed  and  further  developed  the  analytical  approaches  that  are  used.  The  Klotz 
(ref  2)  and  Dole  and  Klotz  (ref  3)  efforts  relate  to  the  removal  of  compounds  of  military 
significance  from  a  contaminated  airstream. 
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7.2 


“The  Adsorption  Wave”  by  I.  M.  Klotz 


“The  Adsorption  Wave”  (1946;  ref  2)  constituted  a  chapter  within  a 
compiled  technical  report  summary  generated  through  the  National  Defense  Research 
Committee  (NDRC)  during  the  1940s.  This  early  work  describes  the  development  of  a 
differential  mass  balance  and  several  of  the  rate  aspects  that  factor  into  the  evaluation  of 
the  resultant  analytical  relationships.  A  modified  version  of  the  NDRC  chapter  by  Klotz 
was  also  published  in  the  technical  journal  Chemical  Reviews  (ref  9).  The  Chemical 
Reviews  version  was  intended  for  general  public  release  and  does  not  contain  the 
compound  specificity  that  appears  in  the  NDRC  version. 

As  indicated  in  Sections  2  and  4  of  this  current  review  effort,  the  analytics 
associated  with  adsorption  are  complicated.  It  is  for  this  reason  that  rigorous  approaches 
to  characterize  the  performance  of  activated  carbon  adsorbent  substrates  against 
contaminated  vapor  streams  are  often  dealt  with  through  simplified  analytical 
expressions.  Some  analytical  modeling  approaches  are  very  dependent  on  actual  data  to 
estimate  parameter  values,  whereas  other  modeling  approaches  are  less  dependent  on 
actual  data  and  make  use  of  correlative  estimation  methods.  In  the  absence  of  actual  data, 
correlative  approaches  can  provide  reasonable  tools  to  estimate  filter  performance. 
Although  it  is  useful  to  appreciate  the  foundations  upon  which  the  analytical  models  are 
based,  the  resultant  relationships  are  often  so  complex  as  to  render  them  impractical  for 
estimation  purposes  by  the  generic  engineering  design  community. 

In  “The  Adsorption  Wave”  (ref  2),  Klotz  provides  a  relationship  based  on 
a  specific  rate-limiting  step  that  has  considerable  utility  to  the  design  community.  The 
resultant  expression  requires  a  minimum  of  auxiliary  data  generation  and  is  commonly 
employed  to  this  day.  On  the  basis  of  the  simplifications  that  went  into  the  developed 
expression,  there  are  limitations  inherent  with  its  applicability  for  all  situations.  These 
inherent  limitations  require  understanding  in  the  application  of  the  modeling  expressions. 


7.2.1  Fundamental  Analytical  Relationships  from  “The  Adsorption  Wave” 


The  Klotz  effort  begins  with  the  formulation  of  a  mass  balance  across  a 
volume  element  of  infinitesimal  thickness.  The  analytical  development  is  based  on  the 
specification  of  a  parameter  q,  which  represents  the  amount  of  toxic  gas  within  the 
volume  element  per  unit  volume  of  bed.  This  q  parameter  is  a  localized  variable  for  the 
amount  of  adsorbate  that  exists  in  both  the  adsorbed  state  and  the  vapor  phase.  In  effect, 
q  represents  a  combination  of  C,  (i.e.,  the  concentration  of  adsorbate  in  the  vapor  phase) 
and  q,  (i.e.,  the  adsorbate  loading  on/within  the  carbon  granules). 

The  mathematical  interpretation  of  q  is  given  as 


q  =  (moles  of  toxic  gas)  /  (volume  of  bed) 


=  Mi  /  [(MW);  Vtot] 


(147) 
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where  Vtot  represents  the  combined  volumes  of  the  adsorbent  and  void  volumes.  Vtot  does 
not  represent  the  total  volume  of  the  adsorbent  column. 

The  amount  and  distribution  of  toxic  gas  (between  the  adsorbed  solid 
phase  and  the  vapor  phase)  is  a  function  of  time  and  position.  Relative  to  a  fixed  spatial 
reference,  the  localized  rate  of  pickup  per  unit  volume  of  bed  can  be  represented  by 
differentiating  eq  147  with  respect  to  time.  This  gives  the  relationship 

(dil Idi)  =  {  1/[(MW);  Vtot] }  (dMi/dt)  (148) 

In  eq  148,  it  has  been  assumed  that  [5(Az)/dt]  is  negligible  relative  to 
(dMj/dt)  for  a  system  having  a  constant  cross-sectional  flow  area. 

With  the  localized  rate  of  change  of  mass  established  by  eq  148,  the 
development  of  a  mass  balance  over  a  differential  volume  element  proceeds  as  follows. 

Consider  the  adsorption  system  depicted  in  Figure  8. 


Figure  8.  A  cylindrical  adsorbent  system  with  primary  flow  in  the 
axial  direction  (from  ref  10). 


Figure  8  depicts  a  cylindrical  adsorbent  system  in  which  the  principle  flow 
is  associated  with  the  axial  direction.  The  figure  illustrates  a  sectional  view  of  the 
differential  volume  element  in  a  cylindrical  configuration  and  is  consistent  with  Figure  2. 
The  gradients  in  the  angular  and  radial  directions  are  neglected  in  the  development  which 
follows. 

7.2. 1.1  Rate  of  Mass  Flow  of  i  into  the  Volume  Element  at  Axial  Position  z: 

Wi)  in 


(tlli)m  =  (MW),-  N,{z,  t)  (8  A) 


(149) 


The  molar  flux  appearing  in  eq  149  is  evaluated  at  axial  position  z  and  time  1. 
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7.2. 1.2 


Rate  of  Mass  Flow  of  i  Exiting  Volume  Element  at  Axial  Position 
(z  +  Az):  (TRi) out 


(771,) out  =  (MW),-  Ni(z  +  A z,t  +  At)  (8  A)  (150) 


As  with  eq  149,  the  molar  flux  of  i  exiting  the  differential  volume  element 
is  evaluated  at  axial  position  (z  +  Az)  and  time  (t  +  At). 

7.2. 1 .3  Rate  of  Change  of  Mass  of  i  within  the  Volume  Element:  ( dMJdt ) 

An  expression  for  the  localized  rate  of  change  of  mass  within  the  volume 
element  can  be  obtained  from  the  relationship  given  by  eq  148.  From  eq  148,  the 
following  relationship  can  be  written: 

(dM-Jdt)  =  (MW),-  Vtot  (or\/dt)  (151) 

7.2. 1.4  Analytical  Mass  Balance  Development 


In  the  development  that  follows,  the  rate  of  consumption  of  species  i 
through  a  reactionary  means  is  included  within  the  q  parameter.  As  such,  the  following 
mass  balance  expression  for  i  can  be  written: 


(dMJdt)  =  (m,)  in 


(152) 


where  ( 1Tl)gen  =  0  because  the  rate  of  removal  of  i  through  a  reactionary  pathway  is 
reflected  in  the  localized  adsorptive  flux  [i.e.,  (dx\!dt)\. 


Substituting  the  relationships  given  by  eqs  149,  150,  and  151  into  eq  152 
yields  the  following  expression: 

Vtot  (dr\/dt)  =  -£  A  [NJz  +  Az,  t  +  At)  -  NJz,  t)\  (153) 


The  quantity  [N,(z  +  A z,t  +  At)  -  NJz,  t)\  represents  the  total  differential  of  the  molar  flux 
of  i  (see  Appendix  B).  The  relationship  given  by  eq  153  can  therefore  be  written  as 


Vtot  (dx\!dt)  =  -8  A  dV,  (154) 

Now,  V0t  represents  the  volume  occupied  by  the  void  spaces  and  the 
volume  occupied  by  the  adsorbent  within  the  differential  volume  element.  Assuming  a 
constant  cross-sectional  area  normal  to  the  principle  flow  direction,  the  relationship 
Vtot  =  A  dz  holds.  Substituting  this  relationship  into  eq  154  gives  the  expression 


(dr\/dt)  =  -8  (dV,/dz) 


(155) 
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Relationships  for  the  molar  flux  and  gradient  in  the  axial  direction  have 
been  used  throughout  this  review.  For  an  adsorbent-based  system,  the  following 
relationship  (eq  91)  is  applicable: 

Ni  =  vz  Ci  -  Bl  ( dCi/dz ) 

Klotz  assumes  that  the  dispersive  effects  associated  with  nonideal  plug 
flow  are  considered  to  be  negligible.  To  this  end,  D l  =  0,  and  the  molar  flux  can  be 
expressed  by 

Ni  =  V;  Ci 

The  total  derivative  of  this  expression  with  respect  to  the  axial  direction  yields 

(d/V,/dz)  =  vz  (dC/dz)  (156) 

The  interstitial  velocity  is  taken  to  be  constant  in  eq  156. 

Substituting  the  relationship  given  by  eq  156  into  eq  155  yields 

-(1/8)  (dr\/dt)  =  vz  (dCi/dz)  (157) 

The  concentration  gradient  of  eq  157  represents  the  total  gradient  of  C, 
with  respect  to  axial  position.  This  derivative  is  dealt  with  as  follows  (see  Appendix  B  for 
information  relating  to  total  and  partial  derivatives). 

Because  the  vapor-phase  concentration  is  a  function  of  axial  location  and 
time,  the  following  functionality  can  be  written: 

Q  =  Cj  (z,  t) 

The  total  differential  of  this  relationship  can  be  expressed  as 

d Ci  =  (dCi/dz)  d z  +  (i dCJdt )  d t  (158) 

or 

(dC,/dz)  =  (dCM)  +  ( dCJdt )  (dt/dz)  (159) 

Because  (dt/dz)  is  taken  to  be  the  reciprocal  of  the  interstitial  velocity,  eq 
159  can  be  expressed  as 

(dC,/dz)  =  ( dCi/dz )  +  ( l/vz)  ( GCJdt )  ( 1 60) 

Substituting  the  relationship  given  by  eq  160  into  eq  157  gives 

-(1/8)  (dr\/dt)  =  (dCj/dt)  +  vz  ( dC-Jdz )  (161) 
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Equation  161  is  the  same  as  eq  1 1  in  the  Klotz  reference  paper  (ref  2). 

It  is  instructive  at  this  point  to  compare  the  Klotz  mass  balance  expression, 
as  given  by  eq  161,  to  the  relationship  that  appears  in  the  Ruthven  text  (ref  6).  A  form  of 
this  Ruthven  reference  relationship  was  provided  by  eq  86. 

The  intent  of  the  following  development  is  to  demonstrate  the  consistency 
between  the  Klotz  mass  balance  and  that  given  in  the  text  by  Ruthven.  The  Ruthven 
(ref  6)  relationship  is  generically  applicable  to  systems  in  which  axially  dispersed  plug 
flow  can  be  assumed. 

A  relationship  between  the  total  volume  of  the  adsorbent  bed  and  the  bulk 
packing  density  of  the  adsorbent  was  reviewed  in  the  development  of  the  relationship 
leading  to  eq  61.  From  eq  61,  the  following  relationship  can  be  written: 

V,o,  =  Mads/p,  (162) 

Substituting  the  relationship  given  by  eq  162  into  eq  147  gives 


q  =  M,/[(MW)i  (AWpzdJ  =  (M,/Mads)  [p,,/(MW),]  (163) 

Now  (M,/Mads)  is  the  specific  loading  of  species  i  on  the  adsorbent  expressed  on  a  mass 
basis,  which  was  previously  defined  as  Wj.  Equation  163  can  therefore  be  equivalently 
expressed  as 


q  =  Wj  [p//(MW)i] 


(164) 


Differentiation  of  eq  164  with  respect  to  time  yields 

(8r\/dt)  =  ( p/A MW),]  (GWJGt)  (165) 

A  relationship  relating  the  mass  based  adsorbate  loading  to  a  volume- 
based  loading  was  previously  developed  as  eq  69,  which  is  given  as 

(dWJdt)  =  [(1  -  £)/ Pb]  ( dqjdt ) 

Substituting  the  relationship  given  by  eq  69  into  eq  165  yields 

(8f]/dt)  =  [(1  -  £)/(MW)i]  (dqjdt)  ( 1 66) 

Substituting  the  relationship  given  by  eq  166  into  eq  161  (which 
represents  eq  11  of  the  Klotz  mass  balance  relationship)  gives  the  following  expression: 

y-  (GCJGz)  +  (dCJdt)  +  [  l/(MW)j]  [( 1  -  s)/s]  (dqjdt)  =  0  (167) 
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Neglecting  the  effects  associated  with  dispersion  (Dt  =  0)  and  assuming 
that  the  interstitial  velocity  in  the  axial  direction  is  constant  (as  Klotz  does),  then  eq  167 
is  observed  to  be  consistent  with  the  Ruthven  relationship  as  given  by  eq  86. 


7.2.2  The  Ruthven  Relationship 

Throughout  this  review  effort,  comparisons  to  a  relationship  appearing  in 
the  Ruthven  text  reference  (ref  6)  have  been  made.  Using  the  approaches  contained  in 
Section  4,  the  development  of  the  Ruthven  relationship  is  briefly  summarized  as  follows. 

7.2.2. 1  Rate  of  Mass  Flow  of  i  into  the  Volume  Element  at  Axial  Position  z: 

Wi)  in 


(7ft()in  =  (MW),  V,(z)  (8  A) 


(168) 


7. 2. 2. 2  Rate  of  Mass  Flow  of  i  Exiting  the  Volume  Element  at  Axial 

Position  (z  +  A z):  (772,)0u t 


(m()out  =  (MW),  N,(z  +  Az)  (8  A) 


(169) 


7 .2.2.3  Rate  of  Change  of  Mass  of  i  within  the  Volume  Element:  ( dMJdt ) 

A  variation  to  the  approach  used  to  develop  eq  76  is  needed  here  because 
of  the  difference  in  the  geometry  of  the  system  being  considered.  In  this  current 
development,  a  cylindrical  system  is  again  considered;  however,  the  principle  flow 
direction  is  now  taken  to  be  in  the  axial  direction  (with  eq  76,  the  primary  flow  direction 
was  assumed  to  be  radially  directed). 

In  accordance  with  the  development  of  eq  76,  the  following  relationships 
describe  the  localized  mass  of  i  within  the  element  depicted  in  Figure  8  (assuming  that 
the  property  characteristics  are  uniform  within  the  differential  slice): 


Mj  —  (A/,')  adsorbed  phase  (A7/)  vapor  phase 

=  Vtot  (fads/btot)  Qi  +  [(Af,')vap/IAoid]  lAoid 

=  [A  Az]  (1  -  8)  qi  +  (MW),  C,  (8  Vtot) 

=  [A  Az]  (1  -  8)  q,  +  (MW),  C,  8  (A  Az) 


(170) 
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Differentiating  eq  170  with  respect  to  time,  assuming  that  the  void  fraction 
and  the  area  normal  to  the  principle  flow  direction  can  be  taken  as  constants,  and 
assuming  that  the  terms  [5(A z)ldt]  are  negligible  yields 

(dMj/dt)  =  [A  Az\  (1  -  8)  {dqjdt)  +  (MW),  8  (A  A z)  (, dQ/dt )  (171) 


7. 2. 2.4  Species  i  Mass  Balance 

The  generic  mass  balance  for  species  i  is  given  by  eq  152.  Substitution  of 
the  relationships  given  by  eqs  168,  169,  and  171  into  eq  152  yields 

[  1/(MW),]  [( 1  -  S)/S]  {dq-Jdi)  +  (dQ/dt)  =  -(, dNJdz )  ( 1 72) 

The  axial  flux  can  be  expressed  in  terms  of  its  convective  and  diffusional 
components  by  the  following  modified  form  of  eq  39  (modified  in  that  an  interstitial 
velocity  is  used  as  opposed  to  a  superficial  velocity): 

(dNi/dz)  =  [d(v  Ci)/8z ]  -  Bl  ( d2Ci/dz 2)  (173) 

Substitution  of  the  relationship  given  by  eq  173  into  eq  172  yields 

-Bl  {oCJoz  )  +  [S(v  Ci)/dz\  +  (dQ/dt)  +  [1/(MW)J  [(1  -  s)/8]  (oqjdt)  =  0  (174) 

Equation  174  is  seen  to  be  identical  to  eq  86,  which  represents  the 
Ruthven-based  reference  mass  balance  relationship. 


7.2.3  Aspects  Associated  with  (dr\/dt)  Appearing  in  the  Klotz  Article  (Ref  2) 

The  dynamics  associated  with  the  adsorptive  loading  within  the  MTZ  of 
an  activated  carbon  column  factors  considerably  into  the  form  of  the  analytical 
expressions  that  would  characterize  the  vapor-phase  concentration  of  the  contaminant  as 
a  function  of  time  and  position.  The  localized  sorptive  flux,  expressed  through  ( dr\/dt )  or 
(dqi/dt)  (or  other  equivalent  relationships),  is  very  much  rate-dependent  and  will 
influence  the  forms  of  the  concentration  expressions. 

The  mechanistic  aspects  associated  with  the  adsorption  process  include  the 

following: 


a.  The  transfer  of  mass  from  the  bulk  fluid  phase  to  the  outer  surface  of 
the  carbon  granules.  This  mechanism  is  referred  to  as  bulk  diffusion. 
In  the  bulk  vapor  phase,  the  contaminant  molecules  are  largely 
interacting  with  diluent  molecules  when  the  contaminant 
concentrations  are  small.  As  the  contaminant  molecules  diffuse 
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through  the  boundary  layer,  their  interactions  among  themselves 
increase  as  their  localized  concentrations  increase  with  decreasing 
distance  from  the  adsorbent  surface.  The  rate  of  transfer  is 
proportional  to  a  transfer  coefficient  and  a  driving  force  gradient. 

b.  The  intraparticle  transfer  of  the  adsorbed  species  into  the  pore 
structure  of  the  individual  carbon  granules.  This  mechanism  is  called 
intraparticle  or  Knudsen  diffusion.  At  this  level,  the  contaminant- 
contaminant  and  contaminant- surface  interactions  occur  at  a  much 
greater  frequency  than  would  have  occurred  as  the  contaminant 
diffused  through  the  boundary  layer.  The  localized  concentration  of 
the  contaminant  increases  significantly  as  the  molecules  are 
transported  within  the  porous  networks  of  the  carbon.  The  rate  of  this 
transport  mechanism  would  be  proportional  to  a  transfer  coefficient  as 
well  as  a  driving  force. 

c.  For  situations  in  which  the  contaminant  removal  is  through  a 
chemisorption  pathway,  a  reactionary  rate  expression  may  very  well 
factor.  Such  rates  are  typically  proportional  to  a  reaction  rate  constant 
as  well  as  concentration(s). 

d.  Other  potential  rate-influencing  factors. 

It  becomes  evident  that  the  analytical  complexity  of  the  Klotz  mass 
balance  expression,  as  given  by  eq  161,  can  become  formidable  with  even  simplified 
forms  of  (dr\ldt).  To  limit  the  complexity  of  the  analytic  rate  relationships,  the  general 
approach  is  to  assume  that  a  singular  mechanistic  pathway  governs  the  overall  removal 
process.  This  singular  mechanistic  pathway  is  referred  to  as  the  rate-limiting  step. 

Compounds  having  a  strong  affinity  for  the  surface  of  the  activated  carbon 
substrate  tend  to  exhibit  strong  physical  adsorption  characteristics.  The  adsorption  wave 
contours  of  this  class  of  materials  tend  to  be  consistent  with  those  appearing  in  Figures  5 
and  7.  With  reasonable  flow  velocities,  the  rate-limiting  step  governing  the  removal  of 
this  class  of  materials  from  a  contaminated  airstream  is  assumed  to  be  bulk  diffusion. 

The  Klotz  article  associates  a  critical  bed  thickness  (termed  /*)  to  each  of 
the  mechanistic  pathways  associated  with  a  particular  removal  process.  Correlative 
relationships  for  the  primary  mechanistic  pathways  are  provided  in  the  article.  These 
critical  bed  contributions  have  an  analogy  to  the  heights  of  transfer  units  in  the  chemical 
engineering  practice  relating  to  mass  transfer  operations.  A  difficulty  in  using  these 
expressions  rests  with  parameter  quantification. 

For  cases  in  which  bulk  diffusion  constitutes  the  predominant  removal 
mechanism,  Klotz  suggests  the  use  of  a  correlation  attributed  to  Gamson,  Thodos,  and 
Hougen  to  calculate  the  critical  bed  depth  (M-).  The  Gamson,  Thodos,  and  Hougen 
(GTH)  correlation  presented  in  the  Klotz  reference  (ref  2)  is  given  as 
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I,  =  (2303/a)  (dp  G/p)0'41  [p/(p  B)f67  log  (CJCe) 


(175) 


in  which  the  dimensionless  grouping  (dp  Gt p)  is  the  Reynolds  number  (Re)  and  the 
dimensionless  grouping  [p/(p  B)]  represents  the  Schmidt  number  (Sc). 

It  is  indicated  that  the  GTH  correlation  was  derived  for  cases  in  which  the 
ratio  of  ( C0/Ce )  is  very  large  and  is  applicable  only  in  cases  of  turbulent  flow.  Klotz 
presents  a  plot  in  which  the  airflow  resistance  at  various  flow  rates  through  carbon 
columns  of  varying  particle  sizes  was  measured.  It  was  found  that  there  was  an 
approximate  transition  from  laminar  to  turbulent  flow  at  Re  =  10  for  the  charcoals 
investigated  (ref  2).  The  transition  region  was  identified  as  that  region  in  which  there  was 
a  notable  change  in  the  slope  of  the  airflow  resistance  versus  flow-rate  plot.  It  is  further 
indicated  that  although  the  correlation  given  by  eq  175  is  more  appropriately  applicable 
to  turbulent  flow  systems,  it  can  also  satisfactorily  approximate  It  in  the  laminar  flow 
regimes. 

7.2.4  Breakthrough  Time  Relationship 

In  situations  in  which  diffusion  in  air  constitutes  the  predominant  rate¬ 
controlling  step,  Klotz  combined  the  simplified  mass  balance  relationship  attributed  to 
Mecklenburg  (eq  146)  with  the  correlation  given  by  eq  175  to  arrive  at  the  following 
breakthrough  time  relationship: 

tb=[(pbWe)/(uC0)](L-Mz) 

=  [(p*  We)/(u  C0)]  {  L  -  { (1/a)  (dp  G/p)0'41  [p/(p  B)]a67  In  (CJCe) }  }  (176) 

=  [(p*  We)/(u  Co)]  {L-  [(1/a)  Re041  Sc067  In  (CJCe)]  } 

The  relationship  given  by  eq  176  forms  the  basis  of  current  approaches  to 
estimate  the  breakthrough  times  for  compounds  that  are  strongly  adsorbed.  Current 
relationships  incorporate  additional  parameters  aimed  at  refining  the  analytical  results  to 
actual  data;  however,  the  fundamental  underlying  principles  that  went  into  eq  176  are 
assumed  to  be  valid  to  this  day.  The  use  of  this  expression  requires  an  understanding  of 
the  underlying  assumptions  that  went  into  its  development.  These  assumptions,  which 
appear  throughout  this  report,  are  compiled  as  follows  for  convenience: 

a.  Bulk  diffusion  is  the  predominant  rate-controlling  step  relating  to  the 
removal  of  contaminants  from  the  bulk  vapor  phase  to  the  carbon 
surface.  Intraparticle  mass  transfer  and  reactionary  processes  (if 
applicable)  are  assumed  to  occur  on  a  timescale  that  is  very  rapid 
relative  to  the  rate  of  bulk  phase  mass  transfer. 


49 


b.  The  contaminant  is  strongly  adsorbed  onto  the  surface  of  the  adsorbent 
through  an  ideal  physical  adsorption  pathway.  The  adsorption  is 
irreversible. 

c.  Propagation  of  the  adsorption  wavefront  is  consistent  with  the 
trendlines  appearing  in  Figures  5  and  7.  The  adsorption  waves  retain  a 
consistent  contour  with  the  progression  of  time.  This  is  inferring  a 
constant  value  of  Mz  (or  I, )  with  time.  The  adsorbent  bed  is  assumed  to 
have  sufficient  depth  to  permit  the  establishment  of  a  uniform  wave 
profile. 

d.  The  ratio  of  the  challenge  concentration  to  the  breakthrough 
concentration  is  very  large.  Klotz  (ref  2)  does  not  clarify  as  to  what 
constitutes  “very  large”. 

e.  The  inlet  challenge  concentration  is  constant,  and  the  equilibrium 
loading  capacity  is  a  constant.  This  assumption  is  consistent  with  that 
proposed  by  Bohart  and  Adams  (ref  7)  with  their  residual  loading 
parameter  (a). 

f.  The  flow  behavior  is  ideal  such  that  the  effects  associated  with  axial 
dispersion  can  be  neglected. 

7.3  Sorption  of  Chloropicrin  and  Phosgene  on  Charcoal  from  a  Flowing  Gas 

Stream  by  M.  Dole  and  I.  M.  Klotz  (1946;  ref  3) 

In  this  work.  Dole  and  Klotz  conducted  an  extensive  investigation  into 
determining  how  various  parameters  influenced  the  breakthrough  characteristics  of 
chloropicrin  (PS)  and  phosgene  (CG)  within  carbon  beds.  PS  has  a  moderate  vapor 
pressure  (26  torr  at  25  C),  and  its  removal  mechanism  on  activated  carbon  is  attributed  to 
moderate  physical  adsorption  (ref  11).  Phosgene  has  a  much  higher  room-temperature 
vapor  pressure  (1425  torr  at  25  C)  than  PS;  its  removal  mechanism  is  attributed  to  a 
combination  of  hydrolysis  and  subsequent  reaction  on  impregnated  activated  carbons 
(ref  11). 


The  parameters  investigated  by  Dole  and  Klotz  included  influent 
concentration  and  temperature,  relative  humidity,  bed  depth,  and  flow  rates.  In  the 
analysis  that  follows,  the  analytical  approaches  pursued  by  Dole  and  Klotz  (ref  3)  in 
quantifying  the  critical  bed-depth  thickness  (or  M-)  are  reviewed. 

The  Dole  and  Klotz  (ref  3)  analytical  approach  can  be  likened  to  a  further 
development  of  the  rate  expressions  suggested  by  Bohart  and  Adams,  with  their  early 
efforts  relating  to  the  breakthrough  of  chlorine  through  carbon  adsorbents.  In  the  Dole 
and  Klotz  effort,  the  rate  expressions  of  Bohart  and  Adams  are  combined  with  the 
Mecklenburg  mass  balance  concept  using  modem  forms  of  mass  transfer  correlative 
relationships. 
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A  review  of  the  effort  of  Dole  and  Klotz  (ref  3)  is  summarized  as  follows. 
The  Mecklenburg  relationship  as  given  by  Dole  and  Klotz  is 


C„  L  tb  =  A  n0  (x  -  h)  (177) 

where  x  is  the  total  bed  length,  L  is  the  volumetric  flow  rate,  C0  is  the  inlet  concentration 
expressed  on  a  molar  basis,  and  n0  is  the  equilibrium  saturation  capacity  of  the  carbon 
expressed  on  a  moles  of  i  per  volume  of  adsorbent  basis.  The  equilibrium  saturation 
capacity,  n0,  is  related  to  We  by  the  following  correspondence  through  dimensional 
considerations: 


We  =  na  (MW);  [(1  -  8)/  pft]  (178) 

Substituting  the  relationship  given  by  eq  178  into  eq  177  and  using 
variable  designations  that  have  been  consistently  used  throughout  this  review  gives  the 
following  expression: 

(MW);  Ca  Q  tb  =  AWe  [p*/(  1  -  8)]  (L  -  h)  (179) 


7.3.1  Linearity  of  Life-Thickness  Plots 

Figure  9  presents  a  life-thickness  plot  for  phosgene  at  various  flow  rates. 


Figure  9.  Life-thickness  curves  for  phosgene  at  various  flow  rates  (from  ref  3). 


The  plot  suggests  a  linear  relationship  between  breakthrough  time  and  bed 
thickness.  On  extrapolation  to  tb  =  0,  the  linear  trend  lines  cross  the  abscissa  at  a  finite 
(i.e.,  nonzero)  thickness.  This  thickness  represents  the  critical  bed  depth.  Dole  and  Klotz 
describe  this  thickness,  x0,  as  that  depth  of  adsorbent  just  sufficient  to  prevent  the 
penetration  of  gas  at  a  concentration  greater  than  that  established  for  the  breakthrough 
concentration  at  zero  time.  An  analogy  is  then  made  in  that  this  critical  bed  thickness, 
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x0  =  Zo,  corresponds  to  the  h  parameter  given  in  the  Mecklenburg  relationships  (eqs  177 
and  179). 


7.3.2  Analytical  Basis  for  the  Dole  and  Klotz  Development 

The  analytical  development  pursued  by  Dole  and  Klotz  continues  with  a 
relationship  appearing  in  the  Bohart  and  Adams  technical  paper  (ref  7).  This  relationship 
(eq  19  of  ref  7)  is  given  as 

(l/C/)=  l-e“J'  +  e(z'“°  (180) 


It  is  noted  that  eq  180  is  the  same  as  eq  124,  which  was  considered  in 
Section  5.2  of  this  review.  Substituting  the  appropriate  variable  transformations  (from 
Bohart  and  Adams)  into  eq  180  yields 

(CJCd  -  1  =  e[" {zJv)  ~kCot]_Q-kCot  ( 1 8 1 ) 


Equation  181  can  be  reexpressed  through  the  following  series  of 


expressions: 


(CJCd  -  1  =  [e 


_  ao  ( zJv )  j  k  Co  t 


']  -(\lckCot) 


YtKavy7JVI-\\l  t 


Jc  ao  (z/v) 

_  [e*ao(z/v)_  X]  e 


k  Co  t 

~k  Co  t 


(182) 


Taking  the  natural  logarithms  of  each  side  of  eq  182  yields  the  following 
series  of  expressions: 


In  [(CJCd  -  1]  =  In  { [ckao  (zhj)  -  1]  eTk  Co ' } 

=  In  [ekao(zlv)-  \\  +  \nekCot  (183) 

=  In  [ekao(z,v)-l]-kC0t 

Bohart  and  Adams  describe  the  variable  a0  as  an  initial  uniform  volume 
capacity  of  the  adsorbent  for  the  adsorbing  component.  In  other  words,  aQ  can  be  likened 
to  an  equilibrium  saturation  capacity  ( na ).  Equating  a„  to  n„,  eq  183  can  be  written  as 
follows: 


In  [(CJCd  -  1  ]  =  In  [ek  no  {zJu)  -  1  ]  -  k  Ca  t  (184) 


Equation  184  is  the  same  as  eq  2  in  the  Dole  and  Klotz  paper  (ref  3). 


7.3.3  Critical  Bed  Depth  Proportional  to  In  (CJCJ 

At  zero  breakthrough  time,  eq  184  gives 
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(CJCj)  -  1  =  e 


1 


k  no  (zo/v) 

or 

z0  =  [v/(k  n0)\  In  (CJCh)  (185) 

Through  this  analysis,  having  a  basis  with  the  Bohart  and  Adams 
development,  the  critical  bed  thickness  ( z0 )  is  observed  to  be  proportional  to  In  ( CJCb ). 

7.3.4  Correlative  Mass  Transfer  Relationships 

In  this  section,  mass  transfer  rates,  incorporating  a  mass  transfer 
coefficient,  are  used  to  model  the  adsorption  process  through  a  diffusion-limiting,  rate¬ 
controlling  mechanism.  Dole  and  Klotz  make  use  of  a  correlative  relationship  involving 
the  dimensionless  Reynolds  and  Schmidt  numbers  in  their  analytical  development.  The 
approach  is  to  equate  the  rate  of  adsorption  onto  the  surface  of  the  adsorbent  within  the 
active  MTZ  to  the  rate  of  introduction  of  contaminant  at  the  influent.  Their  analytical 
development  is  detailed  as  follows. 

The  rate  at  which  i  is  transferring  onto  the  adsorbent  is  given  by 

Wi) ads  =  (MW),  kg  a  A  M-  A pm  (186) 

where  kg  represents  a  localized  mass  transfer  coefficient  expressed  on  a  molar  basis. 
Examination  of  eq  186  reveals  its  similarities  to  modem  mass  transfer  expressions  in 
which  the  rate  of  mass  transfer  is  proportional  to  the  product  of  a  transfer  coefficient  and 
a  gradient,  which  serves  as  the  driving  force.  The  relationship  given  by  eq  186 
corresponds  to  eq  6  of  the  Dole  and  Klotz  reference  (ref  3).  A  narrative  description  of  the 
parameter  groupings  given  by  eq  186  is  as  follows: 

a:  The  effective  area  of  adsorbent  (for  mass  transfer)  per  unit 

volume  of  adsorbent. 

A  Mz :  The  volume  within  the  adsorbent  column  in  which  mass 

transfer  occurs.  This  is  also  the  volume  that  comprises  the 
MTZ  thickness. 

a  A  AT:  The  effective  area  within  the  MTZ  region  available  for  mass 

transfer.  This  represents  the  area  available  for  bulk  transport  of 
i  from  the  vapor  phase  to  the  exterior  surface  of  the  adsorbent 
granules. 

kg.  The  localized  gas-phase  mass  transfer  coefficient  expressed  on 

a  molar  basis. 

(MW),  kg\  The  localized  gas-phase  mass  transfer  coefficient  expressed  on 
a  mass  basis. 

A pm:  The  pressure  gradient  driving  force. 
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Dole  and  Klotz  give  the  following  correlative  relationship  for  the  mass 
transfer  coefficient: 

kg  =  {  [1.82  G/[(MW)m p\m,B\  }  [(dP  G)/p]  0  51  [p/(p  D)]  2/3  (187) 

The  correlative  relationship  given  by  eq  1 87  is  clearly  seen  to  contain  the 
dimensionless  groupings  of  the  Re  and  Sc.  The  characteristic  length  appearing  in  the  Re 
is  the  particle  diameter  of  the  adsorbent.  The  analog  of  eq  187  appearing  in  the  Dole  and 
Klotz  reference  (ref  3)  is  eq  7.  Equation  187  contains  the  logarithmic  mean  partial 
pressure  of  the  nontransferring  component  (pim^).  The  presence  of  this  term  and  how  it 
arises  in  the  analytical  relationships  associated  with  mass  transfer  is  addressed  in 
Appendix  D  of  this  review. 

Provided  that  no  contaminant  has  reached  the  effluent  of  the  bed.  Dole  and 
Klotz  maintain  that  the  rate  of  adsorption  of  i  onto  the  surface  must  be  equal  to  the  rate  at 
which  contaminant  is  introduced  into  the  adsorbent  system  at  the  inlet.  The  rate  of 
introduction  of  contaminant  i  at  the  inlet  on  a  mass  basis  can  be  expressed  as 

(m,)in  =  Q  C0  (MW),  (188) 

Equating  the  rate  at  which  mass  is  introduced  into  the  adsorbent  column 
(eq  188)  to  the  rate  at  which  mass  is  transferred  to  the  external  adsorbent  surface  (eq  186) 
and  solving  for  the  MTZ  thickness  gives 

Mz  =  (Q  Ca)/(kg  a  A  Apm)  (189) 

Substituting  the  mole-based  mass  transfer  coefficient,  as  given  by  eq  187, 
into  eq  189  gives  the  following  relationship: 

AT  =  [1/(1. 82  a)]  {  [(MW)m  Q  C0  plmM\/(G  A  A pm) }  Re0'51  Sc2/3  (190) 

Because  the  mass  velocity  of  the  airstream  goes  as  G  =  u  p  and  the 
volumetric  flow  rate  goes  as  Q  =  u  A,  eq  190  can  be  rewritten  as 

AT  =  [1/(1. 82  a)\  {  [(MW)„,  Cc]/p  }  OWApJ  Re0'51  Sc2/3  (191) 


The  quantity  {  [(MW)m  C„]/p  }  goes  as  the  vapor  phase  mole  fraction  of  i 
at  the  inlet  of  the  adsorbent  column  through  the  following  development; 


{  [(MW),,,  C0]/p  }  =  C,/[p/(MW)J  =  CJC  =  \(pdJp\ 

=  (ydo  =  y0 


(192) 
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Substituting  the  relationship  given  by  eq  192  into  eq  191  yields 


0.51  c„2/3 


Mz  =  [1/(1. 82  a)]  (y0  Pim.u/A/v)  Reu'31  Sc 


(193) 


Assuming  that  the  vapor  phase  can  be  characterized  as  an  ideal  mixture, 
then  ya  =  p„!p,  and  eq  193  becomes 


Mz  =  [1/(1. 82  a)\  [(poPim,B)/(p  A pm)]  Re0'51  Sc2/3 


(194) 


The  relationship  given  by  eq  194  is  very  close  to  eq  9  of  the  Dole  and 
Klotz  report  (ref  3).  The  relationships  differ  in  the  presence  of  the  total  pressure  term 
appearing  in  the  denominator  (in  the  Dole  and  Klotz  relationship,  this  pressure  factor  is 
absent). 


Dole  and  Klotz  indicate  that  the  following  relationship  is  valid  based  on 
engineering  experience: 


l(Po  Pim,«)/Ap,„]  =  In  (CJCh)  (195) 

The  relationship  given  by  eq  195  is  then  substituted  into  eq  194,  which 
gives  the  following  expression: 

Mz  =  [1/(1. 82  a )]  Re0'51  Sc2/3  In  (CJCh)  (196) 

It  is  observed  that  the  relationship  given  by  eq  195  lacks  dimensional 
consistency  in  that  the  collection  of  parameters  on  the  left  side  of  the  equality  would  have 
the  units  of  a  pressure,  while  the  logarithmic  ratio  on  the  right  side  of  the  equality  would 
be  unitless.  It  is  noted  that  the  relationship  given  by  eq  194  retains  dimensional 
consistency.  Some  rationalizing  considerations  for  the  form  of  the  relationship  given  by 
eq  195  are  postulated  in  the  next  section. 


7.3.5  Suggested  Rationalization  for  the  Validity  of  the  Relationship 

Given  by  Equation  194 

Dole  and  Klotz  give  the  following  relationship  (eq  1 1  of  ref  3)  for  the 
mean  partial  pressure  difference  of  the  transferring  gas  as  measured  from  the  main  gas 
stream  to  the  interface: 


A Pm  =  (Po-  pbWn(p0/pb)]  (197) 

Examination  of  the  relationship  given  by  eq  197  infers  that  it  can  be 
likened  to  a  logarithmic  mean  pressure  difference  taken  across  the  full  axial  length  of  the 
adsorbent  column.  It  can  also  be  interpreted  as  the  logarithmic  mean  pressure  difference 
within  the  fully  developed  MTZ  region  of  the  adsorbent  column. 
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Substituting  the  relationship  given  by  eq  197  into  the  pressure  ratio  term 
contained  in  the  relationship  provided  by  eq  194  gives 

[ip0  P\mJi)Kp  A Pm)]  =  [ Pol(Po-pb )]  (p\m.ll/p)  \  Mp,Jpb)\  (198) 

Assuming  that  the  partial  pressure  of  the  contaminant  in  the  inlet  flow  is 
much  greater  than  the  breakthrough  partial  pressure,  such  that  pa  »  pb,  then  the 
following  approximation  can  be  made: 


\PoKPo  -  Pb )]  ~  (pjpo)  =  1 


(199) 


Substituting  this  approximation  into  eq  198  yields 

l(Po  PunJi)/(p  A Pm)]  ~  (p\m,is/p)  [In (pjpb)\  (200) 

The  logarithmic  mean  partial  pressure  of  the  nontransferring  component  is 
given  as  eq  D-21  in  Appendix  D: 

P\m,B  =  (P 112  ~  Pb\  )/[MPB2/Pb\  )\ 

If  an  ideal  binary  vapor  phase  is  assumed  to  be  present  (where  A  is  the 
transferring  component  and  B  is  the  nontransferring  component),  the  following  partial 
pressure  relationships  can  be  written: 


and 


(Pb2-Pbi)  =  (Pai~Pa2) 


( PbtJpbx )  =  ip  ~  Pai)Kp  ~  Pa\) 

Substitution  of  these  expressions  into  the  relationship  given  by  eq  D-21  (in 
Appendix  D)  yields  the  following  equivalency  for  the  logarithmic  mean  partial  pressure 
of  the  nontransferring  component: 


P\m,B  =  (PAI  -PA2)/  {  ln[(p  ~Pa2)/(P  ~  PaO]  }  (201) 

Assuming  that  the  partial  pressure  in  the  bulk  fluid  phase  is  much  greater 
than  the  partial  pressure  of  A  near  the  surface  of  the  adsorbent  such  that  pA  1  »  Pai,  then 
eq  201  can  be  approximated  by 


P\m,B  ~  Pa\! {  ln[/V(p  —  Pai)]  }  (202) 

The  relationship  given  by  eq  202  can  be  reworked  through  the  following 

expressions: 
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(203) 


Pkn,B  ~  PAl/ {ln[p/(p  ~  PAl)]}  =/?Al/{ln{l/[l  ~(pAl/p)]}  } 

=  />Ai/{ln[l/(l  -yAi)]  }  =  -pAi/[ln(l  —  yAl)] 

The  logarithmic  term  of  eq  203  can  be  expressed  as  a  Taylor-series 
expansion  through  the  following  relationship: 

oo 

ln(l  - >’4i)  =  -£  (yAi)"/«  for  -1  <  yAl  <  1  (204) 

n  =  1 


The  first  several  terms  of  the  expansion  yield 


ln(l  -  yA1)  ~  ~yA\  -  (yM)2/2  -  (yA1)3/ 3  - .  (205) 

Truncating  the  series  after  the  first  term  on  the  basis  that  since  yA\  itself  is 
small,  then  the  subsequent  powers  of  yAi  would  be  negligible,  and  eq  205  can  be 
approximated  as 


yields 


ln(l -yAi)~-yAi  (206) 

Substitution  of  the  approximation  given  by  eq  206  into  eq  203  gives 

P\m,B  ~  PAllyM  (207) 

Because  pA  \  =  (p  yA\),  eq  207  further  simplifies  to 

P\m,B~P  (208) 

The  relationship  given  by  eq  208  is  then  substituted  into  eq  200,  which 


[  (Po  P\m,B)l(P  Apm)]  ~  [In (pjpb)]  =  [In (C0/Cb)]  (209) 

The  relationship  given  by  eq  209  is  consistent  with  the  final  expression  as 
proposed  by  Dole  and  Klotz  (ref  3)  by  their  eq  10.  It  is  suspected  that  a  pressure  term  was 
inadvertently  omitted  in  eq  10  of  the  Dole  and  Klotz  reference.  It  is  therefore  surmised 
that  the  relationship  given  by  eq  194  was  the  intended  form  of  the  MTZ  thickness  as 
developed  by  Dole  and  Klotz. 


7.3.6  Relationships  for  We 

Several  analytical  models  characterize  adsorbent  loading  as  a  function  of 
the  partial  pressure  of  the  contaminant  in  the  vapor  phase.  These  relationships  fall  under 
the  broad  heading  of  isotherms.  Isotherm  data  for  the  low- volatility  chemical  agent 
compounds  of  interest  typically  display  a  Type  1  Brunauer,  Emmett,  and  Teller  (BET) 
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isotherm,  which  is  generally  characterized  by  a  rapid  rise  in  the  adsorbent  loading  at  low 
partial  pressures  of  the  adsorbate.  This  rise  in  the  loading  eventually  reaches  a  plateau 
region  at  which  the  adsorbent  has  reached  a  saturation  capacity.  The  Type  1  BET 
isotherm  is  characterized  by  a  downward  concavity  on  a  plot  of  adsorbent  loading  versus 
partial  pressure  over  a  very  wide  range  of  partial  pressures.  The  generic  types  of  the  BET 
isotherms  are  illustrated  in  Figure  10  (ref  6). 


Figure  10.  Illustration  of  BET  isotherm  types  (from  ref  6). 


Other  loading  relationships  attributed  to  Langmiur  and  Dubinin  et  al.  have 
also  been  used.  The  removal  of  low-volatility  compounds  by  activated  carbon-based 
substrates  is  generally  well  characterized  by  the  Dubinin-Radushkevich  isotherm 
relationship. 


The  equilibrium  loading  capacity  of  the  adsorbent  for  the  adsorbate  under 
the  expected  conditions  of  use  (which  may  well  be  dynamic)  is  a  key  parameter  that  must 
be  quantified  for  use  in  breakthrough  estimations.  For  a  given  set  of  conditions,  the 
equilibrium  loading  capacity  places  some  bounds  on  the  amount  of  contaminant  that  can 
be  adsorbed.  Being  equilibrium-dependent,  these  loading  characteristics  are  taken  to  be 
rate-independent  for  fixed  conditions. 


7.3.7  Concluding  Comments:  Dole  and  Klotz  Analytical  Relationships 


The  relationship  given  by  eq  196  can  now  be  substituted  for  the  parameter 
h  in  eq  179  and  solved  for  the  breakthrough  time,  which  yields  the  following  expression: 


tb  ={(WeA  p b)/\Q  C0  (1  -  8)  (MW)/] } 
{L  -  { [1/(1. 82  a )]  Re0'51  Sc2/3  In  (CJCb) }  } 


(210) 


The  similarities  between  eq  210  and  those  expressions  that  are  currently 
used  to  estimate  the  breakthrough  times  for  strongly  adsorbed  components  are  noted.  The 
current  analytical  relationships  generally  incorporate  factors  to  account  for  the  non¬ 
sphericity  of  the  carbon  granules  and  different  exponents  for  the  dimensionless  Re  and  Sc 
numbers;  however,  the  fundamental  forms  are  equivalent. 
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Expressing  the  influent  challenge  concentration  on  a  mass  basis  such  that 
C0  =  C„  (MW),,  and  substitution  of  the  relationship  given  for  a  as  developed  by  eq  59  into 
eq  210,  yields  the  following  equivalent  form: 


tb={{WeApb)/[QC0(l-E)]} 

[L  -  {  { [dp/[6  (1.82)  (1  -  8)] }  Re0'51  Sc2/3  In  (C0/Cb) }  } 


(211) 


The  expression  given  by  eq  21 1  is  even  more  similar  to  the  current 
analytical  approaches  aimed  at  estimating  the  breakthrough  times  of  strongly  adsorbed 
compounds  onto  carbon  filters  for  military  applications. 

For  comparative  purposes,  eqs  212  and  213  were  extracted  from  page  17 
of  Peterson  and  Karwacki  (ref  11): 


tb  =  [(  WeA  pb)/(Q  C0)\  (L  -  AT) 


(212) 


where 


Mz  =  { (cp,  dp)/[6  (1  -  £)] }  [((p,  dp  u  p)/pf41  [p/(p  D)f67  In  (CJCb)  (213) 

A  comparison  of  the  forms  of  eqs  212  and  213  with  those  of  eqs  210  and 
211  indicates  a  strong  correspondence  and  a  historical  lineage  of  these  modem-day 
relationships  to  prior  analytical  efforts. 


8.  ANALYTICAL  EFFORTS  BY  WHEELER  ET  AL.  AND 

JONAS  ET  AL. 


8.1  B  ackground  Information 

Efforts  to  develop  analytical  modeling  relationships  aimed  at  adsorbent 
bed  optimization  were  undertaken  by  the  Department  of  the  Army,  Edgewood  Arsenal 
through  contractual  approaches.  These  efforts  were  conducted  during  the  1960s  under 
contract  DA  18-035-AMC-279(A)  (refs  12  and  15),  which  was  awarded  to  the 
AiResearch  Manufacturing  Company  (a  division  of  the  Garrett  Corporation,  Los  Angeles, 
CA)  by  the  U.S.  Army  Edgewood  Arsenal,  Chemical  Research  and  Development 
Laboratories  (now  known  as  the  U.S.  Army  ECBC). 

An  often-cited  relationship  that  apparently  resulted  from  these  efforts  is 
the  Wheeler  relationship,  including  subsequent  modifications  to  the  fundamental 
relationship.  Wheeler  (and  coauthors)  have  also  published  articles  in  technical  trade 
journals  (e.g.,  ref  14).  For  the  purposes  of  this  review,  Wheeler’s  efforts  toward 
breakthrough  relationships  formed  the  focus  of  the  analysis  described  in  the  next  section. 
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8.2 


Basic  Underlying  Analytical  Relationships 


As  with  previous  developments,  Wheeler  (and  others)  developed 
analytical  models  based  on  the  mass  balance  approach.  It  appears  that  integrated  forms  of 
fundamental  mass  balance  expressions  formed  the  basis  of  the  so-called  Wheeler  (and 
variant  forms)  relationships.  It  is  hypothesized  that  Dr.  Leonard  Jonas  (former  employee 
of  an  ECBC  predecessor  organization)  consulted  with  Wheeler  on  the  development  of  a 
breakthrough  relationship  that  would  have  practical  applicability  to  design  staff 
personnel.  The  source  of  this  hypothesis  rests  with  a  reference  footnote  appearing  within 
a  technical  article  published  in  the  journal  Carbon  by  L.  A.  Jonas  and  J.  A.  Rehrmann 
(ref  13)  that  reads,  “Wheeler,  A.,  private  communication.”  In  their  technical  article,  Jonas 
and  Rehrmann  give  the  following  equivalent  relationship  to  characterize  the  breakthrough 
concentration  for  a  nonselective  catalytic  bed  poisoning: 

ln(C-/C0)  =  [(kv  L)/u ]  { tb  /  [(p6  We  L)/(C0  nj]  }  “‘  (214) 

The  original  relationship  as  given  in  the  Jonas  and  Rehrmann  reference 
used  a  different  nomenclature  convention  than  that  appearing  in  eq  214.  The  relationship 
as  given  by  eq  214  reflects  the  nomenclature  convention  that  has  been  used  consistently 
throughout  this  current  review  effort.  This  relationship  corresponds  to  eq  1  of  the  Jonas 
and  Rehrmann  reference  (ref  13). 

Jonas  and  Rehrmann  cite  Wheeler  and  Robell  (ref  14)  as  the  source  of  the 
relationship  given  by  eq  214. 

The  following  relationship  appears  as  eq  2  in  the  Jonas  and  Rehrmann 
reference  (appropriate  variable  correspondences  were  made  to  be  consistent  with  current 
nomenclature  convention): 

tb  =  [(p/,  We)/(C0  u)\  [L  -  (u/kv)  ln(C0/Cz)J  (215) 

The  relationship  given  by  eq  214  can  be  rearranged  to  yield  eq  215; 
however,  the  methodology  is  not  explicitly  contained  within  the  Jonas  and  Rehrmann 
reference  (ref  13).  The  development  of  the  relationship  given  by  eq  215  is  attributed  to  a 
private  communication  with  Wheeler.  It  is  noted  that  eq  215  is  of  a  form  that  is  consistent 
with  those  relationships  attributed  to  Wheeler. 

8.3  Macroscopic  Mass  Balance  Considerations 

Early  efforts  undertaken  by  the  AiResearch  Manufacturing  Company 
under  contract  to  the  U.S.  Army  Edgewood  Arsenal  during  the  1960s  provide  the  basis 
for  how  the  relationship  given  by  eq  215  was  developed.  In  the  development  that  follows, 
the  underlying  principles  upon  which  the  Wheeler  expressions  are  based  are  described  in 
part.  This  development  is  based  on  Wheeler  et  al.  (ref  15). 

The  total  mass  of  contaminant  i  that  enters  the  reactor  vessel  over  a  time 
interval  t  is  given  by 
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(Mi) in  =  (MW)f  C„  Q  t 


(216) 


given  by 


The  total  mass  of  i  that  exists  in  the  adsorbed  state  within  the  reactor  is 


L  L 

(M/U  =  /  Wt(z,  t)  pb  A  Az  =  pb  A  /  Wfe  0  dz  (217) 

0  0 

The  cross-sectional  area  and  the  bulk  density  have  been  taken  as  constants 
in  eq  217.  The  functional  dependence  of  the  specific  mass  loading  on  both  axial  position 
and  time  has  been  explicitly  retained  in  eq  217. 

The  total  mass  of  i  that  exits  the  reactor  vessel  at  axial  position  z  =  L  is 

given  by 


t 

(AT)en  =  /  C,(z  =  L,t)Q  (MW),-  d t  (218) 

0 

The  variable  C,{z  =  L,  t)  is  interpreted  as  the  vapor  phase  concentration  of 
contaminant  i  evaluated  at  axial  position  z  =  L.  In  the  subsequent  analytical  relationships, 
the  following  variable  correspondence  is  used: 


Q(z  =  L,t)  =  Ci(L)  (219) 

When  the  variable  correspondence  given  by  eq  219  is  inserted  into  eq  218 
and  the  volumetric  flow  rate  is  assumed  to  be  constant,  eq  218  can  be  reexpressed  by 

t 

(M/)e ff  =  (MW)/  Q  /  Ci(L)  At  (220) 

0 

The  macroscopic  mass  balance  over  the  entire  reactor  length  (L)  is  given 
by 


(Mp  in  —  (Mj)  ads  +  {Mj)e  ff 


Substitution  of  the  relationships  given  by  eqs  216,  217,  and  220  into  the 
macroscopic  mass  balance  relationship  yields 

L  t 

(MW)/  C0  Q  t  =  pbA  /  Wiiz,  t)  Az  +  (MW),-  Q  /  C,(L)  At  (221) 

0  0 
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Wheeler  et  al.  (ref  15)  indicate  that  at  very  high  inlet  concentrations,  an 
additional  term  would  have  to  be  incorporated  into  eq  221  to  account  for  the  mass  of 
agent  that  exists  within  the  vapor  phase  along  the  length  of  the  reactor.  This  additional 
mass  term  would  have  the  form 


L 

(Mi), ap  =  (MW);  A  /  Ci(z)  dz  (222) 

0 

Wheeler  et  al.  infer  that  at  reasonably  expected  inlet  concentrations,  the 
bulk  of  the  agent  mass  will  reside  on  the  adsorbent  along  the  reactor  length.  Given  that 
this  approximation  is  valid,  the  mass  contribution  described  by  eq  222  would  be 
negligible  relative  to  the  adsorbed  phase  mass  as  given  by  eq  217. 

The  saturation  loading  capacity  of  the  carbon  contained  within  the  system 

is  given  by 


LApbWe  (223) 

Adding  and  subtracting  the  saturation  loading  capacity  of  the  carbon  (the 
relationship  given  by  eq  223)  to  the  right  side  of  the  equality  of  the  relationship  given  by 
eq  221  gives 


L 

(MW);  C0  Q  t  =  pb A  / Wi(z,  t)dz  -  LA  pb  We 

0 


t 

+  (MW);  Q  / C,(L)  d t  + LA  pb  We 
0 


(224) 


The  relationship  given  by  eq  223  can  be  equivalently  expressed  through 
the  following  correspondence: 


L 

LApb  We=pbAjwedz  (225) 

0 

Substitution  of  the  relationship  given  by  eq  225  into  eq  224  yields 
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L  L 

(MW);  Co  Q  t  =  ph A  / Wi(z,  t)  dz  ~  pb A  J Wedz  +  LA  pb  We 
0  0 


L 

=  phAf[Wi(z,t) 
0 


t 

+  (MW),  Q  /  C,<L)  d t 
0 


We]dz  +  LAphWe 


t 

+  (MW);  Q  /  CiL)  dt 
0 


(226) 


The  following  mathematical  correspondences  are  noted  (where  the  explicit 
functional  dependence  of  the  localized  specific  loading  has  been  dropped  for 
convenience): 


Wi  -We  =  -(We-  Wj)  = -We  [  1  -  (Wi/We)\  (227) 

Substitution  of  the  relationship  given  by  eq  227  into  eq  226  and 
recognizing  that  We  is  a  constant,  yields  the  following  expression: 

L  t 

(MW);  C0  Q  t  =  pbA  We  L~  pb  A  We  /[I  -  (Wi/We)]  d z  +  (MW),  Q  / C,(L)  dt  (228) 

0  0 

Dividing  eq  228  by  the  quantity  [ Q  C0  (MW),]  yields 
t=[(pbA  We)/(Q  C0  (MW),)]  L 

L 

-  [(p bA  We)/(Q  C0  (MW);)]  /[I  -  (WJWe)\  dz 

0  (229) 

t 

+  (i /c0)  /  cm  dt 
0 


Recognizing  that  (Hu)  =  A/Q  and  that  Ca  is  a  constant,  eq  229  can  be 
expressed  in  the  following  modified  form; 
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L 

t  =  [(p*  We)/(u  Co  (MW),)]  {  L  -  J[1  -  (Wi/We)]  dz 

0 


t 

+  [(«  Co  (MW)f)/(pi  We)]  / [Q(L)/Co]  dr} 

0 


(230) 


The  relationship  given  by  eq  230  is  eq  12-23  in  the  Wheeler  et  al. 
reference  (ref  15).  The  development  of  this  relationship  was  based  on  the  assumption  that 
the  mass  of  agent  present  in  the  vapor  phase  along  the  reactor  length  was  negligible 
relative  to  the  mass  of  agent  in  the  adsorbed  state.  It  is  observed  that  many  other 
commonly  accepted  mass  balance  relationships  account  for  the  presence  of  the  adsorbate 
in  the  vapor  phase  within  the  interstitial  volume  of  the  adsorbent  packing. 

The  Wheeler  et  al.  reference  (ref  15)  gives  the  following  interpretations  to 
the  parameter  groupings  contained  in  eq  230: 

[(p*  We)/(u  Co  (MW),)]  L 

represents  the  time  required  to  saturate  the  adsorbent  column  with  infinitely  rapid 
adsorption.  The  relationship 


L 

Up/,  We)/(u  Co  (MW),)]  /[I  -  (WJWe)]  dz 

0 

represents  a  correction  factor  to  account  for  that  segment  of  the  bed  which  cannot  be  fully 
used;  and 


t 

/  [Ci(L)/C0]dt 

0 

represents  the  amount  of  agent  exiting  the  reactor. 

For  an  adsorbent  bed  having  prior  contaminant  loading,  the  localized 
contaminant  loading  is  reduced  by  the  specific  loading  of  the  contaminant  (W0). 
Therefore, 


( IF,') contain  mated  bed  —  (  W/)  clean  bed  W„  (231) 

Substituting  the  expression  given  by  eq  231  and  the  following  total 
differential  relationships 

dz  =  {dCi!dz)~l  dCj  and  d t  =  (dCJdl)  1  dC, 
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into  eq  230  gives  the  following  expression  for  an  adsorbent  bed  having  prior  contaminant 
loading: 


Cb  L 

tb  =  Kpb  We)/(u  Cn  (MW),)]  {L  -  /[I  -  C Wi/We)\  0 dQ/dzT 1  d Q  -  f  (WJWe)  d z 

Co  0 


Cb 

+  [(«  c„  (MW),)/(p b  we)]  /  [c,<l)/c0]  (dCi/dty1  dc,} 


(232) 


Through  a  series  of  mathematical  manipulations,  rationalization 
arguments,  the  targeting  of  resultant  relationships  to  a  form  consistent  with  the  basic 
Mecklenburg  mass  balance,  and  introduction  of  an  adsorption  rate  functionality,  Wheeler 
et  al.  (ref  15)  give  the  following  breakthrough  relationship  for  low-break  concentrations 
on  an  initially  clean  adsorbent  bed: 


Co 


tb  =  [(Pi  We)/(U  C, 


(MW)/)]  {  L  -  (/ u/kv )  /  { [1  -  (Q/CoWiCi  V|/(C,)] }  dQ 


Cb 


Cb 

+  (u/kv)  f  {  1/[C0  V|/(C,)] }  dCi } 

0 


(233) 


Further  assumptions  with  this  expression  are  the  direct  proportionality 
between  the  adsorbent  mass  loading  ratio  to  the  vapor-phase  concentration  ratio  and  the 
incorporation  of  a  steady-state  function  to  characterize  the  adsorption  kinetics.  This 
steady-state  functionality,  depicted  by  the  symbol  \|/(C,),  is  dependent  on  the 
concentration  reduction  ratio,  defined  as  F  =  C/Ca  such  that 


\|/  =  fcn(C,)  =  \| I(Q)  =  \| ){F)  (234) 

The  development  of  the  relationship  given  by  eq  233  (for  an  initially  clean 
bed)  is  rather  formidable  and  was  not  well  understood  by  this  author. 

Substituting  the  defined  relationship  for  F  into  eq  233  yields 


1 

h  =  [(p*  We)/(u  C0  (MW),)]  {  L  -  (i u/kv )  /  { (1  -  F)/[F  \|/(F>] }  dF 

Fh 


Fb 

+  (u/kv)  j [  1/V|/(F)|  dF} 

0 


(235) 
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Upon  comparison  of  eq  235  to  the  basic  Mecklenburg  relationship 
(eq  146),  the  depth  of  the  MTZ  thickness  is  given  by 


1  Fb 

Mz  =  ( ulkv )  {/(  (1  -  F)/[F  V| 1(F)] }  d F  -  J[1/V|/(F>]  d F} 
Fb  0 


(236) 


The  relationship  given  by  eq  236  is  given  by  eq  12-48  of  the  Wheeler  et  al. 
reference  (ref  15). 

For  the  special  case  in  which  \| f(F)  =  (1  -  F),  eq  236  yields 

Mz  =  (; ulkv )  In [(C0/Cb)  -  1]  (237) 

The  relationship  given  by  eq  237  is  termed  as  the  “modified 
Mecklenburg”  form  of  the  MTZ  thickness.  Wheeler  et  al.  indicated  that  the  modified 
Mecklenburg  expression  is  valid  only  if  the  velocity  of  propagation  of  the  adsorption 
wave  is  constant  from  t  =  0  onward.  This  is  pointed  out  to  be  a  poor  assumption,  because 
the  velocity  of  propagation  of  the  adsorption  wavefront  is  somewhat  less  than  the  steady- 
state  propagation  associated  with  a  uniformly  developed  wavefront.  The  uniformity  of  the 
wavefront  develops  with  time  and  distance  into  the  column.  Prior  to  the  establishment  of 
a  uniformly  characterized  adsorption  wavefront,  the  contours  of  the  adsorption  wave  are 
dynamic. 


Upon  consideration  of  several  assumed  analytical  forms  to  characterize 
the  adsorption  kinetics  (see  relationship  given  by  eq  234),  Wheeler  et  al.  were  able  to 
ascertain  that  the  ratio  (u/kv)  naturally  arose.  This  occurrence  suggested  a  form  of  the 
Mecklenburg-based  breakthrough  relationship  consistent  with  the  following  relationship: 

tb  =  [(pb  We)/(u  C0  (MW),)]  {  L  -  ( u/kv )  \n[K  {CJCb)\  }  (238) 

where  K  is  a  multiplicative  factor  that  is  dependent  on  the  functional  form  assumed  for 
the  concentration  reduction  ratio  (F).  The  relationship  given  by  eq  238  appears  as 
eq  12-54  of  the  Wheeler  et  al.  reference  (ref  15).  The  ratio  ( C„/Cb )  has  been  assumed  to 
be  much  larger  than  1  in  the  relationship  given  by  eq  238. 

Figure  1 1  contains  a  table  appearing  in  the  Wheeler  et  al.  reference  that 
denotes  various  values  for  the  K  parameter  corresponding  to  assumed  forms  of  the 
adsorption  rate  function.  The  symbol  La  is  the  depth  of  the  MTZ,  which  has  been 
indicated  by  the  symbol  M-  throughout  this  review  effort. 
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Figure  11.  Functional  forms  of  \| f(F)  and  resultant  MTZ  thicknesses  (from  ref  15). 


The  variation  in  the  value  of  the  parameter  K  with  respect  to  the  assumed 
functionality  of  \| f(F)  is  readily  apparent  in  Figure  11.  When  used  in  conjunction  with  the 
Wheeler  et  al.  relationship  given  by  eq  238,  this  can  cause  quite  a  bit  of  variance  with  the 
estimated  breakthrough  time.  The  evaluation  of  the  rate  constant  kv  would  be  determined 
experimentally. 

8.4  Jonas  and  Rehrmann  Breakthrough  Relationship 

In  their  investigative  efforts  to  characterize  the  breakthrough  of 
organophosphorous  compounds  through  activated  carbon  beds,  Jonas  and  Rehrmann 
(ref  13)  used  the  Wheeler-based  relationship  as  given  by  eq  238  with  K=  1.  This  yields 
the  following  relationship: 

h  =  [(p b  We)/(u  C0  (MW),)]  {  L  -  ( u/kv )  In (C0/Cb) }  (239) 

From  the  tabular  entries  contained  in  Figure  1 1,  a  parameter  value  of 
K  =  1  corresponds  to  the  following  functional  form  for  \| f(F): 

\| /(F)  =  1  -F  (240) 
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Wheeler  et  al.  (ref  15)  indicated  that  the  functional  breakthrough 
relationship  given  by  eq  239  contains  some  rather  significant  assumptions,  the  major  one 
being  the  assumption  of  a  steady-state  propagation  velocity  of  the  adsorption  wave  at 
time  t  =  0.  Wheeler  et  al.  indicated  that  an  initial  time  interval  and  adsorbent  bed  depth 
are  required  to  establish  a  consistent  wave  contour  for  an  initially  clean  adsorbent 
column. 

8.5  Closing  Commentary  on  the  Wheeler-Based  Relationships 

As  indicated  earlier  in  this  section,  significant  portions  of  the  analytical 
manipulations  and  rationalizations  pursued  by  Wheeler  et  al.  were  not  well  understood. 
Although  the  development  of  the  basic  mass  balances  (given  by  eqs  230  and  232  are 
relatively  straightforward,  much  of  the  later  analytical  development  is  of  a  complex 
nature. 


An  assumption  relating  to  the  basic  Wheeler  and  Jonas  breakthrough 
relationships  is  that  the  amount  of  contaminant  that  exists  in  the  vapor  phase  (i.e.,  the 
interstitial  volume  of  the  adsorbent  packing)  is  negligible  relative  to  the  amount  that 
exists  in  the  adsorbed  state.  Wheeler  et  al.  qualify  this  approximation  in  their 
development  (ref  15).  It  is  noted  that  other  mass  balance  approaches  also  include  the 
mass  of  contaminant  that  is  contained  in  the  interstitial  volume. 

A  potential  weakness  of  the  Wheeler  et  al.  and  Jonas  et  al.  breakthrough 
relationships  is  their  reliance  on  having  access  to  actual  data  to  estimate  the  rate  constant 
kv.  This  data  is  often  not  readily  available  to  the  typical  practicing  engineer  in  making 
filter  performance  estimations. 

The  Klotz  et  al.  mass  balance  relationships  appear  to  be  more  user-friendly 
in  that  the  expressions  require  the  evaluation  of  specific  parameters  through  readily 
accepted  correlative  techniques.  The  currently  used  Klotz-based  relationships  assume  that 
bulk  diffusion  is  the  rate-limiting  step. 


9.  CONSIDERATIONS  WITH  CYANIDE-BASED  CONTAMINANTS 

Sections  2  through  8  of  this  review  focused  on  the  adsorption  of 
compounds  that  are  assumed  to  be  very  favorably  removed  from  a  contaminated 
airstream  through  an  idealized  physical  adsorption  pathway.  The  notable  exceptions  were 
chlorine  (ref  7)  and  phosgene  (ref  3),  in  which  hydrolysis  and  chemisorption  mechanistic 
pathways  contribute  considerably  to  the  removal  of  these  materials  from  a  contaminated 
air  flow.  With  compounds  having  a  room  temperature  vapor  pressure  of  less  than  10  torr, 
the  primary  removal  mechanism  is  often  attributed  to  physical  adsorption.  The  rate  of 
removal  is  often  assumed  to  be  diffusion-limited  in  the  bulk  phase. 

The  cyanide  materials  having  elevated  room- temperature  vapor  pressures 
require  the  presence  of  specific  impregnant  materials  within  the  pore  structure  of  the 
activated  carbon  provide  capacity  against  this  class  of  compounds.  A  requisite 
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characteristic  for  activated  carbons  that  can  be  impregnated  to  increase  their  capacity 
against  the  cyanide-containing  compounds  is  that  the  basic  substrate  possess  a  sufficient 
macropore  volume.  This  macroporous  structure  is  needed  to  provide  ready  access  to  the 
vapor  contaminants  into  the  porous  structure  of  the  carbon  and  to  the  intraporous  regions 
that  contain  the  impregnants. 

As  the  analytical  expressions  for  the  strongly  adsorbed  materials  are  based 
on  diffusive-based  rate  relationships,  these  expressions  are  not  appropriate  to  use  with  the 
cyanide -based  compounds. 

Figure  12  (ref  16)  depicts  a  series  of  life-thickness  curves  for  Cyanogen 
Chloride  (CK)  on  a  current  impregnated  carbon  substrate.  This  plot  is  based  on  data 
obtained  from  the  developmental  effort  for  ASZM-TEDA  impregnated  activated  carbon 
with  a  particle  size  range  within  the  12  x  30  mesh  envelope. 
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Figure  12.  CK  life-thickness  curves  for  12  x  30  ASZM-TEDA  carbon  (from  ref  16). 


Figure  12  is  analogous  to  Figure  9  of  this  review  effort.  The  fitted 
trendlines  exhibit  slight  curvature.  Linear  trendlines  could  also  be  fitted  (as  in  the  case  of 
Figure  9)  that  would  clearly  cross  the  abscissa  at  definitive  values  corresponding  to  a 
critical  bed  depth. 

With  an  elevated  room-temperature  vapor  pressure,  CK  exhibits  greatly 
reduced  physical  adsorption  capacity  as  compared  to  low-volatility  organics  under  near- 
atmospheric  conditions.  The  addition  of  impregnant  materials  (to  the  activated  carbon 


69 


substrate)  is  needed  to  provide  capacity  to  the  carbon  against  high- volatility  and  reactive 
chemical  contaminants. 

The  fundamental  mass  balance  relationships  developed  in  the  various 
sections  of  this  report  can  be  used  to  characterize  the  vapor  phase  concentrations  of  CK 
throughout  the  adsorbent  bed  as  a  function  of  position  and  time  provided  that  rate 
parameters  associated  with  the  destruction  of  CK  (i.e.,  the  r|  and  R,  terms)  are  known  or 
can  be  determined  through  experimentation.  As  has  been  demonstrated,  the  complexity  of 
the  resulting  analytical  relationships  is  formidable.  The  energetic  aspects  associated  with 
the  adsorption  and  heats  of  reaction  should  also  be  considered  for  reactionary  removal 
processes  (as  well  as  for  removal  mechanisms  that  are  characterized  by  ideal  physical 
adsorption). 


The  rate-limiting  step  associated  with  the  removal  of  the  cyanide 
compounds  is  believed  to  be  due  to  a  combination  of  initial,  limited  physical  adsorption 
followed  by  a  sequence  of  chemical  reactionary  mechanisms.  This  is  compared  to  a  mass 
transfer  rate-limiting  mechanistic  sequence.  It  is  generally  accepted  that  the  basic 
activated  carbon  substrate  must  process  a  broad  distribution  of  pore  sizes  to  have 
sufficient  capacity  against  CK  under  a  range  of  exposure  conditions.  The  large  pores  are 
needed  to  provide  rapid  access  of  the  contaminated  airstream  to  the  meso-  and 
micropores  where  the  impregnants  typically  reside.  It  is  within  these  meso-  and 
micropores  that  the  majority  of  the  porous  volume  resides. 

A  simplistic  approach  at  estimating  the  CK  performance  of  a  filter  is  to 
relate  the  actual  flow  characteristics  of  the  filter  to  the  closest  corresponding  velocity 
trend  line  appearing  on  a  life-thickness  plot  such  as  provided  in  Figure  12.  The 
normalizing  parameter  is  taken  to  be  the  residence  time.  Once  a  pseudo-bed  depth  is 
determined  (through  the  residence-time  normalization),  a  breakthrough  time  can  be 
estimated  through  a  correlative  fit  to  the  data  set. 

Filter  performance  estimations  against  hydrogen  cyanide  are  generally 
conducted  in  an  analogous  manner.  The  key  requirement  would  be  to  have  actual  life¬ 
thickness  data  at  the  flow  conditions  that  mimic  the  actual  system. 
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NOMENCLATURE  AND  ABBREVIATIONS 


a 

8 

8 

n 

e 

0  (subscript) 
0 
P 
Pi 

V 


Cp,' 

¥ 


dy  dp 

d 

AB  (subscript) 
ads  (subscript) 
At 

b  (subscript) 


C 

CG 

Q 

Co 


cyl  (subscript) 
(l  *  dp 
D 


porosity  or  fraction  of  voids  per  unit  gross  volume  of  bed 

fraction  of  the  total  adsorbent  system  volume  that  can  be  utilized 
for  adsorptive  capacity 

volume -based  void  fraction 

moles  of  toxic  gas  on  or  in  the  granules  contained  in  a  unit  volume 
of  bed 

circumferential  direction 
circumferentially  directed  flow 
fluid  phase  viscosity 

density  of  airstream  or  density  of  gas  film 
density  characteristic  associated  with  i 

analogous  to  a  superficial  flow  velocity  related  to  the  rate  of 
supply  (as  used  in  Section  5)  of  chlorine  per  length  of  charcoal 
saturated  per  time 

sphericity  factor 

symbol  used  for  collection  of  parameters  =  [(<9C,/o0)  (1/r)];  also, 
function  to  characterize  the  kinetics  of  adsorption.  As  used  in 
further  developments,  the  assumed  form  of  this  function  is  taken  to 
be  a  function  the  vapor  phase  concentration  reduction  ratio  (as 
used  in  Section  8) 

ratio  of  particle  surface  area  to  volume  of  adsorbent  bed 
residual  volume-based  chlorine  capacity  of  the  charcoal  (as  used  in 
Section  5) 

component  A  with  respect  to  B 
of  or  relating  to  the  adsorbent 

surface  area  associated  with  characteristic  i  descriptor 

of  or  relating  to  a  bulk  system  characteristic  such  as  a  bulk  density 
or  when  used  with  a  molar  concentration,  to  designate  a 
breakthrough  condition 
total  molar  concentration 
phosgene 

vapor  phase  molar  concentration  of  component  i 

vapor  phase  concentration  of  challenge  flow  expressed  on  a  mass 

basis 

of  or  pertaining  to  a  cylinder 

particle  diameter  assumed  to  be  an  ideal  sphere 

binary  diffusivity 
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B, 

% 

Bab 

Bl 

e  (subscript) 

F 

F 

G 

h 

i  (subscript) 
in  (subscript) 

h 

j  (subscript) 

Ja 

K 

k 

kv 

kg 

L 

MTZ 

Mt 

Mz 

(MW), 

(MW)m 

Vdi 

Nij 

net  (subscript) 
na 

o  (subscript) 
out  (subscript) 


generalized  dispersion  coefficient  of  component  i  in  a  binary 
system 

dispersion  coefficient  of  component  i  in  the  vapor  phase  associated 
with  direction  j 

binary  diffusion  coefficient  of  component  A  in  component  B 
uniform,  nondirectional  dispersion  coefficient 
effluent  value  or  evaluated  at  equilibrium  condition 

generic  mass  transfer  coefficient,  {  moles  of  i  transferred/[area 
time] } 

concentration  reduction  ratio,  CJCn  (as  used  in  Section  8) 
mass  velocity  of  airstream,  product  of  a  density  and  velocity 
height  of  cylinder  or  critical  bed  depth 
designator  for  component  i  or  generic  designator 
associated  with  a  flow  or  flux  entering  into  volume  element 

critical  bed  thickness  contribution  due  to  specific  adsorptive 
mechanistic  pathway  x 

designator  associated  with  a  specific  direction  in  space 
molar  flux  of  component  A 

multiplicative  factor  that  is  a  function  of  the  concentration  ratio  F 
of  the  Wheeler  et  al.  modeling  development  (as  used  in  Section  8) 

rate  constant 

first-order  rate  constant  (as  used  in  Section  8) 
gas-phase  mass  transfer  coefficient, 

{  moles  of  i  transferred/[area  time  (driving  force)]  } 

bed  depth  or  length  of  an  adsorbent  column  also  used  as  a 
volumetric  flow  rate 

mass  transfer  zone 
mass  of  component  i 

mass  transfer  zone  or  critical  bed  depth/thickness 

molecular  weight  of  component  i 

mean  molecular  weight  of  binary  vapor-phase  system 

rate  of  mass  flow  associated  with  component  i 
molar  flux  of  component  i  associated  with  direction  j 
net  characteristic 

saturation  capacity  of  carbon  for  contaminant  i,  expressed  in  moles 

of  i  per  volume  of  adsorbent 

influent  concentration  or  initial  condition 

associated  with  a  flow  or  flux  exiting  from  a  volume  element 
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p 

Po 
PAi 
Plm,B 
A Pm 


Q 

Qi 

r  (subscript) 

Re 

Ri 

r 

S 

Sc 

t 

t  (subscript) 
tot  (subscript) 
total  (subscript) 

Uj 

V 

V 

Vi 

void  (subscript) 


Wi 

w 

X 

X0 


x,  y,  z  (subscript) 


z 

Zo 


total  pressure 

partial  pressure  of  gas  entering  adsorbent  system 

partial  pressure  of  component  A  at  position  i 

log  mean  partial  pressure  of  nontransferring  species  in  gas  film 

mean  partial  pressure  difference  of  transferring  compound 
measured  from  the  bulk  vapor  phase  to  a  concentration  measured 
at  the  interface 
volumetric  flow  rate 

mass  loading  of  adsorbate  i  per  volume  of  adsorbent 
radially  directed  flow 
Reynolds  number 

consumption  rate  of  component  i  on  a  volumetric  basis 
radius  or  radial  direction 
arc  length  (r  9) 

Schmidt  number 
time 

critical  bed  thickness  contribution  due  to  bulk  diffusion 
of  or  relating  to  a  total  system  characteristic 
of  or  relating  to  a  total  parameter  characteristic 
averaged  superficial  velocity  associated  with  direction  j 
volume 

volume  associated  with  characteristic  i  descriptor 

fluid  flow  velocity  associated  with  direction  i 

of  or  relating  to  the  void  space  associated  with  a  granular  packing 

of  solids 

mass  loading  of  adsorbate  i  per  mass  of  adsorbent 
rate  of  mass  transfer 

characteristic  edge  dimension  of  FCC  unit  cell 

bed  depth  thickness  over  which  active  mass  transfer  is  occurring, 

critical  bed  depth,  mass  transfer  zone  thickness 

three  mutually  perpendicular  axis  in  an  orthogonal  coordinate 
system  or  associated  with  a  directional  flow 

vapor  phase  mole  fraction  of  component  A 
directional  axis 

bed  depth  thickness  over  which  active  mass  transfer  is  occurring, 
critical  bed  depth,  mass  transfer  zone  thickness 
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Blank 
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APPENDIX  A 

FLUX  GRADIENT  IN  THE  RADIAL  DIRECTION 


The  intent  of  this  appendix  is  to  detail  the  analytical  development  of  the 
radial  flux  gradient  given  by  eq  38.  Per  the  format  provided  by  the  eq  1 1  relationship,  the 
radial  flux  gradient  would  be  expressed  in  terms  of  an  arc  length  ( 5 ).  This  appendix 
provides  the  mathematical  details  in  relating  the  flux  gradient  expressed  in  terms  of  an 
arc  length  to  an  angle  (9). 

Dimensional  consistency  requires  that  the  gradient  of  the  molar  flux 
associated  with  the  radial  direction  have  the  units  of  moles  of  //(volume  time).  Using  the 
form  of  eq  11,  the  flux  gradient  in  the  radial  direction  would  formally  take  on  the 
following  relationship: 


(dNn/dS)  =  [d(Ui)  Q)/dS ]  -  Bl  (52C,/552)  (A- 1) 

where  D y  has  been  assumed  to  be  constant-valued  and  equal  to  BL. 

The  arc  length  is  related  to  its  angular  and  dimensional  components 
through  the  following  expression: 


5  =  r  0  (A-2) 

The  total  differential  of  eq  A-2  is  given  by 

dS  =  (55/50)  d9  +  (55/5r)  dr  (A-3) 

As  an  expression  for  the  change  in  5  with  respect  to  angular  position  (9)  is 
sought  at  a  constant  radial  position,  eq  A-3  simplifies  to 

(d5/d0)  =  (55/50)  (A-4) 

From  the  relationships  given  by  eqs  A-2  and  A-4,  the  following 
derivatives  can  be  written  for  a  constant  radius: 

(55/50)  =  r  or  (50/55)  =  1/r  (A-5) 

With  the  relationship  given  by  eq  A-5,  the  following  correspondence 

holds: 


[5(m0  CO/55]  =  [5(h0  Q/50]  (50/55)  =  (1/r)  [5(«0  C0/50]  (A-6) 

Similarly,  the  following  correspondence  holds: 

(5AW55)  =  (5AW50)  (50/55)  =  ( 1/r)  (5/VlO/50)  (A-7) 
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2  2 

With  the  relation  given  by  eq  A-5,  the  derivative  (d  Cj/dS  )  can  then  be 
related  to  the  angular  position  9  through  the  following  development; 


0 d2Q/dS 2)  =  [0(0C,/0S)/0S]  =  0[(0C,/0 0)  (00/0S)]/0S 

=  o[(oCi/o0)  (l/r)]/as 


(A- 8) 


Letting  \| /  =  [(0C,/S0)  (1/r)],  eq  A-8  can  be  further  developed  as 


0[(0C,/00)  (l/r)]/0S  =  (5\[/  IdS)  =  (Sv|//c0)  (50/05)  =  (0Y|//00)  (1/r) 
=  (1/r)  (0y/00)  =  (1/r)  {  0[(0C;/00)  (l/r)]/00  } 


(A- 9) 


Since  the  radius  is  a  constant,  it  can  be  removed  from  the  derivative  of 
eq  A-9,  which  yields  the  following  series  of  relationships: 


(< d2Q/dS 2)  =  0[(0C,/00)  (l/r)]/50  =  (1/r2)  {  0[(0C//00)]/00  } 

=  (1/r2)  (02C,/002) 


(A- 10) 


The  molar  flux  gradient  in  the  radial  direction  can  therefore  be  related  to  9 
through  the  relationships  given  by  eqs  A-6,  A-7,  and  A- 10.  This  correspondence  is  as 
follows: 


(0AW0S)  =  (1/r)  (0AW00) 

=  (1/r)  [0(m9  C,)/S0]  -  Dl  (1/r2)  (02C//002) 


(A- 11) 


APPENDIX  A 


78 


APPENDIX  B 

CONTINUITY  RELATIONSHIP  FOR 
CYLINDRICAL  COORDINATE  SYSTEM 


This  appendix  provides  the  rationale  for  the  relationship  given  by  eq  43. 
The  analytical  development  contained  in  this  appendix  is  based  on  the  cylindrical 
configuration  depicted  in  Figure  2.  It  is  further  assumed  that  the  differential  volume 
element  is  spatially  fixed  within  its  coordinate  reference  frame. 

With  respect  to  the  z-direction,  the  rate  at  which  mass  enters  the 
differential  volume  element  at  position  z  is  given  by 

p  uz(z )  r  A9  A r  (B-la) 

The  rate  at  which  mass  exits  the  differential  volume  element  in  the 
--direction  at  position  (z  +  A z)  is  given  as 

p  uz(z  +  Az)  r  A9  A r  (B-lb) 

With  respect  to  the  9-direction,  the  rate  at  which  mass  enters  the 
differential  volume  element  at  position  9  is  given  by 

p  «e(9)  A r  Az  (B-2a) 

The  rate  at  which  mass  exits  the  differential  volume  element  in  the 
9-direction  at  position  (9  +  A9)  is  given  as 

p  we(9  +  A9)  A r  Az  (B-2b) 

The  rate  at  which  mass  enters  the  differential  volume  element  in  the 
r-direction  at  position  r  is  given  by 


p  u,(r)  r(r)  A9  Az  (B-3a) 

The  rate  at  which  mass  exits  the  differential  volume  element  in  the 
r-direction  at  radial  position  (r  +  Ar)  is  given  as 


p  u,(r  +  Ar)  r(r  +  Ar)  A9  Az  (B-3b) 

The  total  mass  of  material  contained  within  the  differential  volume 
element  is  given  by 


M  =  p  AV  =  p  r  A9  Ar  Az 


(B-4) 
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by 


The  rate  of  change  of  mass  within  the  differential  volume  element  is  given 


(dM/dt)  =  A9  A r  A z  [o(p  r)/dt\  (B-5a) 

If  the  mass  density  within  the  differential  volume  element  is  taken  to  be  a 
constant  (mass  conservation),  eq  B-5a  becomes 

(dM/dt)  =  p  A0  A r  A z  ( dr/dt )  (B-5b) 

The  mass  balance  relationship  written  over  the  differential  volume 
element  is  expressed  as 


(dM/dt)  =  £(m()in  -  2(771;) out  (B'6) 

Substitution  of  the  relationships  given  by  eqs  B-l,  B-2,  B-3,  and  B-5  into 
the  relationship  given  by  eq  B-6  gives  the  following  expression: 

p  A9  Ar  Az  (dr/dt)  =  -  p  r  A9  Ar  \uz(z  +  Az)  -  uz(z)] 

-  p  A r  A z  [mo(9  +  A0)  -  mo(9)]  (B-7) 

-  p  A9  Az  [u,-(r  +  Ar)  r(r  +  Ar)  -  ur(r )  r(r)] 

Division  of  eq  B-7  by  the  quantity  (p  A9  Ar  Az)  and  taking  the  limit  as  A9, 
Ar,  and  Az  collapse  to  zero  gives  the  following  relationship: 

(dr/dt)  =  ~  r  (dujdz)  -  (6uo/dB)  -  [d(ur  r)/dr]  (B-8a) 

With  a  fixed  spatial  reference,  (dr/dt)  =  0,  and  eq  B-8a  simplifies  to 

r  (dujdz)  +  (due/dO)  +  [ d(ur  r)/dr\  =  0  (B-8b) 

or 

(dujdz)  +  (Hr)  (duJdO)  +  (dujdr)  +  (ujr)  =  0  (B-8c) 

The  relationship  given  by  eq  B-8c  is  eq  43. 
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APPENDIX  C 

ON  DIFFERENTIALS,  PARTIAL  DERIVATIVES, 
AND  TOTAL  DERIVATIVES 


The  differential  mass  balance  developments  in  this  effort  rely  heavily  on 
the  concepts  of  differentials,  partial  derivatives,  and  total  derivatives.  The  intent  of  this 
appendix  is  to  clarify  their  use  and  interpretation. 

Consider  the  molar  flux  which  has  a  functional  dependence  on  both  time 
and  spatial  position.  Time  and  the  positional  parameters  are  taken  to  be  independent 
variables.  For  a  cylindrical  coordinate  system,  the  functional  form  of  the  molar  flux  can 
be  represented  as  follows; 

V  =  V,  it,  r,  0,  z)  (C-l) 

The  molar  flux  corresponding  to  an  incremental  change  in  time  and 
position  can  be  represented  by  the  following  relationship: 

Ni  =  Ni  it  +  A t,r  +  A r,  9  +  A0,  z  +  Az)  (C-2) 

The  differential  or  total  change  of  molar  flux  in  making  this  time  and 
positional  transition  is  the  difference  between  eqs  C-2  and  C-l.  Therefore, 

A Ni  =  N/  it  +  A t,r  +  A r,  9  +  A0,  z  +  A z)  ~  Nt  (; t ,  r,  0,  z)  (C-3) 

The  relationship  given  by  eq  C-3  represents  the  total  differential  of  the 
molar  flux.  It  accounts  for  changes  in  the  flux  due  to  spatial  translations  in  each  vector 
direction  as  well  as  time. 


An  expression  that  relates  the  change  in  the  molar  flux  in  terms  of  its 
independent  variables  is  then  sought.  The  approach  is  detailed  as  follows. 

Addition  and  subtraction  of  V,-  (t,  r  +  Ar,  9  +  A0,  z  +  A z)  to  the  right  side  of 
eq  C-3  gives  the  following  expression: 


A/V,-  =  {  N,  {t  +  A t,r  +  Ar,  0  +  A9,  z  +  Az)  ~  Ni  it,  r  +  Ar,  9  +  A9,  z  +  Az) } 
+  Ni  it,  r  +  Ar,  0  +  A0,  z  +  Az)  -  Ni  it,  r,  0,  z) 


(C-4) 


Addition  and  subtraction  of  Ni  it,  r,  0  +  A9,  z  +  Az)  to  the  right  side  of 
eq  C-4  gives  the  following  expression: 
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AA,-  =  {  A  ( t  +  A t,r  +  A r,  9  +  A9,  z  +  A z)  ~  Nj  ( t ,  r  +  A r,  9  +  A9,  z  +  Az)  } 

+  {  A  ( t ,  r  +  Ar,  9  +  A9,  z  +  Az)  -  A  ( t ,  r,  9  +  A9,  z  +  Az) }  (C-5) 

+  Nj  ( t ,  r,  9  +  A9,  z  +  Az)  -  A  (t,  r,  9,  z) 

Addition  and  subtraction  of  A  (t,  r,  9,  z  +  Az)  to  the  right  side  of  eq  C-5 
gives  the  following  expression; 

A Nj  =  {Nj  ( t  +  At,  r  +  Ar,  9  +  A9,  z  +  Az)  -  A,-  ( t ,  r  +  Ar,  9  +  A9,  z  +  Az) } 

+  {Nj  ( t ,  r  +  Ar,  9  +  A9,  z  +  Az)  -  Nt  ( t ,  r,  9  +  A9,  z  +  Az) } 

(C-6) 

+  {  Nj  ( t ,  r,  9  +  A9,  z  +  Az)  -  Nt  ( t ,  r,  9,  z  +  Az) } 


+  {Nj  ( t ,  r,  9,  z  +  Az)  -  Nj  ( t ,  r,  9,  z) } 

Multiplying  the  first  bracketed  expression  on  the  right  side  of  eq  C-6  by 
(At/ At),  the  second  bracketed  expression  by  (Ar/Ar),  the  third  bracketed  expression  by 
(A9/A9),  and  the  fourth  bracketed  expression  by  (Az/Az),  gives  the  following  relationship: 

AA;-  =  { [Nj  ( t  +  At,  r  +  Ar,  9  +  A9,  z  +  Az)  -  A;  ( t ,  r  +  Ar,  9  +  A9,  z  +  Az)]/Af }  At 


+  { [A,-  (/,  r  +  Ar,  9  +  A9,  z  +  Az)  -  A,-  (/,  r,  9  +  A9,  z  +  Az)]/Ar  }  Ar 
+  { [A,-  (t,  r,  9  +  A9,  z  +  Az)  -  A,-  (t,  r,  9,  z  +  Az)]/A9  }  A9 


(C-7) 


+  { [A,-  (t,  r,  9,  z  +  Az)  -  A  ( t ,  r,  9,  z)]/Az  }  Az 

In  the  limit  as  the  differentials  of  the  independent  variables  and  the 
dependent  variable  approach  zero,  the  differentials  become  total  differentials,  and  the 
bracketed  expressions  become  partial  derivatives.  Performing  this  operation  on  eq  C-7 
gives  the  following: 

dA,-  =  ( dNj/dt )  dt  +  ( dNj/dr )  dr  +  (5A/39)  d9  +  ( dNj/dz )  dz  (C-8) 

Dividing  each  term  of  eq  C-8  by  the  total  time  differential  gives 
(dNj/dt)  =  (dNj/dt)  +  (dNj/dr)  (dr/d/)  +  (dNj/dQ)  (d9/df)  +  (dNj/dz)  (dz/dt)  (C-9a) 

The  (d/dt)  terms  represent  total  derivatives  with  respect  to  time.  As  given 
by  eq  C-9a,  they  correspond  to  directionally  oriented  velocities  associated  with  a  moving 
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reference  frame.  Following  the  path  of  a  specific  volume  element,  these  oriented 
velocities  can  be  represented  as  follows: 

vr  =  (dr/d/) 

ve  =  (d9/dt)  (C-lOa) 

vz  =  (dzJdt) 

Substituting  the  relationships  given  by  eq  C-10  into  eq  C-9a  gives  the 
following  relationship: 

(dNi/dt)  =  ( dNi/dt )  +  vr  ( dNj/dr )  +  ve  ( dNi/dd )  +  ( dNJdz )  (C-9b) 

The  relationship  given  by  eq  C-9b  would  be  applicable  if  one  were 
interested  in  the  change  in  molar  flux  as  a  given  fluid  element  traversed  the  adsorbent 
column.  In  many  situations  of  practical  importance,  one  of  the  directional  velocities  is 
taken  to  be  dominant  relative  to  the  others  and  it  is  the  principle  velocity  that  is 
considered.  The  gradients  and  velocities  associated  with  the  non-principle  directions  are 
then  assumed  to  be  negligible. 

For  instances  in  which  one  is  interested  in  the  molar  flux  relative  to  a  fixed 
spatial  position,  the  differential  volume  element  is  stationary  such  that  vr  =  v@  =  vz  =  0, 
and  eq  C-9b  gives  the  following  relationship: 

(dNi/dt)  =  (dNi/dt)  (C-ll) 

The  relationship  given  by  eq  C- 1 1  would  have  applicability  where  one  is 
interested  in  the  rate  of  change  of  the  molar  flux  with  time  with  respect  to  a  fixed  spatial 
reference. 
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APPENDIX  D 

ASPECTS  ASSOCIATED  WITH  MASS  TRANSFER  COEFFICIENTS 

IN  BINARY  SYSTEMS 


(a)  Development  of  Fundamental  Steady-State  Molar  Flux  Relationship  for  a 

Binary  System 

The  Dole  and  Klotz  developmental  efforts  introduce  a  mass  transfer 
coefficient  (kg)  and  a  logarithmic  mean  partial  pressure  of  the  nontransferring  component 
(pim,b)  in  their  analytical  relationships.  The  intent  of  this  appendix  is  to  provide  a  basis  for 
how  these  terms  arise.  The  effort  that  follows  is  based  on  a  treatment  given  in  a  text 
written  by  R.  E.  Treybal,  titled  Mass  Transfer  Operations  (ref  17). 

In  the  analysis  that  follows,  a  binary  fluid  system,  comprised  of  species  A 
and  B ,  is  considered.  The  transfer  of  mass,  relative  to  a  fixed  reference  plane,  will  be 
assumed  to  occur  in  a  singular  direction  only.  This  singular  direction  is  taken  to  be  the 
z-axis.  Figure  D-l  gives  a  diagram  of  the  system  that  was  used  in  this  analysis. 


Bulk  Fluid 
Phase 


- ►  Na 

Film 

Thickness 


Exterior  Surface  of 
Carbon  Particle 


z  =  0  z  =  5 

Ca\!C  ~  Cai/C 


ypa 


y  A2 


Figure  D-l.  Depiction  of  molar  flux  through  a  film  thickness. 


Through  Fick’s  law  of  diffusion,  the  molar  flux  in  the  z-di rection  (JA)  is 
proportional  to  a  diffusion  coefficient  (or  diffusivity),  and  its  concentration  gradient. 
Mathematically,  this  relationship  is  expressed  as  follows; 

JA  =  -Bab  (dCA/dz)  (D-l a) 

If  the  total  molar  concentration  (C)  is  taken  as  a  constant,  then  eq  D-l  a 
can  be  expressed  in  terms  of  its  mole  fraction  equivalent  as  follows: 

JA  =  -CBAB(8yA/dz )  (D-lb) 

The  net  total  molar  flux  of  this  binary  system  relative  to  a  fixed  reference 
plane  is  given  as 


N  =  Na  +  Nb 


(D-2) 


APPENDIX  D 


85 


The  transfer  of  A  in  this  binary  system  is  comprised  of  that  portion  of  the 
total  flux  which  is  comprised  of  component  A  and  that  resulting  from  its  diffusion. 
Therefore,  the  transfer  of  A  can  be  described  by  the  following  relationships: 


Na  —  N  yA  +  Ja 


=  (Na  +  Nb)  yA  ~  Bab  (dCA/dz.)  (D-3) 

=  (Na  +  Nb)  (Ca/C)  -  B Ab  (dCA/dz) 

By  analogy,  a  similar  expression  can  be  written  for  component  B  which 
yields  the  following; 

Nb  =  (. Na  +  NB)  (CB/Q  -  Bba  ( dCB/dz )  (D-4) 

Since  the  total  molar  concentration  is  given  by  C  =  (CA  +  CB ),  eq  D-3  and  eq  D-4  can  be 
added  together  which  gives  the  following  expression; 

-Bab  ( dCJdz )  =  Bba  ( dCB/dz. )  (D-5) 

If  the  total  concentration  can  be  taken  as  a  constant,  then  Dab  =  D BA.  This 
argument  demonstrates  the  equivalency  between  the  diffusivities  of  this  binary  system 
provided  that  the  total  molar  concentration  can  be  taken  as  a  constant. 

The  following  analysis  focuses  on  eq  D-3.  Separation  of  eq  D-3  yields  the 
following  relationship: 

~[C  NA  -  ( Na  +  NB)  CA]_1  8CA  =  (C  Dab)"1  dz  (D-6) 

With  the  total  molar  concentration  (C)  and  steady-state  localized 
component  fluxes  (i.e.,  NA  and  NB  as  constants),  eq  D-6  can  be  integrated  between  the 
concentration  limits  of  CA\  and  CA 2,  and  the  position  limits  z\  and  zi,  respectively.  The 
integration  is  carried  out  as  follows: 

Ca  2  Zi 

-/  [C  Na  ~  (Na  +  NB)  CA]~l  dCA  =  f  (C  D ab)"1  dz  (°-7) 

Cai  z  1 

With  C  and  Dab  as  constants,  the  integral  on  the  right  hand  side  of  eq  D-7  is  readily 
integrated  to  yield 

Z2 

he  Dab)'1  dz  =  (z2  -  Zi)/(C  Bab) 

Z 1 
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The  following  variable  substitutions  are  made  to  the  integral  on  the  left  side  of  eq  D-7: 


Letting  u  =  [C  NA  -  (N,\  +  NB)  Ca ],  the  following  differentials  result: 

(Gu/oCa)  =  [G(C  Na )/  GCA]  -  {  d[(NA  +  Nb)  Ca]/  GCa  } 

=  -(Na+Nb ) 
or 


dCA  =  ~(Na  +  Nb)  1  du 


Changing  the  limits  of  integration  (to  reflect  the  variable  substitution)  and 
making  the  variable  substitutions  into  the  integral  on  the  left  side  of  eq  D-7  yields  the 
following  expressions: 

Ca2  CNa-{Na  +  Nb)Ca2 

-/ [C  Na  ~  (i Na  +  NB)  CA]~l  dCA  =  ~  / ul  [-  (NA  +  NB)~l]  du 
Cai  CNa-(Na  +  Nb)Cai 

CNa-(Na  +  Nb)Ca2  (D-9) 

=  [1/(Na  +  Nb)]  fu~ldu 

CNA-(NA  +  NB)CA1 

=  [1/(Na  +  Nb)]  in  { [  C  NA  -  (NA  +  NB)  CAi]l\C  NA  -  (NA  +  NB)  CA i]  } 

Substituting  the  relationships  given  by  eqs  D-8  and  D-9  into  eq  D-7  yields 


[1/(Na  +  Nb)]  in  { [C  NA  -  (NA  +  NB)  CA2]/[C  NA  -  (NA  +  NB)  Cai]  } 

=  [(Z2-Zl)/(CDAB)] 

Equation  D-10  can  equivalently  be  expressed  as 


(D-10) 


[Na/(Na  +  Nb)]  in  {[CNA~  (NA  +  NB)  CA2]I[C NA  -  (NA  +  NB)  CA i]  } 


=  [Na(z2-za)/(CDab)] 


(D- 11) 


Letting  z  =  z2~Zi,  the  right  side  of  eq  D-l  1  can  be  solved  for  NA,  which  yields  the 
following  relationship: 


Na  =  [Na/(Na  +  Nb)]  [(C  D ab)/z\ 

x  In  {[CNa-  (Na  +  NB)  CA2]/[C NA  -  (NA  +  NB)  CAl]  } 


(D-12) 
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The  argument  of  the  logarithm  in  eq  D-12  can  equivalently  be  expressed 
through  the  following  relationships: 


[C  na  -  C na  +  Nb )  cA2]/[c  Na  -  (Na  +  Nb )  CAl] 


=  (Na  +  Nb)  { [(C  NA)/(NA  +  NB)\  -  C,v2}/{(Mi  +  NB)  { [(C  NA)/(NA  +  NB)]  -CA i }  } 
=  C  {  [Na/  (Na  +  Nb)]  ~  (CA2/C)  }/{c{  [Na/  (Na  +  NB)]  -  (CAl/C) }  } 

=  {  [Na/  (Na  +  Nb)]  -  (CA2/C) }/{ [NA/  (NA  +  NB)]  -  (CAl/C) } 

Substituting  the  relationship  given  by  eq  D-13  into  eq  D-12  yields 
Na  =  [Na/(Na  +  Nb)]  [(C  D ab)/z] 

x  in  {  {  [Na/  (Na  +  NB)]  -  (CA2/C) }  /  { [NJ  (NA  +  NB)]  -  (CAl/C) }  } 


(D-13) 


(D-14a) 


With  the  quantity  [(C  D AB)lz\  representing  a  mass  transfer  coefficient  such 
as  F,  this  last  expression  can  be  expressed  as 


Na  =  [Na/(Na  +  Nb)]F 


x  in  {  {  [Na/  (Na  +  NB)]  -  (CA2/C) }  /  { [NA/  (NA  +  NB)]  -  (CA  [/C) }  } 


(D-14b) 


The  development  of  eq  D-14b  was  based  on  the  assumptions  of  a  constant 
total  molar  concentration  and  steady-state  molar  fluxes  for  each  of  the  mobile  species. 


(b)  Background  Information  on  Analytical  Basis  of  Mass  Transfer 

Coefficients 

Some  background  information  on  the  use  and  concepts  associated  with 
mass  transfer  coefficients  is  called  for  at  this  stage.  Much  of  the  information  that  follows 
is  extracted  from  Mass  Transfer  Operations ,  by  Robert  E.  Treybal  (ref  17). 

Molar  fluxes  are  generally  written  in  the  following  format: 

Na  =  (mass  transfer  coefficient)  (driving  force)  (D-15a) 

Combinations  of  mass  transfer  coefficients  and  driving  forces  can  be  expressed  as 
follows. 
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Table  D-l.  Mass  Transfer  Coefficients  and  Corresponding  Gradients 
_ (variables  used  for  NB  =  0;  see  Table  D-2) _ 


Driving  Force  Basis 
(Gradient) 

Gas  Phase 

Mass  Transfer  Coefficient 

Pressure  (i.e.,pAi  ~Pai ) 

kG 

Mole  fraction  (i.e.,  yAi  -  yAi) 

Ity 

Concentration  (i.e.,  CA i  -  CA 2) 

kc 

Mass  ratio  (mass  A/mass  B) 

ky 

The  units  on  the  k-type  mass  transfer  coefficients  must  be  consistent  with 
the  units  of  the  driving  force  gradient.  While  the  molar  flux  of  a  component  (e.g.,  NA)  is 
assigned  with  the  units  of  [moles  transferred/(area  time)],  dimensional  consistency 
requires  some  consideration  on  the  specification  of  the  units  for  the  k- type  mass  transfer 
coefficient.  Table  D-2  is  based  on  an  extraction  from  Treybal  (ref  17;  p  49). 


Table  D-2.  Relations  between  Mass  Transfer  Coefficients  (ref  17) 


Equimolar 
Counterdiffusion 
(Na  =  -Nb) 

Diffusion  of  A 
through 
Nondiffusing  B 

C Nb  =  0) 

Units  of  Mass  Transfer  Coefficient 

Gases 

Gases 

Gases 

Na  =  kG'  A PA 

Na  =  kG  A  Pa 

Moles  transferred/[(area)  (time)  (pressure)] 

Na  =  ky'  Ay A 

Na  =  ky  Ay  a 

Moles  transferred/[(area)  (time)  (mole  fraction)] 

3 

<1 

11 

£ 

Na  =  kc  Ac A 

Moles  transferred/[(area)  (time)  (mol/vol)] 

WA  =  ky  AY  A 

Mass  transferred/[(area)  (time)  (mass  fraction)] 

With  this  terminology  approach,  eq  D-15a  can  be  expressed  in  the 
following  analytical  form  for  a  binary  gas-phase  mass  transfer  system  in  which 
component  A  is  diffusing  through  a  mass  transfer  film  thickness  and  component  B  is  not 
diffusing  through  the  mass  transfer  film  thickness: 


Na  =  kG  A pA  =  kG(pA\- Pai)  (D- 1 5b) 

where  A pA  =  (pai  ~  Pai)  is  a  pressure -based  driving  force,  and  kG  is  a  gas-phase  mass 
transfer  coefficient. 

Mass  transfer  coefficients  are  generally  correlated  through  a  Reynolds 
(Re)  number  and  a  Schmidt  (Sc)  number  through  analytical  fits  of  data.  The  correlating 
parameter  that  contains  the  mass  transfer  coefficient  is  known  as  the  Sherwood  (Sh) 
number.  Such  correlations  generally  take  on  the  following  form 

Sh=/(Re,  Sc)  (D-l 6) 
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The  parameters  used  in  the  Sherwood  number  depend  on  the  basis  of  the 
mass  transfer  coefficients  and  the  driving  force.  For  instance,  the  mass  transfer 
coefficient  (F)  as  depicted  in  eq  D-14b  can  be  shown  to  be  related  to  the  gas-phase  mass 
transfer  coefficient  in  eq  D-15b  (for  gas-phase  mass  transfer  in  which  there  is  the 
diffusion  of  component  A  through  nondiffusing  component  B)  through  the  following 
expression: 


F  =  kGp\m,B  (D-17) 

where  pimp  is  the  logarithmic  mean  partial  pressure  of  component  B.  Treybal  (ref  17) 
gives  the  following  correspondences  among  the  transfer  coefficients  for  binary  diffusion: 


F  =  kc,  p\m,B  =  ky  (pun.li  IP,)  =  kc  (p\m,B  IRT) 

=  kY/MB  =  kG  p,  =  ky  =  kc'  ip, IRT)  =  kc'  c 


(D-18) 


(c)  Ideal  Physical  Adsorption  Considered  as  the  Transfer  of  A  through 

Nontransferring  B 

Two  situations  are  now  considered  for  the  transfer  of  mass  through  the 
film.  These  are  described  as  follows: 

1 .  Equimolar  Counter  Diffusion:  In  this  situation,  the  net  molar  flux  in  the 
z-direction  is  zero  such  that  NA  =  ~NB.  By  the  relationship  given  by 

eq  B-14b,  this  would  require  that  [NA/(NA  +  NB)]  approach  infinity. 

2.  Transfer  of  A  through  Nontransferring  B:  In  this  situation,  component  B 
is  essentially  motionless  such  that  NB  =  0.  Through  the  relationship 
given  by  eq  B-14b,  this  would  require  that  [NA/(NA  +  NB)]  =  1. 

Assuming  the  process  of  ideal  physical  adsorption  to  be  one  in  which  there 
is  the  transfer  of  a  unitary  component  (i.e.,  the  vapor-phase  contaminant)  through  a 
medium  largely  composed  of  a  nontransferring  component  (i.e.,  air),  then  the  situation 
becomes  one  in  which  NB  =  0  and  eq  D-14b  simplifies  to  the  following  relationship: 

Na  =  F  ln{  [1  -  (CA2/C)]/[1  -  (CA,/C)J }  (D-14c) 

The  concentration  ratios  appearing  in  eq  D-14c  are  mole  fractions  and  are 
related  to  pressures  through  the  following  expressions 

yA,=  iCA\IC)  =  IpmHRT)]  /  \pl(RT)\  =pM/p 

and 

V/12=  (Cai/C)  =  [ paiURT )]  /  \p/(RT)\  =  pA2/p 
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Substituting  these  relationships  into  eq  D-14c  gives 

Na  =  F  ln{  [1  -  (pA2/p)V[  1  -  (pai/p)]  }  (D-14d) 

Development  of  the  p\m,b  factor  that  appears  in  the  Dole  and  Klotz 
reference  (ref  3)  and  is  given  by  eq  201  is  based  on  the  following  approach. 

The  molar  flux  relationships  as  given  by  eqs  D-14d  and  D-15b  are 
equated,  which  gives  the  following  expression: 

Na  =  F  ln{  [1  -  (pA2/p)]/[l  -  ( Pai/p )]  }  =kG  (pai  ~Pai)  (D-19) 


Therefore, 


F  =  kG(pA  1  -  PA2)/  In  { [  1  -  (pai/pM  1  -  (pai/p)]  }  (D-20a) 

For  a  binary  system,  the  following  partial  pressure  relationships  would 

hold: 


and 


(Pa\-Pa2)  =  (Pb2~Pb\) 

Pai  =P  ~Pb\ 

PA2=P  ~PB2 

Substituting  these  relationships  into  eq  D-20a  gives 

F  =  kG  (pH1  -  pm)/  \n[p,rJpm\  (D-20b) 

The  definition  of  the  logarithmic  mean  partial  pressure  of  component  B  is 


given  as 


/I I m, W  =  (P 112-  Pill)/  MP/pJPlU  I 


(D-21) 


With  the  definition  of  the  logarithmic  mean  partial  pressure  now  defined, 
eq  D-20b  can  be  expressed  as 


F  =  k(j  p\ 


m,/j 


(D-20c) 


Equation  D-20c  is  consistent  with  the  interpretation  of  F  as  given  by  eqs  D-17  and  D-18. 
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